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The  material  in  this  dissertation  can  be  effectively 
divided  into  two  subtopics:   philosoohy  of  optimal  mechanism 
design,  and  optimization  of  dyad-based  spatial  mechanisms. 

The  first  subtopic,  philosophy  of  optimal  mechanism 
design,  is  intended  to  be  general  in  nature,  applying  to  all 
types  of  mechanisms,  both  higher  and  lower  pair,  and  both 
planar  and  spatial.   This  is  covered  in  Chapters  One  through 
Three.   Chapter  One  examines  past  approaches  to  mechanism 
optimization.   Chapter  Two  is  a  brief  review  of  optimization 
theory,  particularly  as  it  applies  to  mechanism  optimization. 
Chapter  Three  draws  upon  the  insights  gained  in  the  first 
two  chapters  to  formulate  a  general  approach  to  the  mechanism 
optimization  problem. 

The  second  subtopic  of  this  dissertation,  optimization 
of  dyad-based  spatial  mechanisms,  is  covered  in  Chapters 


Four  through  Seven.   This  is  actually  a  rather  limited 
example  of  applying  the  philosophy  developed  in  the  first 
three  chapters.   Nevertheless,  the  mechanisms  treated  in 
this  section  are  believed  to  represent  some  of  the  most 
useful  motion  generating  spatial  mechanisms,  and,  therefore, 
those  for  which  improved  design  theories  are  most  urgently 
needed.   In  Chapter  Four,  closed-form  synthesis  equations 
are  derived  for  dyads  containing  revolute  (R) ,  spheric  (S) 
and  cylindric  (C)  pairs.   Chapters  Five  and  Six  present 
detailed  examples  of  the  optimization  of  the  four-link 
RCCC  and  five-link  RSSR-SC  and  RSSR-SS  mechanisms.   Finally, 
Chapter  Seven  outlines  procedures  for  the  optimization  of 
other  dyad-based  spatial  mechanisms,  and  offers  suggestions 
for  further  research. 


CHAPTER  1 
BACKGROUND  AND  MOTIVATION  FOR  THIS  RESEARCH 

1 . 1   Introduction 

The  design  of  mechanisms,  like  almost  all  design 
problems,  is  an  iterative  process.   Generally,  a 
mechanism  can  be  synthesized  to  meet  some  of  the  design 
requirements,  but  then  it  must  be  analyzed  or  tested  to 
determine  whether  or  not  it  satisfies  the  remaining 
requirements.   Most  often  it  does  not,  and  the  designer 
must  "go  back  to  the  drawing  board"  and  synthesize  a 
new  mechanism  which  again  must  be  analyzed  or  tested. 
This  second  design  may  also  not  satisfy  all  the  design 
requirements,  but  it  is  likely  to  be  somewhat  of  an 
improvement  over  the  first  design,  since  the  designer 
had  the  "experience"  of  the  first  design  on  which  to 
base  this  second  synthesis.   The  third  design  attempt 
will  likely  be  better  still,  since  now  two  points  of 
reference  are  available.   This  iterative  process  is 
continued  until  a  satisfactory  solution  is  obtained. 
It  is  probably  the  oldest  and  most  often  used  form  of 
ootimization. 


For  the  purpose  of  this  dissertation,  optimization 
may  be  considered  to  be  the  process  of  seeking  the  best 
result  under  a  given  set  of  circumstances.   The  questions 
of  what  optimization  processes  are  available,  what  is 
meant  by  the  "best"  result  and  how  these  "circumstances" 
are  incorporated  into  the  solution  will  be  topics  of 
considerable  importance.   The  following  paragraphs 
briefly  introduce  these  concepts. 

As  has  already  been  stated,  optimization  may  be  any 
process  which  seeks  to  find  the  best  result  to  a  given 
problem.   Available  optimization  methods  include  the 
intuitive  successive  improvement  approach  described 
previously;  classical  optimization  methods,  such  as  the 
Lagrange  multiplier  method,  based  on  the  principles  of 
variational  calculus;  and  a  large  and  rapidly  expanding 
body  of  knowledge  known  as  mathematical  programming. 
The  last  of  these,  mathematical  programming,  is  most 
important  with  respect  to  this  work.   In  a  sense, 
mathematical  programming  methods  are  an  extension  and 
formalization  of  the  age-old  successive  improvement 
process.   These  newer  methods,  though,  have  become  quite 
sophisticated,  employing  a  variety  of  numerical  techniques 
and  almost  always  requiring   programming  on  a  digital 
computer.   Also,  by  formulating  some  measurement  of  the 
quality  of  a  given  design,  the  burden  of  decision- 
making after  each  iteration  can  be  removed  from  the 


designer.   The  most  general  form  of  mathematical 
programming  aims  at  solving  nonlinear  types  of  problems 
and,  not  surprisingly,  is  called  nonlinear  programming. 
It  will  soon  become  evident  that  mechanism  optimization 
generally  involves  nonlinear  programming.   A  detailed 
discussion  of  optimization  methods  is  given  in  Chapter 
Two. 

The  second  question  to  be  addressed  is  what  is  meant 
by  the  "best"  result.   Suppose,  for  example,  that  an 
airplane  strut  is  to  be  designed  for  a  given  design  load 
with  maximum  stiffness,  but  with  minimum  weight.   It  is 
apparent  that  these  represent  competing  design  require- 
ments, and  that  some  compromise  between  the  two  must  be 
reached.   Here,  the  designer's  judgment  must  enter  the 
process  and  he  must  decide  what  importance  or  "weight" 
to  assign  to  each  requirement.   Of  course,  mechanism 
design  problems  are  rarely  this  simple,  and  will  usually 
involve  a  much  larger  number  of  design  requirements.   As 
a  result,  finding  the  best  design  usually  involves 
some  amount  of  subjective  opinion.   Even  if  the  importance 
of  the  various  design  requirements  were  exactly  known, 
the  best  mathematical  solution  would  probably  not 
correspond  to  the  best  physical  solution.   This  is  because 
there  are  inherent  errors  and  approximations  present  in 
the  modeling  and  analysis  of  any  real  system.   Thus, 
although  mathematicians  may  speak  of  a  "globally  optimum" 


or  best  solution,  the  best  that  can  be  hoped  for  in  an 
engineering  setting  is  to  find  some  approximation  to  the 
best  solution. 

The  final  topic  of  this  section  addresses  the 
question  of  what  role  the  "circumstances"  play  in 
determining  a  best  design.   Will  the  machine  designed 
to  run  at  60  revolutions  per  minute  (rpm)  perform 
acceptably  at  600  rpm?   Perhaps  not,  since  the  circum- 
stances or  "constraints"  under  which  the  design  was 
developed  have  been  altered.   Constraints  may  be  cast  in 
the  form  of  equality  or  inequality  conditions,  and  they 
may  be  explicit  or  implicit.   In  any  case,  it  is  evident 
that  the  circumstances  surrounding  a  design  will  also 
have  an  effect  on  which  solution  is  best. 

1. 2   The  Elements  of  Mechanism  Optimization 

Generally  speaking,  the  optimization  of  any  system 
or  component  comprises  four  principal  elements.   These  are 

(1)  the  choice  of  a  conceptual  design; 

(2)  the  development  of  an  analyzable  model 
of  the  physical  system  from  which  the 
design  variables  and  the  governing 
equations  may  be  extracted,  both  for 
synthesis  and  analysis  ; 

(3)  the  setting  up  of  a  scalar  function 
(called  the  objective  function)  of  the 
design  variables  which  serves  as  a 


measure  of  the  effectiveness  of  the 
system,  accompanied  by  specifying 
all  design  constraints;  and  finally, 
(4)   finding  the  set  of  values  of  the  design 
variables  which  produce  the  best  value 
of  the  objective  function  of  part  (3) 
and  are  consistent  with  all  design 
constraints . 

The  first  three  of  these  elements  are  primarily 
dependent  on  the  specific  problem  at  hand;  they  will 
collectively  be  called  the  problem  formulation.   The 
fourth  element  is  purely  an  application  of  optimization 
theory.   One  of  the  primary  objectives  of  this  dissertation 
is  to  develop  a  method  of  mechanism  problem  formulation 
in  such  a  way  as  to  accurately  reflect  the  problem 
requirements,  while  maintaining  a  simple  enough  structure 
to  allow  efficient  and  sufficiently  accurate  solutions 
by  one  of  the  available  optimization  methods. 

The  first  element  of  the  problem  formulation  to  be 
addressed  is  choosing  a  conceptual  design.   The  designer 
must  decide  what  types  of  mechanisms,  such  as  linkages, 
cams,  belts,  gears,  cam-modulated  linkages,  etc.,  are 
best  suited  to  solve  a  particular  problem.   Following 
this,  he  must  decide  the  particular  configuration  of  the 


mechanism  chosen.   If  a  linkage  has  been  selected,  should 
it  have  four  links  or  six  links?   Should  it  be  planar  or 
spatial?   The  above  decisions  are  known  as  type  and 
number  synthesis,  respectively.   Although  these  are 
extremely  important  decisions  in  the  optimization  process, 
a  detailed  discussion  of  them  is  beyond  the  scope  of  this 
work.   Throughout  the  remainder  of  this  dissertation,  it 
will  be  assumed  that  a  specific  mechanism  type  has  been 
selected,  and  optimization  must  proceed  from  this  point. 

The  second  element  of  problem  formulation  consists 
of  synthesis,  modeling  and  analysis  of  the  system  concep- 
tualized above.   In  the  present  work,  this  stage  will 
draw  heavily  upon  the  classical  theories  of  mechanism 
synthesis  and  analysis.   Here  the  term  analysis  includes, 
among  other  things,  application  of  Grashof's  condition 
(or  a  similar  spatial  condition) ,  branch  and  order  deter- 
mination and,  of  course,  position,  velocity,  and 
acceleration  analysis.   Typically,  the  mechanism  to  be 
optimized  will  be  modeled  as  having  rigid  links  and  zero 
clearance  in  the  joints,  although  other  assumptions  such 
as  elastic  links  and  joint  clearances  and/or  compliances 
may  be  applied  where  appropriate. 

The  third  and  final  element  of  the  problem  formulation 
is  the  setting  up  of  the  objective  function  and  the  design 
constraints.   The  objective  function  is  a  scalar  function 
of  the  design  variables  whose  numerical  value  reflects 


the  quality  of  the  system  being  designed.   In  designing 
a  mechanism  for  a  specified  performance,  it  may  be 
desired  to  minimize  the  sum  of  the  squares  of  the 
structural  errors   at  a  large  number  of  positions.   Or 
perhaps  it  is  desired  to  minimize  the  maximum  deviation 
of  the  transmission  angle  from  the  ideal  value  of  90 
throughout  the  motion  cycle.   Maybe  some  weighted 
combination  of  these  two  requirements  in  the  objective 
function  is  best.   It  is  evident  that  the  designer  must 
possess  a  thorough  understanding  of  the  requirements 
of  the  system  being  designed  if  he  is  to  choose  an 
effective  objective  function. 

1. 3   History  and  Literature  Review 

This  dissertation  brings  together  two  fields,  namely, 
optimization  theory  and  mechanism  science,  which  histor- 
ically have  a  related  and  remarkably  similar  development. 
Both  disciplines  have  their  roots  in  early  civilization, 
although  neither  was  formally  recognized  until  quite 
recently.   Both  disciplines  were  topics  of  interest  and 
were  significantly  advanced  by  the  great  mathematicians 
of  the  seventeenth  and  eighteenth  centuries.   Newton, 
Euler  and  Lagrange  are  but  a  few  of  the  men  whose  names 


Inherent  deviation  from  the  prescribed  performance, 
expressed  as  a  scalar  or  vector  quantity. 


are  common  to  both  subjects.   Even  with  the  contributions 
of  these  and  many  other  great  men,  development  of  both 
subjects  remained  slow  and  sporadic  well  into  the 
twentieth  century. 

By  about  1945,  both  mechanism  science  and  optimization 
theory  had  reached  somewhat  of  a  plateau.   Graphical 
methods  of  planar  mechanism  design  and  analysis  were  well 
established,  although  their  usefulness  was  limited  by  the 
sometimes  inaccurate  and  tedious  nature  of  graphical 
constructions.   Classical  optimization  methods,  based  on 
calculus,  were  also  well  established,  but  were  useful  for 
solving  only  a  limited  class  of  problems.   Both  fields, 
mechanism  science  and  optimization  theory,  were  advancing 
at  a  relatively  slow  pace.   The  possibility  of  solving 
new  types  of  problems  or  more  difficult  problems  seemed 
distant  indeed. 

Then  came  the  dawning  of  the  computer  age.   Suddenly, 
the  potential  arose  for  handling  large  amounts  of  informa- 
tion, and  for  performing  tedious  calculations  quickly  and 
accurately.   The  growth  rate  of  both  fields  virtually 
exploded.   Kinematicians  adopted  algebraic  methods  which 
were  soon  translated  into  computer  code.   It  became 
possible  to  generate  hundreds  or  even  thousands  of 
mechanisms  as  possible  candidates  to  solve  a  given  problem. 
It  soon  became  apparent  that  the  computer  could  generate 
more  solutions  than  the  designer  could  reasonably  assess. 


The  next  logical  step  was  to  let  the  computer  examine 
the  relative  merits  of  each  solution,  presenting  to 
the  designer  only  those  mechanisms  which  met  some 
preestablished  criteria.   It  was  quickly  recognized 
that  the  rapidly  advancing  field  of  optimization  theory 
was  exactly  the  tool  needed  for  this  purpose.   The 
premise  was  a  simple  one:   to  allow  the  computer  to 
perform  a  logical  search  for  the  mechanism  or  mecha- 
nisms which  best  satisfied  the  designer's  requirements 
under  the  given  circumstances.   This  is  the  foundation 
of  practically  all  modern  mechanism  optimization. 

The  previous  section  gave  a  very  brief  summary  of 
the  early  history  of  kinematics  and  optimization. 
This  section  outlines  the  developments  in  kinematics 
since  about  1950,  particularly  as  they  relate  to 
synthesis  and  optimization  of  spatial  mechanisms. 
Many  of  the  references  cited  here  deal  exclusively 
with  planar  mechanisms.   This  is  to  be  expected  since 
often,  although  not  always,  the  development  of  planar 
methods  of  synthesis  and  optimization  have  been,  or 
can  be,  extended  directly  to  spatial  mechanisms. 
Also,  the  inclusion  of  these  planar  references  will 
help  to  explain  how  the  trends  toward  many  of  the 
currently  used  methods  were  established. 
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Optimization  theory,  as  a  subject  unto  itself,  is 
not  directly  treated  in  this  literature  review.   Rather, 
in  the  later  section  on  optimization  methods,  a  number 
of  standard  references  are  cited. 

The  pioneer  of  modern  kinematics  in  the  United  States 
is  generally  considered  to  be  Freudenstein.   His  1955 
paper,  "Approximate  Synthesis  of  Four-Bar  Linkages,"  (1) 
probably  marks  the  beginning  of  the  shift  in  emphasis 
from  graphical  to  analytical  methods.   The  expression, 
"Approximate  Synthesis,"  was  widely  used  during  this  era 
to  denote  precision  position  synthesis  to  approximate  a 
given  function.   It  should  not  be  confused  with  its  more 
recent  use  to  denote  synthesis  of  a  mechanism  which 
approximately  satisfies  a  large  number  of  precision 
positions.   Freudenstein ' s  later  work  with  Sandor  (2)  in 
1959  employed  complex  number  theory  and  a  programmed  IBM 
650  digital  computer  in  the  synthesis  of  path-generating 
mechanisms.   This  work,  perhaps  more  than  any  other,  marks 
the  marriage  of  the  kinematician  to  the  digital  computer. 

Other  important  contributions  from  this  time  period 
dealing  with  planar  mechanisms  include  the  four-bar 
linkage  atlas  of  Hrones  and  Nelson  (3) ,  the  additional 
works  of  Freudenstein  (4,5)  and  the  combined  work  of  Roth 
and  Freudenstein  (6) ,  to  name  a  few.   In  the  latter  work, 
path  generating  geared  five-bar  linkages  are  synthesized 
using  numerical  methods  and  a  digital  computer.   The  work 
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of  A.S.  Hall  (7),  along  with  his  organizational  involve- 
ment in  the  early  "Conferences  on  Mechanisms"  from  1953 
to  1962,  gave  inspiration  and  insight  to  many  of  the 
researchers  who  followed  him.   McLarnan  (8)  extended 
the  earlier  work  of  Freudenstein  to  include  the  synthesis 
of  planar  six-link  function-generating  mechanisms. 

Analytical  studies  of  spatial  mechanism  synthesis 
in  the  United  States  also  received  attention  during  the 
1950 '  s.   Denavit  and  Hartenberg  (9)  provided  what  has 
become  standard  symbolic  notation  for  the  description 
of  the  kinematic  properties  of  lower-pair  mechanisms. 
In  a  later  paper  (10) ,  Denavit  and  Hartenberg  extended 
the  precision  position  approach  of  Freudenstein  (1)  to 
the  synthesis  of  spatial  RSSR  and  RCCC  mechanisms,  and 
showed  the  synthesis  equations  to  be  linear  up  to  a 
limited  number  of  precision  positions. 

It  should  be  noted  that  a  number  of  German  and 
Russian  scholars  were  also  contributing  to  the  advancement 
of  knowledge  on  spatial  mechanisms  at  this  time.   Notable 
among  these  are  Eeyer  (11) ,  Dimentberg  (12) ,  Novodvorskii 

(13) ,  Stepanov  (14) ,  and  Levitskii  and  Shakvazian  (15) . 
For  a  more  detailed  review  of  these  works,  the  reader  is 
referred  to  the  survey  articles  of  Beyer  (16) ,  Harrisberger 

(17)  and  Yang  (18)  . 

Until  about  1965,  the  only  spatial  mechanism 

synthesis  problems  covered  in  the  literature  involved 


coordinating  motions  of  input  and  output  links,  or  as 
it  is  commonly  known,  the  function  generation  problem. 
Wilson  (19)  changed  this  trend  by  introducing  the 
problem  of  spatial  rigid-body  guidance.   He  also  showed 
the  function  generation  problem  to  be  convertible  to  a 
rigid-body  guidance  problem  by  inversion  about  the 
input  or  output  links. 

The  advancing  computer  technology  of  the  late  1950 's 
gave  rise  to  the  first  numerically-based  attempts  to 
optimize  planar  mechanisms.   Freudenstein  (20)  devised 
an  iterative  scheme  for  respacing  the  precision  points 
in  order  to  find  the  linkage  which  best  approximates  a 
given  function.   Starting  with  the  well-known  Chebychev 
spacing  of  the  precision  points,  a  mechanism  was 
synthesized.   An  analysis  routine  then  evaluated  the 
resulting  structural  error  over  the  range  of  the  desired 
function.   Following  this,  a  new  "modified  Chebychev" 
spacing  of  the  precision  points  was  created,  and  the 
analysis  step  repeated.   The  resulting  errors  of  the  two 
trials  were  then  compared,  and,  based  on  the  result,  a 
new  spacing  could  be  projected. 

Roth,  Freudenstein  and  Sandor  (21)  used  a  digital 
computer  to  synthesize  planar  four-bar  path  generating 
linkages  with  optimum  transmission  characteristics.   This 
was  done  by  synthesizing  linkages  to  satisfy  four  precision 
conditions  in  addition  to  requiring  the  maximum  and  minimum 


values  of  the  transmission  angle  to  deviate  from  90   by 
the  same  numerical  amount.   Synthesized  mechanisms  were 
then  compared  by  means  of  a  "quality  index"  formula  and 
iterations  were  performed  until  the  mechanism  with  the 
best  value  of  the  index  could  be  found. 

Until  the  work  of  Chi-Yeh  (22)  was  published  in 
1966,  synthesis  methods  invariably  required  the  exact 
satisfaction  of  some  number  (usually  3,4  or  5)  of 
precision  conditions.   Chi-Yeh 's  novel  approach  was  to 
minimize  the  structural  error  in  the  least-squares  sense 
at  a  much  larger  set  of  prescribed  positions.   Chi-Yeh' s 
method  involved  taking  the  partial  derivatives  of  the 
objective  function  with  respect  to  the  design  variables 
(the  linkage  dimensions) ,  and  setting  these  equal  to 
zero.   This  resulted  in  a  system  of  nonlinear  equations 
which  were  then  solved  numerically. 

Nonlinear  programming  was  first  introduced  to 
mechanism  optimization  by  Fox  and  Willmert  (23)  in  1967. 
Their  design  objective  was  to  synthesize  a  four-bar 
linkage  whose  coupler  point  would  generate,  as  closely 
as  possible,  a  given  curve,  and  whose  crank  rotations 
would  be  as  close  as  possible  to  the  desired  values. 
Constraints  were  imposed  which  limited  forces  and  torques 
within  the  linkage,  restricted  fixed  pivot  location, 
limited  link  length  ratios,  assured  the  desired  Grashof 
type  and  required  the  output  positions  to  be  in  a 
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specified  order.   Also,  the  possibility  of  limiting 
velocities  and  accelerations  below  a  certain  value  was 
discussed.   This  paper  is  remarkable  in  that  it  introduced 
nonlinear  programming  to  kinematic  synthesis  in  addition 
to  developing  a  number  of  the  constraint  conditions  which 
are  still  in  common  use  today.   In  fact,  many  of  the  more 
recent  works  in  this  area  have  been  in  an  effort  to  find 
more  efficient  optimization  algorithms;  however,  the 
objective  function  and  constraint  equations  have  been 
virtually  the  same  as  those  suggested  by  Fox  and  Willmert. 
Freudenstein,  in  his  discussion  of  the  paper  (23), 
described  the  nonlinear  programming  method  as  a  "natural 
tool"  in  the  area  of  kinematic  synthesis. 

Closely  following,  but  apparently  independent  from, 
the  work  of  Fox  and  Willmert,  was  a  paper  by  Tomas  (24) 
which  also  treated  linkage  synthesis  as  a  nonlinear 
optimization  problem.   Noteworthy  about  the  work  of  Tomas 
was  the  first  treatment  of  the  function  generation  problem 
by  nonlinear  programming,  and  the  first  published  example 
of  optimization  with  respect  to  dynamic  properties  of 
the  mechanism. 

A  different  tack  to  the  planar  linkage  optimization 
problem  was  presented  by  Garrett  and  Hall  (25) .   They 
chose  to  generate,  by  means  of  a  digital  computer,  a  set 
of  random  four-bar  linkages.   Each  of  these  linkages  was 
then  analyzed  to  find  the  one  which  best  suited  the 
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design  requirements.   Following  this,  an  expanded  set 
of  random  linkages  was  generated  in  the  neighborhood 
of  this  "best"  design,  and  repeating  the  analysis  step, 
a  new  best  linkage  was  selected.   This  refinement  process 
was  repeated  until  the  desired  accuracy  was  obtained,  or 
until  the  method  converged.   Although  it  is  usually 
inefficient,  this  method  is  quite  simple  and  is  almost 
guaranteed  to  converge.   Furthermore,  it  may  have  a 
better  chance  of  finding  the  global  optimum  (26), 
although  no  proof  of  this  is  known  to  exist. 

Eschenbach  and  Tesar  (27)  contributed  to  the  theory 
of  planar  linkage  optimization  on  two  important  fronts. 
First,  they  treated  the  formerly  unsolved  problem  of 
optimum  design  of  a  coplanar-motion  generating  four-bar 
linkage.   Second,  they  were  the  first  to  recognize  the 
simplification  which  resulted  from  using  precision 
position  equality  constraints  to  reduce  the  number  of 
free-choice  variables. 

Least-squares  error  synthesis  of  planar  four-link 
mechanisms  was  a  popular  subject  during  the  1960 's,  and 
attracted  the  interest  of  such  authors  as  Lewis  and  Gyory 
(28) ,  Levitskii  and  Sarkisian  (29)  and  McLarnan  (30) 
among  others.   Although  none  of  these  papers  directly 
employ  nonlinear  programming  methods,  they  do  point  out 
the  importance  given  to  the  problem  of  minimizing 
structural  error.   Among  these  works,  only  the  paper  of 
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Lewis  and  Gyory  (28)  considers  the  need  to  satisfy 
additional  constraints.   Here,  a  test  is  performed  to 
be  sure  the  synthesized  mechanism  is  of  the  crank-and- 
rocker  type.   Structural  error  minimization  was  also 
the  objective  of  Sandor  and  Wilt  (31)  in  their  novel 
paper  on  the  optimal  synthesis  of  a  geared  four-link 
mechanism. 

Tomas  (32)  discounts  the  importance  of  minimizing 
structural  error  in  many  practical  problems,  and  again 
presents  the  idea  of  optimization  with  respect  to  the 
dynamic  properties  of  the  mechanism.   Benedict  and  Tesar 
(33)  reinforce  this  concept  by  demonstrating  the  design 
of  a  complex  stamping  and  indexing  machine  with  optimal 
torque  balance.   Improved  dynamic  properties  and 
balancing  of  planar  mechanisms  were  also  the  objective 
of  a  number  of  other  papers,  including  those  by  Berkof 
and  Lowen  (34,35),  Huang,  Sebesta  and  Soni  (36),  Sadler 
and  Mayne  (37)  and  Elliot  and  Tesar  (38) .   Kaufman  and 
Sandor  (39)  developed  a  complete  force  balancing  method 
for  spatial  mechanisms. 

Variations  of  the  planar  mechanism  optimization 
problems  already  mentioned  have  been  the  subject  of 
numerous  other  papers.   Most  of  these  papers  attempt  to 
improve  the  efficiency  of  the  optimization  process, 
either  by  using  slightly  different  problem  formulations, 
or  by  using  one  of  the  recently  improved  optimization 
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methods.   These  include  the  works  of  Golinski  (40) , 
Alizade,  Novruzbekov  and  Sandor  (41),  Rose  and  Sandor 
(42),  Savage  and  Suchora  (43),  Nolle  and  Hunt  (44), 
Bagci  and  Brosfield  (45)  and  Bagci  and  Lee  (46). 

As  work  in  the  area  of  mechanism  optimization 
progressed,  the  importance  of  good  problem  formulation 
combined  with  an  efficient  numerical  optimization  method 
became  apparent.   Kinematicians  began  to  modify  existing 
optimization  methods,  or,  in  some  cases,  create  new 
methods,  which  would  be  tailored  to  the  peculiarities 
of  mechanism  design.   Huang  (47)  developed  a  mechanism 
optimization  scheme  based  on  sensitivity  coefficients. 
Lee  and  Freudenstein  (48)  introduced  hueristic  combina- 
torial optimization  in  the  kinematic  design  of  mechanisms. 
A  hueristic  method  for  sorting  mechanism  variables  into 
independent  groups  was  proposed  by  Datseris  and 
Freudenstein  (49) .   Sutherland  and  Siddall  (50)  presented 
a  method  for  optimization  using  what  they  called  "inverse 
utilities"  as  a  basis  for  comparing  the  undesireable 
characteristics  of  a  mechanism.   Spitznagel  (51)  and 
Spitznagel  and  Tesar  (52)  developed  a  very  effective 
technique  for  optimizing  planar  mechanisms  based  on 
Burmester  synthesis  followed  by  sequential  filtering. 

Other  interesting  and  important  problems  concerning 
planar  mechanisms  were  also  being  addressed  in  the 
literature.   Optimal  synthesis  of  six-link  and  other 


multi-loop  mechanisms  was  studied  by  Chen  and  Dalsania 
(53) ,  Prasad  and  Bagci  (54) ,  Sallam  and  Lindholm  (55) , 
Mariante  and  Willmert  (56)  and  Spitznagel  (57).   Dhande 
and  Chakraborty  (58)  were  the  first  to  use  stochastic 
methods  to  optimize  mechanisms  considering  the  effects 
of  tolerances  and  clearances.   Sevak  and  McLarnan  (59) 
considered  the  problem  of  optimizing  mechanisms  with 
flexible  links.   Huey  and  Dixon  (60)  studied  the 
optimization  of  cam-modulated  linkages  for  path  and 
function  generation.   Kramer  and  Sandor  (61)  and  Kramer 

(62)  developed  an  optimization  technique  based  on  what 
they  termed  "selective  precision,"  wherein  structural 
errors  of  the  mechanism  are  held  within  a  prescribed 
tolerance  at  each  of  the  prescribed  positions.   Sutherland 

(63)  presented  a  method  for  designing  four-bar  linkages 
where  some  of  the  prescribed  positions  are  satisfied 
exactly  and  the  remaining  positions  are  approximated  in 
the  least-squares  error  sense.   Chouby  and  Rao  (64) 
suggest  a  direct  method  for  minimizing  the  structural 
error  together  with  the  mechanical  error  due  to  manu- 
facturing tolerances  on  the  link  dimensions. 

Application  of  numerically  based  optimization  methods 
to  the  design  of  spherical  and  spatial  mechanisms  became 
a  topic  of  interest  in  the  early  1970 's.   Stridher  and 
Torfason  (65)  minimized  the  structural  error  in  a 
spherical  four-link  path-generating  mechanism.   Bagci 
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and  Parekh  (66)  demonstrated  the  optimal  design  of 
spherical  four-bar  and  six-bar  linkages.   Bagci  (67) 
designed  spherical  four-link  mechanisms  to  have  optimal 
transmission  properties.   Rao  and  Ambekar  (68)  considered 
the  design  of  spherical  four-link  function-generating 
mechanisms  with  minimum  structural  error  subject  to 
link  length  and  transmission  angle  constraints. 

Early  efforts  to  optimize  spatial  mechanisms  were 
concerned  mainly  with  either  minimizing  structural  error 
or  maximizing  force  or  load  transmission,  or  some  combi- 
nation of  the  two.   Research  in  this  category  includes 
the  work  of  Sutherland  and  Siddall  (50) ,  Hamid  and  Soni 
(69)  and  Shoup,  Steffen  and  Weatherford  (70) .   Gupta 
(71)  demonstrated  the  synthesis  of  RSSR  and  RSRC  spatial 
linkages  with  minimum  structural  error,  subject  to 
branching,  mobility  and  transmission  constraints.   Suh 
and  Mechlenburg  (72)  synthesized  path-generating  spatial 
mechanisms  with  minimum  structural  error  in  the  least- 
squares  sense.   In  a  unique  work,  Bagci  (73)  optimized 
screw-generating  spatial  mechanisms.   Alizade,  Rao  and 
Sandor  (74,75)  demonstrated  the  synthesis  of  spatial 
function-generating  mechanisms  with  minimum  error  subject 
to  transmission  angle  constraints.   Suh  and  Radcliffe 
(76)  used  nonlinear  programming  methods  for  the  precision 
position  synthesis  of  spatial  mechanisms.   Although 
closed-form  solutions  were  possible  for  some  of  the 
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mechanisms  synthesized  by  Suh  and  Radcliffe,  their  work 
demonstrates  that  it  is  possible  to  circumvent  much  of 
the  labor-intensive  manipulation  of  algebraic  equations 
if  one  is  willing  to  pay  for  more  computer  time. 
Bagci  and  Falconer  (77)  optimized  the  transmission 
characteristics  of  RSSR  and  RSSP  function-generating 
spatial  mechanisms.   Soylemez  and  Freudenstein  (78) 
found  the  optimum  dimensions  of  an  RSSR  function- 
generating  mechanism  when  the  extreme  positions  of  the 
input  and  output  links  are  given.  Karelin  (79)  obtained 
formulas  for  determining  the  optimum  link  sizes  and 
slider  offset  for  the  RSSP  mechanism  given  the  required 
stroke.   His  design  was  subject  to  constraints  on  the 
pressure  angle. 

A  number  of  other  papers  which  do  not  fall  directly 
in  the  category  of  optimization  are,  nevertheless, 
important  to  the  present  study.   These  can  be  classified 
into  two  groups :   papers  dealing  with  determining  the 
quality  of  an  existing  mechanism  and  papers  dealing 
with  precision  position  synthesis  of  spatial  mechanisms. 

The  first  group  is  important  to  this  study  because 
all  mechanism  optimization  techniques  are  based  on 
quality  comparisons.   Perhaps  the  simplest  and  best 
known  quality  criterion  is  the  Grashof  mobility  test 
(80)  for  planar  four-bar  linkages.   Filemon  (81)  extended 
the  usefulness  of  this  test  by  determining  the  regions 
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of  mobility  along  Burmester's  centerpoint  curve. 
Jenkins,  Crossley  and  Hunt  (82)  studied  the  gross 
motion  attributes  of  certain  spatial  mechanisms  based 
on  the  intersections  of  surfaces  generated  by  points 
within  the  mechanism.   Duffy  and  Gilmartin  (83,84,85) 
generalized  Grashof ' s  work  by  determining  the  limit 
positions  and  the  mobility  of  four-link  spatial  and 
spherical  mechanisms  using  the  laws  of'  spatial  and 
spherical  triangles.   Gupta  and  Radcliffe  (86)  used 
geometric  methods  and  design  charts  to  determine  the 
mobility  of  planar  and  spatial  mechanisms.   Waldron 
and  Stevensen  (87)  worked  on  the  development  and 
application  of  conditions  on  branching,  mobility  and 
order  of  positions  in  planar  four-bar  linkages.   Strong 
and  Waldron  (88)  used  joint  displacements  to  determine 
the  mobility  regions  of  four-bar  linkages.   Sutherland 
(89)  developed  an  index  for  determining  the  quality  of 
force  and  motion  transmission  of  planar  and  spatial 
mechanisms.   Gupta  (90,91)  developed  theories  for 
synthesizing  crank-type  planar  four-bar  linkages  with 
transmission  angle  control.   Gupta  and  Tinubu  (92) 
proposed  a  theory  for  synthesizing  planar  and  spatial 
bimodal  function  generating  mechanisms  which  were  free 
from  branching  problems.   Zhuang  and  Sandor  (93,94) 
determined  the  branching  condition,  for  a  variety  of 
spatial  mechanisms  containing  spheric  joints. 
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Several  early  works  dealing  with  the  precision 
position  synthesis  of  spherical  and  spatial  mechanisms 
have  already  been  discussed.   Although  the  subject  is 
too  broad  to  fully  cover  here,  a  few  additional  references 
are  cited  to  give  the  reader  an  idea  of  the  methods 
currently  available. 

Dimentberg  (95)  introduced  the  screw  method  to  the 
study  of  lower-pair  spatial  mechanisms.   This  elegant  and 
powerful  tool  provided  a  concise  method  for  formulating 
and  solving  spatial  mechanism  problems.   Often  the  screw 
method  provided  closed-form  solutions  to  problems  which 
otherwise  require  numerical  solutions  when  formulated 
using  vector  or  matrix  methods.   Sandor  (96)  and  Sandor 
and  Bisshopp  (97)  introduced  methods  of  dual  number 
quaternions  and  stretch-rotation  tensors  to  find  the 
loop-closure  equations  of  spatial  mechanisms.   Beran  (98) 
used  dual  complex  numbers  to  synthesize  the  RCCC 
mechanism  for  multiply  separated  positions.   Tsai  and 
Roth  (99)  used  screw  triangle  geometry  to  synthesize 
open-loop  kinematic  chains.   Kohli  and  Soni  (100) 
developed  spatial  mechanism  synthesis  procedures  based 
on  pair  geometry  constraints  and  successive  screw 
displacements.   Suh  (101)  employed  4x4  matrices  in  the 
synthesis  of  spatial  mechanisms. 

Recently,  Sandor,  Kohli,  Reinholtz  and  Ghosal  (102) 
developed  an  analytical  closed-form  method  for  the 
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synthesis  of  a  spatial  motion-generating  mechanism 
for  three  prescribed  positions  based  on  vector  geometry. 
Sandor,  Kohli,  Zhuang  and  Reinholtz  (103)  extended  this 
work  to  four  prescribed  precision  positions;  their 
solution  involves  finding  the  simultaneous  solutions 
of  two  cubic  equations  in  three  unknown  variables.   The 
value  of  one  of  the  variables  is  assumed  arbitrarily  in 
solving  the  equations. 

While  fairly  extensive,  this  literature  review  is 
by  no  means  complete,  and  the  reader  is  encouraged  to 
peruse  the  excellent  review  articles  by  Fox  and  Gupta 
(104)  and  Root  and  Ragsdell  (105)  . 

1 . 4   Conclusions  of  the  Literature  Review 

The  foregoing  literature  review  reveals  a  number  of 
important  points: 

(1)  In  recent  years,  optimization  theory 
has  been  widely  recognized  as  an 
important  and  natural  tool  to  be 
used  in  the  design  of  mechanisms. 

(2)  No  single  mechanism  optimization 
procedure  has  received  widespread 
acceptance.   This  is  partially  due 
to  the  newness  of  the  optimization 
methods  themselves,  and  partially 
due  to  the  wide  variety  of  problems 
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which  occur  in  mechanism  design. 

(3)  In  a  majority  of  the  papers  reviewed 
dealing  with  optimization,  the 
objective  function  was  in  some  way 

a  measure  of  the  structural  error, 
although  no  reason  for  this  choice 
was  generally  given.   It  is  believed 
that  this  approach  arose  as  a 
"natural  extension"  of  the  methods 
of  precision  position  synthesis. 

(4)  The  authors  who  chose  to  optimize  with 
respect  to  properties  other  than 
structural  error  most  often  incorporated 
either  transmission  quality  or  some 
dynamic  property  of  the  mechanism  in 
the  objective  function.   Quite  often 

in  these  works  equality  constraints 
were  used  to  force  the  solution 
mechanism  to  pass  through  a  small 
number  (usually  3  or  4)  of  precision 
positions. 

(5)  Very  few  designers  have  attempted  to 
devise  a  general  philosophy  or  strategy 
for  mechanism  optimization.   As  a  result, 
many  of  the  methods  developed  have 
limited  application  or  are  relatively 
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inefficient  or  both. 
(6)   The  literature  contains  a  number  of 

very  fine  examples  of  planar  mechanism 
optimization.   Also,  several  useful 
techniques  have  been  developed  for  the 
synthesis  and  analysis  of  spatial 
mechanisms.   However,  the  extension  of 
these  techniques  to  spatial  mechanism 
optimization  has  been  quite  limited. 


CHAPTER  2 
OPTIMIZATION  THEORY 

2  . 1   Introduction  to  Optimization 

Optimization  has  previously  been  defined  as  the  act 
of  seeking  the  best  result  under  a  given  set  of  circum- 
stances.  Finding  the  best  result  ultimately  means 
minimizing  something  (such  as  the  required  effort) ,  or 
maximizing  something  (such  as  the  desired  benefit) .   In 
order  to  have  a  solvable  optimization  problem,  the 
desired  benefit  or  the  required  effort  must  be  express- 
ible as  a  function  of  a  set  of  variables  over  which  the 
designer  has  control.   These  variables  are  called  the 
design  variables.    Limits  on  the  values  of  design 
variables  may  result  from  such  things  as  limited 
material  supplies  or  limited  production  capabilities. 
In  either  case,  these  limiting  factors  are  called 
constraints,  and  they,  too,  must  be  expressible  as 
functions  of  the  design  variables. 

From  figure  2.1  it  can  be  seen  that  the  maximum 
value  of  a  function,  f  (x)  ,  corresponds  to  the  minimum 
value  of  the  negative  of  the  function,  -f(x).   Therefore, 
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iigure  2.1   Example  showing  the  maximum  of  f (x)  to  be  the 
same  as  the  minimum  of  -f (x) 
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no  generality  is  lost  by  assuming  the  optimization  problem 
to  always  be  one  of  function  minimization. 

2 . 2   Formal  Definition  of  the  Optimization  Problem 

An  optimization  problem  or  a  mathematical  programming 
problem  can  be  stated  in  the  following  general  form 
(104,106)  : 


Find        X  =  <  *  >  (2.1) 


which  minimizes  f(X,0)   (where  0  is  an  independent  set 
of  input  motion  variables  having  a  predetermined  range) , 
subject  to  inequality  constraints 

g. (X,0)  <  0,      j  =  1,2, ...,m  (2.2) 

and  equality  constraints 

£,(X,0)  =0,      j  =  m+1,  m+2,  ...,  p  (2.3) 

Here  X  is  called  the  design  vector,  and  f(X,0)  is  called 
the  objective  function.   The  design  vector,  X,  is  composed 

of  the  design  variables  x, ,  x- ,  ...,  x  ,  which  may  be 

■,  T 
written  in  transposed  lxn  matrix  form  as  {x  , x2 ,  .  . . , x^}  , 

where  the  superscript  T  denotes  the  transpose.   The  vector 

T 

0  =  {S,,6_,...,8  }   of  a  set  of  independent  parameters 
12       q 
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may  typically  comprise  quantities  such  as  time  or 
position.   Most  types  of  optimization  problems  do  not 
contain  this  independent  parameter  vector,  and,  hence, 
most  textbooks   do  not  include  it  in  defining  the 
standard  optimization  problem.   However,  in  mechanism 
design,  the  independent  parameter  (s)  representing  the 
mechanism  input (s)  often  appear  in  the  design  equations. 

When  the  vector  9  is  present  in  the  objective  function 
or  in  any  of  the  constraint  equations,  the  problem  will 
be  called  a  parametric  programming  problem.   This  type 
of  problem  is  generally  much  more  difficult  to  solve  than 
nonparametric  problems,  because  for  each  new  design  vector, 
X,  that  is  tried,  the  objective  function  and  the  constraints 
must  be  evaluated  at  every  possible  value  of  the  independent 
input  parameter (s) .   In  this  dissertation,  a  great  deal 
of  emphasis  will  be  placed  on  avoiding  the  need  to  involve 
these  independent  parameters,  while  assuring  optimal 
performance  throughout  their  range. 

Finally,  it  should  be  pointed  out  that  the  components 
of  the  design  vector,  X,  may  also  be  functions  of  inde- 
pendent parameters.   This  occurs,  for  example,  when  link 
shapes  are  considered  variables  in  the  design  of  a 
mechanism  which  is  to  have  maximum  stiffness.   This  type 
of  optimization  problem  is  known  as  a  dynamic  programming 
problem.   No  dynamic  programming  problems  have  been 
included  among  the  examples  of  this  dissertation. 
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2. 3   The  Mechanism  Optimization  Problem 

To  put  the  preceding  definition  of  the  optimization 
problem  in  a  clearer  light,  it  will  now  be  cast  in  the 
framework  of  application  to  the  design  synthesis  of  a 
simple  mechanism. 

Suppose  it  is  desired  to  design  a  function-generating 
planar  four-bar  linkage  for  a  number  of  arbitrarily 
prescribed  positions  greater  than  five.   Since,  in  this 
case,  five  is  the  maximum  number  of  positions  that  can 
be  satisfied  exactly  (unless  input  and  output  scale 
factors  are  made  components  of  the  design  vector) ,  the 
best  that  can  be  hoped  for  is  to  minimize  the  structural 
error  at  the  prescribed  positions.   A  logical  and  popular 
approach  to  this  problem  has  been  to  minimize  the  sum  of 
the  squares  of  the  resulting  structural  errors  at  the 
prescribed  positions.   Thus,  if  the  output  angle,  cf> .  ,  is 
some  known  function,  F,  of  the  input  angle,  0.,  at  each 
of  the  n  prescribed  positions,  then 

<j>i  =  F(6i)  ,       i  =  1,2,.  ..n,  (2.5) 

and  the  objective  function  (O.F.)  to  be  minimized  is 

n  2 

O.F.  =   E  {F(8.)  -  G(9.) }  (2.6) 

i=l     X  1 

where  G(8.),  i  =  1,2,.. . ,n  are  the  output  positions 
actually  generated  by  the  mechanism.   The  mechanism  to 
be  optimized  is  shown  in  figure  2.2. 
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figure  2.2   The  function-generating  planar  four-bar 
linkage 
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The  design  vector  consists  of  the  design  variables 

which  are  the  link  lengths  A,B,C,D  and  the  starting 

angles  8n  and  <j>  .   For  a  given  value  of  6=8.,  the  values 

of  these  variables  determine  the  values  of  F(0.)  and 

i 

G(9.).   The  values  of  the  variables  A,B,C,D,9Q  and  t}>Q 
are  to  be  found  which  minimize  the  function  2.6. 

The  reader  familiar  with  kinematics  will  know  that 
the  relative,  rather  than  the  absolute,  proportions  of 
this  mechanism  determine  the  input/output  functional 
relationship.   Therefore,  D,  the  length  of  the  grounded 
link,  can  be  set  equal  to  unity  without  loss  of  generality. 
The  design  vector  now  becomes 


X  -    /  :c,  \    -    /  C   S   =  {A,B,C,80,d>0}T     (2.7) 


While  finding  the  value  of  X  which  minimizes  the 
function  2.6  may  be  quite  an  interesting  problem  in 
itself,  in  practice  there  are  usually  a  number  of  other 
requirements  the  linkage  must  satisfy.   For  example,  the 
transmission  angle  ,  labeled  \i    in  figure  2.2,  must  be 


The  transmission  angle  in  a  four-bar  mechanism  is  defined 
as  the  acute  angle  between  the  coupler  link  and  the  output 
link. 
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held  reasonably  close  to  90°  to  ensure  effective  conver- 
sion of  the  force  in  the  coupler  to  torque  about  the 
output  crank  pin.   In  practice,  a  value  of  u  less  than 
about  30   will  generally  be  unacceptable.   This  require- 
ment can  be  expressed  in  the  form  of  an  inequality 
constraint  equation 


M    >  30( 


y  -  30   >  0 


(2.8) 


However,  this  is  not  an  acceptable  form  of  the  equation 
because  the  transmission  angle,  u,  is  not  a  design 
variable  or  an  independent  parameter.   However,  u  can 
be  expressed  in  terms  of  the  design  variables  and  the 
independent  parameter  0  as  follows  (107,pp319): 


H  =  arccos 


A2  -  B2  -  C2  +  1  -  2Acos(9  +  9  ) 
2BC  U 


(2.9) 


Upon  substituting  this  result  into  the  transmission  angle 
constraint,  equation  2.8  becomes 


A2  -  32  -  C2  +  1  -  2Acos(9  +  9  J 

0 


23C 


-  30"  >  0 


(2.10) 


Unfortunately,  this  is  a  parametric  constraint  since 
it  contains  the  independent  parameter  9.   It  therefore 
seems  that  equation  2.10  must  be  evaluated  at  every  possible 
value  of  9,  from  9=0°  to  9=360°,  for  every  set  of  design 
variables  tried,  to  be  sure  the  constraint  is  satisfied. 
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Fortunately,  this  is  not  usually  necessary  because,  as 
shown  by  Roth,  Freudenstein  and  Sandor  (21) ,  the  trans- 
mission angle,  y ,  is  a  maximum  when  8  +  6_  =  180°  and  a 
minimum  when  9  +  8„  =  0° ,  provided  these  positions  are 
real,  i.e.,  if  the  linkage  closes  in  these  positions  of 
the  input  link.   This  demonstrates  an  extremely  important 
concept:   the  need  for  parametric  constraints  can 
sometimes  be  eliminated  by  determining  the  critical 
values  of  the  independent  parameter. 

Nonparametric  constraints  can  result  from  a  number 
of  requirements.   For  example,  it  is  usually  desirable  to 
limit  the  ratio  of  link  lengths  within  a  mechanism; 
otherwise  solutions  will  result  with  extremely  long 
coupler  links  and  short  cranks  (for  practical  purposes, 
these  become  slider-crank  mechanisms) .   To  eliminate  these 
unwanted  solutions,  the  following  constraints  are  specified: 

A„  <  A  <  A 
I   —        —   u 

B,  <  B  <  B 

I   —        —     u 

C„  <  C  <  C  (2.11) 

I   —   —   u 

or,  expanding  these  into  the  standard  form  of  separate 
inequalities , 


Ai  ~  A  £  °  A  "  Au  -  ° 

B£    -    B    <  0  B-Bu<0 

r      -   C    <  0  C-C<0  (2.12; 

L  —  U      — 
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where,  again,  the  length  of  the  fixed  link,  D,  is  arbi- 
trarily set  equal  to  unity,  and  the  subscripts  i   and  u 
refer  to  the  lower  and  upper  limits  of  the  link  length. 

Equality  constraints  in  linkage  synthesis  problems 
most  often  result  from  precision  position  requirements. 
Referring  to  the  notation  used  in  equation  2.6,  suppose 
that,  at  the  three  positions  q,r  and  s,  the  actual 
mechanism .output  is  required  to  be  the  same  as  the 
specified  output,  then 


F(9i)  =  G(8i) ,      i  =  q,r,s  (2.13) 


F(8i)  -  G(9i)  =0,      i  =  q,r,s  (2.14) 

are  the  needed  equality  constraints.   Positions  q,r  and  s 
may  now  be  removed  from  consideration  in  equation  2.6, 
although  this  is  not  mandatory  since  they  do  not  contribute 
to  the  value  of  the  objective  function.   An  important 
point  which  will  be  studied  in  detail  in  section  2.3.2 
is  that  equality  constraints  can  very  often  be  used  to 
eliminate  variables  from  consideration,  without  any  loss 
of  generality  in  the  problem  formulation. 

Summarizing  the  four-bar  linkage  optimization  problem 
at  hand , 
find 

X  =  {A,B,C,G0,(»0}T  (2.15) 


which  minimizes 


Z{F(8i)  -  G(9i)  }  , 
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i  =  1,2, . . . ,q, . . . ,r, . . . ,s, 


(2.16; 


subject  to 


2     2     2 
A"  -  B   -  C 


1  +  2A 


2BC 


-  30°  >  0 


(2.17; 


arccos 


2     2     2 

A   -  B   -  C   +  1  -  2A 

2BC 


-  30°  >  0 


(2. IS 


A,  -  A  <  0 
x,      — 


A  -  A   <  0 
U  — 


B„  -  B  <  0 


B  -  B   <  0 

u  — 


cz  -   C  <  0 


c  -  c  <  o 

u  — 


(2.19; 


and 


F(9i)  -  G(6i)  =  0, 


l  =  q,r,s 


(2.20) 


Notice  that  the  critical  values 


+  8  )  =  0°  and 


(9  +  9n)  =  180°  have  been  substituted  into  constraint 
equation  2.10  to  yield  equation  2.17  and  2.13.  Also, 
the  positions  at  which  the  function  must  be  exactly 
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satisfied,  6.,  i  =  q,r,s,  have  not  been  removed  from 
consideration  in  the  objective  function. 

A  number  of  other  constraints  could  be  added  to 
this  problem.   For  example,  it  is  often  desirable  to 
have  a  crank-rocker  type  linkage,  which  can  be  driven 
from  a  continuously  rotating  prime  mover.   The  inequality 
constraint  which  ensures  this  type  of  mechanism  is 
derived  from  the  so-called  Grashof  criterion.   Other 
constraint  conditions  will  be  discussed  in  Chapter 
Three.   The  present  problem  is  sufficiently  general  for 
use  in  discussing  the  various  methods  available  for 
solving  the  general  nonlinear  optimization  problem  which 
typically  results  from  mechanism  optimization;  these 
solution  methods  are  presented  in  the  following  section. 

2. 4   Solving  the  Mechanism  Optimization  Problem 

Optimization  problems  may  be  either  linear  or  non- 
linear.  Nonlinear  problems  occur  when  the  objective 
function  or  any  of  the  constraint  equations  are  nonlinear 
in  any  of  the  design  variables.   It  will  not  surprise  the 
experienced  kinematician  to  know  that  almost  all 
mechanism  optimization  problems  fall  into  this  category. 
Even  judging  from  the  simple  example  of  the  previous 
section,  it  becomes  evident  that  these  types  of 
problems  are  generally  nonlinear,  and  usually  involve 
a  large  number  of  constraints.   For  this  reason,  the 
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remainder  of  this  section  will  deal  exclusively  with 
methods  for  solving  constrained,  nonlinear  optimization 
problems . 

2.4.1   Constrained  Nonlinear  Optimization 

A  number  of  methods  are  currently  available  for 
solving  constrained  nonlinear  optimization  problems. 
While  some  of  these  methods  are  more  widely  used  than 
others,  no  single  method  is  best  suited  to  solve  every 
type  of  problem.   In  fact,  even  small  changes  in  the 
way  a  problem  is  formulated  can  grossly  alter  the  effec- 
tiveness of  the  optimization  procedure  being  used.   This 
explains,  in  part,  the  large  number  of  past  approaches 
which  have  been  taken  in  mechanism  optimization,  and 
the  claims,  often  conflicting  or  confusing,  about  the 
efficiency  and  the  effectiveness  of  these  approaches 
(104,105)  . 

Although  it  is  possible  to  solve  some  constrained 
nonlinear  optimization  problems  using  the  classical 
techniques  of  variational  calculus,  the  complexity  of 
most  mechanism  optimization  problems  renders  this  approach 
impractical.   Therefore,  the  remainder  of  this  section 
will  be  devoted  to  the  iterative  solution  methods 
known  as  constrained  nonlinear  mathematical  programming, 
or  simply  constrained  nonlinear  programming. 
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Constrained  nonlinear  programming  techniques  can 
be  divided  into  two  distinct  groups:   direct  methods  and 
indirect  methods.   Table  2.1  shows  the  various  techniques 
that  come  under  these  headings  (106) . 

Table  2.1 
Constrained  Nonlinear  Programming  Techniques 

Direct  Methods  Indirect  Methods 

1)  Heuristic  search  1)   Transformation  of 

variables 

2)  Constraint  approximation 


3)   Feasible  directions 


2)   Penalty  functions 


Classification  is  based  on  the  manner  in  which  the 
constraints  are  handled.   The  direct  methods  deal  with 
the  constraints  explicitly,  whereas  the  indirect  methods 
first  transform  the  constrained  problem  into  an  uncon- 
strained problem,  and  then  solve  this  new  easier  problem 
using  one  of  the  unconstrained  nonlinear  programming 
methods . 

All  of  the  methods  listed  in  Table  2.1  are  poten- 
tially useful  in  mechanism  optimization.   However,  it 
will  not  be  practical  to  discuss  each  of  these  methods 
in  detail.   Since  the  indirect  methods  are  used  in 
working  the  examples  of  this  dissertation,  these  will 
be  discussed  in  greater  detail.   Complete  details  on 
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all  of  the  methods  listed  are  available  in  the 
references  (106,108,109). 

2.4.2   Indirect  Constrained  Nonlinear  Programming  Techniques 

As  can  be  seen  from  Table  2.1,  two  methods  come  under 
this  heading:   transformation  of  variables  and  penalty 
functions. 

In  the  transformation  of  variables  technique,  the 
design  variables  are  changed  in  such  a  way  as  to  auto- 
matically ensure  constraint  satisfaction.   There  are 
two  cases  where  this  is  possible:   (1)  when  the  constraints 
are  simple,  explicit  functions  of  the  decision  variables, 
and  (2)  when  equality  constraints  can  be  used  to  eliminate 
variables . 

As  an  example  of  the  first  case,  recall  the  first 
of  the  link  length  constraint  equations  2.11  of  section 
2.2,  repeated  here  as  equation  2.21. 

A0  <  A  <  A  (2.21) 

I   —   —  u 

where  A  was  the  length  of  the  input  crank  of  the  four- 
bar  linkage  of  figure  2.1,  and  A„  and  A   were  the  upper 
and  lower  limits  placed  on  this  length.   These  constraints 
can  be  automatically  satisfied  by  transforming  the 
variable  A  to  the  form 

A  =  A,  +  (A   -  AJsinV,  (2.22) 

I  u     I 
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where  A*  is  the  new  variable  which  can  take  on  any 

numerical  value.   This  technique  is  called  change  of 

variables.   Notice  the  A  will  always  be  between  A,  and 

A   for  any  value  of  A*  in  eauation  2.22. 
u       2 

While  this  approach  seems  quite  promising  at  first, 
experience  has  shown  that,  unless  all  the  constraints 
can  be  transformed  in  this  way,  it  is  probably  better 
not  to  use  the  transform  at  all  (102).   This  is  because 
substitution  of  the  right-hand  side  of  equation  2.22 
into  the  objective  function  may  distort  it  to  the  point 
where  it  is  more  difficult  to  minimize  than  the 
original  function,  when  other  constraints  are  present. 
As  a  result,  this  approach  appears  to  be  impractical 
for  most  mechanism  optimization  problems,  although  the 
author  believes  further  study  in  this  area  is  warranted. 

The  second  case  for  which  transformation  of  variables 
is  sometimes  possible  occurs  when  equality  constraints 
are  present.   For  example,  in  section  2.2,  equation  2.13 
expresses  a  set  of  three  precision  position  requirements 
for  the  four-bar  linkage  of  figure  2.1 


F(9i)  =  G(9i) ,      i  =  q,r,s  (2.23) 

where  the  functions  F(9.)  and  G(9.)  express  the  desired 
and  the  generated  output  angular  positions  of  the  mechanism, 
respectively,  at  a  given  input  position  9..   Recall  from 
equation  2.5  that  F(9.)  =  cj> .  is  a  known,  prescribed 
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function.   Equation  2.23  may  be  rewritten  in  the  form  of 
the  well-known  Freudenstein  equation  (1) ,  which  gives 
the  input/output  angular  relationship  of  the  four-bar 
function  generator  in  terms  of  the  linkage  dimensions 

■ 

K1cos(ei  +  9Q)  -  K2cos((})i  +  <j>  )  +  K.,  = 

cos(6i  +  eQ  -  <()i  -  <j>0)  (2.24) 

where 


Kn  =  i,     K9  =  1;     K   =  A2  -  B2  +  C2  +  1      {2>25, 
1    C  A       J  2AC 


and  Gn  and  cf>n  are  the  starting  angles  (see  figure  2.2). 
Equation  2.24  is  linear  in  the  coefficients  K, ,K„  and  K.., 
and  can  be  written  three  times  corresponding  to  the  three 
values  of  i  (i.e.  i  =  q,r  and  s) .   It  is  therefore  a 
relatively  easy  matter  to  solve  for  the  values  of  K,,  K„ 
and  K,  in  terms  of  the  angles  9.,  <b .  ,   i  =  q,r,s,  and 
8Q,  <)>0  (101).   With  K,,  K2  and  K3  known,  the  link  lengths 
may  be  determined  from  the  relations 


A  =  -  ;   C  =  -  ;    B  =  {A2  +  C2  +  1  -  2ACK3}*5     (2.26) 
K2        K1 


Notice  that  the  only  design  variables  remaining  in  the 
right-hand  sides  of  equations  2.26  are  the  starting  angles 
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6   and  tj>Q.   Anywhere  A,B  and  C  appear  in  the  remaining 
constraint  equations  or  in  the  objective  function  they 
may  be  replaced  by  the  equivalent  expressions  given  in 
equation  2.26.   This  technique  is  called  elimination  of 
variables. 

As  in  the  previously-described  change  of  variables 
technique,  the  objective  function  and  the  remaining 
constraints  will  undoubtedly  be  distorted  when  the  non- 
linear expressions  for  A,B  and  C  of  equation  2.26  are 
substituted  into  them.   However,  experience  has  shown 
that  it  is  generally  beneficial,  in  this  case,  to  make 
the  substitution  (26,51,52,63,106).   Obviously,  the 
greater  the  degree  of  nonlinearity  of  the  expressions 
used  to  eliminate  variables,  the  less  beneficial  this 
type  of  substitution  becomes. 

It  has  already  been  noted  that  the  role  of  the 
indirect  methods  of  constrained  nonlinear  programming  is 
to  transform  the  constrained  problem  into  an  equivalent 
unconstrained  problem.   Occasionally,  the  transformation 
of  variables  method  can  accomplish  this  goal  by  itself. 
More  often,  however,  some  of  the  constraint  equations 
are  too  complex  for  this,  and  the  so-called  penalty 
function  approach  must  be  employed. 

To  demonstrate  the  penalty  function  approach, 
consider  the  following  simple  problem,  where  the  parametric 
vector,  9,  has  been  omitted  for  clarity. 
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Find  X  which  minimizes  f (X) 
subject  to 


g.  (X)  <  0,      j  =  1,2,... ,m  (2.27) 

and 


I.  (X)  =  0,      j  =  m+l,m+2, . . . ,p 


This  constrained  problem  is  converted  into  an  unconstrained 
problem  by  constructing  a  new  function  to  be  minimized  of 
the  form 


U  =  U(X,r)  =  f(X)  +  r  Z  G.{g.(X)} 

j=l  J       J    " 

P 
+  r   £   L.U.  (X)  }  (2.28) 

j=m+l  J      3    ~ 


where  G.{g.(X)}  and  L.{£.(X)}  are  functions  of  the  constraint 
functions  g.(X)  and  I .  (X) ,  respectively  and  r  is  a  positive 
constant  called  the  penalty  parameter.   The  solution  of 
the  unconstrained  problem  of  equation  2.28  can  be  made  to 
converge  to  the  solution  of  the  original  problem  of 
equation  2.27  by  repeating  the  minimization  process  for 
a  progressively  larger  series  of  values  of  the  penalty 
parameter,  r.   For  this  reason,  the  penalty  function 
methods  are  often  referred  to  as  "sequential  unconstrained 
minimization  techniques"  or  simply  SUMT . 
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Two  categories  of  penalty  function  methods  exist, 
namely,  interior  methods  and  exterior  methods.   The 
interior  methods  must  be  supplied  with  a  feasible 
starting  vector,  X,  (i.e.  g  (X. )  <  0  for  all  j).   As 
the  parameter  r  is  varied  over  successive  minimizations, 
the  solution  of  the  unconstrained  problem  converges  to 
the  solution  of  the  constrained  problem,  always  remaining 
within  the  feasible  region.   Since  the  search  is  conducted 
within  the  feasible  region,  these  are  called  interior 
methods.   The  exterior  methods  do  not  require  a  feasible 
starting  vector,  and  generally  converge  to  the  constrained 
minimum  from  outside  the  feasible  region,  hence  the 
term  exterior.   The  exterior  methods  have  been  judged 
to  be  generally  superior  to  the  interior  methods  (109)  ; 
therefore,  these  will  be  reviewed  in  greater  detail. 

A  typical  exterior  penalty  function  form  of  equation 
2.28  is 


m 
U(X,r)  =   f(X)  +  r 
j 


2  (g  (X)) 

=1  v  J  ~  ' 


+  r   E   U.  (X)}z  (2.29 

j=m+l  3    ~ 


where,  again,  r  is  a  positive  penalty  parameter,  q  is  a 
constant  greater  than  one,  and  the  singularity  function 
/g.(X)y   is  defined  by 
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g.  (X)   if   g.  (X)  >  0 


(gjW)  =   ^  (2.30) 

if   g  .  (X)  <  0 


It  can  be  seen  that  the  effect  of  this  formulation 
is  to  assess  an  increasingly  severe  penalty  on  the  value 
of  U(X,r)  as  the  constraints  become  violated  by  larger 
amounts.   The  most  successful  way  to  find  the  true 
constrained  minimum  of  the  original  function  has  generally 
been  to  minimize  equation  2.23  using  a  small  value  for 
the  parameter  r  for  the  first  minimization.   Subsequent 
minimizations  use  successively  larger  values  of  this 
parameter,  until  the  solution  is  essentially  forced  to 
converge  in  the  feasible  region.   This  is  necessary 
because,  if  r  remained  small,  very  small  positive  values 
of  g.  (X)  ,  even  though  infeasible,  would  not  contribute 
much  of  a  penalty  to  U(X,r)  and  the  solution  might  remain 
infeasible.   On  the  other  hand,  if  r  were  made  initially 
very  large,  the  contribution  of  f(X)  would  be  negligible, 
and  the  solution  may  not  converge  to  the  minimum  f (X) 
within  the  constrained  region.   In  other  words,  the 
effects  should  be  balanced,  so  that  the  solution  is 
urged  toward  the  minimum  of  f(X)  at  the  same  time  it  is 
being  forced  toward  the  feasible  region.   The  simple 
example  which  follows  will  help  to  demonstrate  some 
of  these  concepts. 
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Find  X  =  {x, }  which  minimizes  f(X)  =  ^x, 
subject  to  the  constraint 


3  -  x1   £  0  (2.31) 

The  objective  function  f(X)  and  the  constraint 
boundary  are  plotted  in  figure  2.3.   The  constrained 
minimum  is  clearly  at  x.  =3.   Now  construct  the  exterior 
penalty  function 

U(X,r)  =  Hx1    +    r/3  -  xx\2  (2.32) 

Table  2.2  gives  various  values  of  U(X,r)  versus  x.  for 
several  values  of  r. 

The  resulting  curves  are  plotted  in  figure  2.3.   It 
is  clear  from  this  figure  that,  as  r  tends  toward  infinity, 
the  solution  to  the  unconstrained  objective  function  of 
equation  2.32  will  approach  the  solution  to  the  constrained 
problem  of  equation  2.31.   It  may  be  noted  that,  in  this 
case,  the  solution  will  reach  the  feasible  region  only  in 
the  limit  as  r  approaches  infinity.   This  usually  is  not 
troublesome  for  practical  problems  because  the  constraints 
are  rarely  known  exactly,  and  some  allowance  must  be  made 
for  errors  when  formulating  them. 

The  preceding  discussion  focused  on  the  use  of 
penalty  function  methods  and  change  of  variable  techniques 
to  transform  constrained  optimization  problems  into 
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Table  2.2 
Example  showing  the  effects  of  the  penalty  parameter,  r. 


U(X,r) 


xl 

r=0.25 

r=0.5 

r=l 

0 

2.25 

4.50 

9.00 

0.25 

2.02 

3.90 

7.69 

0.50 

1.81 

3.38 

6.50 

0.75 

1.64 

2.91 

5.44 

1.00 

1.50 

2.50 

4.50 

1.25 

1.39 

2.16 

3.69 

1.50 

1.31 

1.86 

2.75 

1.75 

1.27 

1.66 

2.44 

2.00 

1.25 

1.50 

2.00 

2.25 

1.27 

1.40 

1.69 

2.50 

1.31 

1.38             j 

1.50 

2.75 

1.39 

1.41 

1     i-«  ! 

3.00 

1.50 

1.50 

1.50 

lold    squares  indicate  the  tabulated  minimum  of  U(X,r). 
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f  (X)  , 

U(X,r) 

infeasible  region 
,r  =  0.5    \r  =  1 


x,  =  3 


feasible  region 


figure  2.3   Penalty  function  example 
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unconstrained  ones.   The  question  still  remains  of  how 
these  unconstrained  problems  may  be  solved;  this  is  the 
topic  of  the  next  section. 

2.4.3   Unconstrained  Nonlinear  Optimization 

Even  though  the  vast  majority  of  nonlinear  optimiza- 
tion problems  involve  constraints,  most  of  the  available 
nonlinear  optimization  techniques  have  been  developed 
for  solving  unconstrained  problems.   This  does  not 
represent  a  serious  limitation,  however,  because  most  of 
these  methods  can  be  extended  to  handle  constrained 
problems,  either  by  directly  considering  the  constraints 
or  by  transformation  to  an  unconstrained  problem  as 
discussed  in  the  previous  section. 

Unconstrained  minimization  methods  may  be  divided 
into  two  groups:   direct  search  methods  and  descent  (or 
gradient)  methods.   The  gradient  methods  require  either 
an  analytical  or  a  numerical  derivative  of  the  objective 
function  with  respect  to  the  design  variables,  whereas 
the  direct  search  methods  do  not.   Some  of  the  commonly 
available  techniques  in  both  groups  are  listed  in  Table 
2.3  (106)  . 

Of  the  methods  listed  in  Table  2.3,  the  random 
search  and  the  grid  search  are  known  to  be  quite  ineffi- 
cient.  However,  these  methods  tend  to  be  reliable  when 
minimizing  discontinuous,  sharply  varying  or 


51 


nondif ferentiable  functions.   They  may  also  be  useful 
for  finding  feasible  solutions  to  initiate  some  of  the 
more  efficient  methods. 

Table  2.3 
Unconstrained  Minimization  Techniques 


Direct  Search  Methods 

1)  Random  search 

2)  Grid  search 

3)  Univariate  search 

4)  Pattern  search  (Powell's 

method,  Hooke  and  Jeeves' 
method) 

5)  Method  of  rotating  coor- 

dinates (Rosenbrock ' s 
method) 

6)  Simplex  method 


Descent  Methods 

1)  Steepest  descent 

method 

2)  Conjugate  gradient 

method  (Fletcher- 
Reeves  ) 

3)  Newton's  method 

4)  Variable  metric 

method  (Davidon- 
Fletcher-Powell) 


An  excellent  comparison  of  many  of  the  numerical 
optimization  methods  commonly  used  to  solve  mechanical 
design  problems  was  made  by  Eason  and  Fenton  (10  9) .   They 
point  out  that  the  ideal  computer  code  for  design  optimi- 
zation should  solve  any  problem  conveniently  and  at 
moderate  cost.   No  code  tested  by  them  fulfilled  this 
requirement,  but  one  method  did  stand  out  above  the 
others,  namely,  the  pattern  search  method  of  Hooke  and 
Jeeves.   A  number  of  other  important  conclusions  were 
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reached  in  this  study;  some  of  these  are  summarized 
below: 

(1)  Derivatives  of  the  objective  function  are 
often  difficult  or  impossible  to  calculate 
analytically  for  many  mechanical  design- 
type  problems,  and  must  be  approximated 

by  numerical  methods,  if  needed. 

(2)  If  derivatives  must  be  calculated  numeri- 
cally, the  direct  methods  (which  do  not 
require  derivatives)  are  generally  superior 
to  the  gradient  methods . 

(3)  Automatic  scaling  of  the  design  variables 
within  the  computer  program  generally 
increases  the  efficiency  of  an  algorithm. 

(4)  The  most  general  methods  (those  which 
could  solve  the  most  types  of  problems) 
were  not  necessarily  slow,  nor  did  they 
require  the  greatest  amount  of  computer 
code  to  program. 

(5)  The  cost  of  preparing  a  problem  for 
computer  solution  may  be  greater  than  the 
execution  cost.   The  algorithm  should, 
therefore,  be  convenient  to  use. 

(6)  A  computer  package  containing  an  assortment 
of  optimization  methods  would  generally  be 
preferred  to  any  single  method.   This 
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would  allow  cross-checking  of  results, 
and  obviously  would  allow  more  types  of 
problems  to  be  solved  than  would  any  of 
the  methods  individually. 

Development  of  a  computer  optimization  package  as 
described  in  conclusion  number  (6)  above  would  be  a 
major  undertaking.   Also,  such  a  package  would  require 
a  relatively  large  amount  of  computer  storage.   Still, 
this  would  be  the  preferred  approach  in  terms  of 
efficiency  and  generality.   Since  no  such  computer 
package  is  readily  available,  and  since  computer  storage 
will  often  be  an  important  consideration,  it  was  decided 
to  use  the  pattern  search  method  of  Hooke  and  Jeeves  in 
working  the  examples  of  this  dissertation.   Later  results 
will  show  this  method  to  be  quite  effective  in  solving 
mechanism  optimization  problems.   Good  results  using 
this  method  for  optimizing  planar  mechanisms  were  also 
reported  by  Kramer  and  Sandor  (61)  and  Kramer  (62). 

2.4.4   Hooke  and  Jeeves'  Nonlinear  Programming  Method 

Hooke  and  Jeeves'  search  method  is  a  direct,  sequen- 
tial stepping  technique  consisting  of  alternating 
exploratory  and  pattern  moves.   The  exploratory  move 
seeks  to  determine  the  local  behavior  of  the  objective 
function,  and  the  pattern  move  uses  this  information  in 
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an  attempt  to  "leapfrog"  to  an  improved  position.   Figure 
2.3  demonstrates  this  procedure  for  a  two-dimensional 
problem,  and  the  algorithm  is  outlined  in  detail  below. 

(1)   Starting  from  an  arbitrarily  selected  point 

0  0  T 

X  ,  located  at  X  =  {x  ,x2>  ,  try  an  explor- 
atory move  by  changing  x,  by  a  predetermined 
positive  step  Ax, .   Evaluate  the  objective 
function  at  this  new  point.   If  its  value  is 
improved,  this  step  is  retained,  and  becomes 
the  new  base  point.   If  the  value  of  the 
objective  function  is  not  improved,  a  negative 
step,  -Ax,,  is  taken  and  the  objective  function 
is  reevaluated.   If  the  move  is  successful, 
the  point  is  retained  as  the  new  base  point. 
If  both  steps  fail,  no  move  is  made.   At 
this  point  it  is  often  beneficial  to  adjust 
the  value  of  Ax, ,  and  save  this  information 
for  future  exploratory  searches  in  the  x, 
direction.   A  successful  step  move  would 
suggest  an  increase  in  the  value  of  Ax, .   If 
neither  step  were  successful,  the  value  of 
Ax,  should  be  decreased.   Starting  from  the 
best  point  found  in  the  x,  search  (labeled 
X   in  figure  2.3)  a  similar  search  is  made 

in  the  x~  direction.   The  best  point  found 

2 
in  this  search  is  labeled  X  . 
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figure  2.3   The  pattern  search  method  of  Hooke  and  Jeeves 
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(2)  A  pattern  move  is  now  attempted  by  repeating 

all  the  successful  moves  of  the  exploratory 

2 
search  from  point  X  .   After  the  initial 

search,  this  move  may  also  include  the 

previous  exploratory  moves  and  the  previous 

pattern  search. 

(3)  If  the  pattern  move  of  part  (2)  is  successful, 

3 
it  is  retained  as  point  X  ,  and  the  exploratory 

search  begins  at  this  point.   If  the  pattern 

move  fails,  the  next  exploratory  search  begins 

2 

at  point  X  . 

(4)  This  process  is  repeated  until  the  values  of 
the  Ax-s  are  below  a  certain  preset  limit. 
At  this  point  the  minimum  is  presumed  to 
have  been  reached. 


A  flow  chart  of  this  procedure  is  shown  in  figure  2.4. 

While  simple  in  concept,  the  Hooke  and  Jeeves' 
method  is  extremely  powerful.   It  should  be  cautioned, 
however,  that  none  of  the  optimization  methods  discussed 
can  discriminate  between  a  local  and  a  global  minimum. 
Geven  ten  different  starting  points,  the  Hooke  and  Jeeves' 
method  may  converge  to  ten  different  local  minima.   As 
discussed  at  length  in  the  next  chapter,  this  will 
seldom  represent  a  serious  problem  in  mechanism 
optimization. 
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CHAPTER  3 
PHILOSOPHY  OF  MECHANISiM  OPTIMIZATION 

3. 1   The  Need  for  a  General  Philosophy 

Although  theories  for  the  kinematic  synthesis  and 
analysis  of  many  types  of  spatial  mechanisms  are  well 
developed  and  readily  available,  complete  design  theories 
for  even  the  simplest  four-link  spatial  mechanisms  are  not 
yet  developed  to  the  point  where  they  are  practical  for 
use  by  the  industrial  machine  designer.   This  "gap" 
exists  because  of  the  difficulty  involved  in  selecting 
a  single  mechanism  which  satisfies  all  of  the  various 
real-world  design  requirements.   This  dissertation 
advocates  optimization  theory  as  the  proper  tool  for 
bridging  this  gap.   But  while  the  application  of 
optimization  theory  to  planar  mechanism  design  has  risen 
at  a  seemingly  exponential  rate,  there  has  been  little 
application  to  spatial  mechanism  design.   The  author 
believes  the  reason  is  the  relatively  greater  complexity 
and  correspondingly  greater  number  of  parameters 
typically  encountered  in  spatial  mechanism  design. 
For  example,  only  six  design  variables  are  needed  to 
specify  the  planar  four-bar  function-generating 
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mechanism  discussed  in  section  2.2,  whereas  twenty-six 
such  design  variables  may  be  needed  to  completely 
specify  a  spatial  four-link  RCCC  motion-generating 
mechanism.    Even  when  a  problem  is  rather  poorly 
formulated,  optimization  can  often  be  performed  in  a 
reasonable  amount  of  time  when  dealing  with  six  or  ten 
design  variables.   However,  the  situation  becomes  decidedly 
different  when  dealing  with  twenty-six  design  variables. 
Since  each  variable  may  generally  take  on  an  infinite 

number  of  values,  a  twenty-six  parameter  problem 

2  Pi 

represents  on  the  order  of  °°    possible  solutions! 

However,  this  does  not  mean  that  the  application  of 
optimization  theory  to  problems  involving  a  large  number 
of  design  variables  should  be  avoided.   Quite  the  contrary, 
optimization  theory  offers  the  only  real  hope  of  solving 
such  difficult  problems.   The  point  is  that,  as  the 
number  of  design  variables  increases,  so  does  the  need 
for  an  efficient  problem  formulation  coupled  with  an 
efficient  optimization  method  and  a  means  for  reducing 
the  number  of  parameters.   Accordingly,  this  chapter 
attempts  to  define  the  goals  of  mechanism  optimization, 
and  to  develop  guidelines  for  obtaining  these  goals  in 


See  the  derivation  of  the  RC  and  CC  dyad  synthesis 
equations  in  Chapter  Four. 
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the  most  efficient  manner;  collectively,  this  will  be 
called  the  philosophy  of  mechanism  optimization. 

3. 2   Objectives  and  Constraints  of  Mechanism  Optimization 

There  are  typically  a  number  of  requirements  which 
must  be  taken  into  account  when  designing  a  mechanism. 
Some  of  these  are  listed  and  described  below. 

(1)  Motion  specification.   Most  often  a  mechanism 
is  required  to  generate    (a)  a  functional 
relationship  between  input  and  output  members; 
(b)  a  point  path,  sometimes  coordinated  with 
input  motion,  or   (c)  a  rigid-body  motion. 
Other  motion  requirements  are  possible,  but 
these  three  have  proven  to  be  adequate  for 
the  vast  majority  of  problems. 

(2)  Branch  avoidance.   Most  mechanisms  can  be 
assembled  in  two  or  more  distinct  configura- 
tions while  keeping  the  links  connected  in 
the  same  order.   Each  distinct  configuration 
is  called  a  branch.   The  mechanism  will 
generally  operate  in  one  branch.   Dissassembly 
is  required  to  move  it  into  a  different  branch. 
Branch  avoidance  implies  the  design  of 
mechanisms  such  that  the  entire  motion 
specification  lies  on  one  branch  only. 
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(3)  Order  of  positions.   Precision  positions 
must  occur  in  the  prescribed  sequence  and 
sense.   For  example,  if  the  prescribed 
positions  were  given  in  the  order  1,2,3,4, 

a  mechanism  generating  these  same  positions, 
but  in  the  order  1,3,2,4,  would  be  unaccept- 
able. 

(4)  Grashof's  condition.   This  refers  to  the 
relative  rotatability  of  links  within  a 
mechanism.   Often  it  is  desirable  to  drive 
the  input  link  of  the  mechanism  from  a 
continuously  rotating  source.   Such  an  input 
is  called  a  crank. 

(5)  Transmission  characteristics.   This  refers 
to  the  effectiveness  of  the  mechanism  in 
transforming  work  at  the  input  to  work  at 
the  output.   For  some  three-  and  four-link 
mechanisms,  this  can  be  expressed  in  terms 
of  a  transmission  angle. 

(6)  Link-length  ratio  restrictions.   Quite 
often,  when  trying  to  optimize  the  trans- 
mission characteristics  of  a  mechanism,  some 


See  Strong  and  Waldron  (88)  for  a  complete  discussion 
of  the  order  problem. 
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of  the  links  within  the  mechanism  become 
relatively  much  longer  than  others.   This 
may  cause  manufacturing  difficulties  as  well 
as  other  problems.   To  avoid  this  effect, 
limits  may  be  placed  on  the  allowed  ratio 
of  link  lengths. 

(7)  Fixed-pivot  location  restrictions.   Fixed 
pivots  must  be  located  so  they  do  not  inter- 
fere with  other  components.   They  must  also 
be  placed  where  they  can  be  connected  to  the 
machine  frame. 

(8)  Workspace  restrictions.   The  operating  space, 
or  workspace,  of  the  mechanism  must  often  be 
restricted  to  avoid  interference  with  other 
components. 

(9)  Dynamic  and  elasto-dynamic  property  restrictions, 
Velocities,  accelerations,  etc.  and  link 
deflections  must  remain  within  prescribed 
limits.   This  catagory  also  includes  balancing 
requirements. 

(10)   Tolerance  and  clearance  effect  restrictions. 

The  errors  in  mechanism  output  due  to  tolerances 
and  clearances  must  be  held  within  prescribed 
limits.   Since  these  are  not  usually  deter- 
ministic quantities,  stochastic  methods  must 
often  be  employed. 
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Of  course,  all  of  these  requirements  will  not  apply  to 
every  problem.   Quite  often  the  designer  can  use  his 
experience  and  judgment  to  eliminate  several  of  these 
from  consideration.   For  example,  dynamic  properties 
would  probably  not  be  a  restrictive  factor  in  a  linkage 
designed  for  the  purpose  of  guiding  open  an  automobile 
hood.   The  reader  can  undoubtedly  think  of  many  other 
examples . 

Classical  mechanism  design  procedures  generally 
begin  by  considering  the  motion  specification.   Typically, 
this  results  in  a  set  of  precision  position  requirements 
for  function,  path  or  rigid-body  motion  generation. 
These  requirements  are  expressed  as  equality  constraints 
which  are  then  solved  to  yield  the  dimensions  of  one 
or  more  mechanisms  which,  at  least  mathematically, 
satisfy  the  precision  positions.   This  procedure, 
commonly  known  as  mechanism  synthesis,  is  fascinating 
because  of  its  unusual  combination  of  geometric  and 
mathematical  complexities  coupled  with  a  concrete 
physical  phenomenon.   Unfortunately,  this  fascination 
has  sometimes  led  to  the  feeling  in  academics  that 
precision  position  synthesis  is  tantamount  to  mechanism 
design.   This  is  not  true,  since  each  of  the  previously 
discussed  requirements,  and  perhaps  more,  must  be 
considered  when  designing  a  mechanism. 
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Precision  position  synthesis  is  so  routinely  used 
as  the  first  step  in  mechanism  design  that  most  designers 
do  not  stop  to  question  why.   Are  the  precision  conditions 
of  greater  importance  than  the  other  requirements  listed? 
Is  a  mechanism  with  branching  problems  more  acceptable 
than  a  mechanism  that  does  not  satisfy  the  prescribed 
motion?   The  answer  to  both  of  these  questions,  of  course, 
is  "no,"  since  a  mechanism  must  satisfy  all  of  the 
design  requirements  if  it  is  to  be  a  workable  solution. 
Why  then  should  mechanism  design  always  begin  with 
precision  position  synthesis?   Why  not  begin  designing 
a  mechanism  by  first  considering,  for  instance,  the 
Grashof  condition?   From  the  set  of  all  possible 
mechanisms  of  a  given  type,  the  subset  consisting  of 
only  those  of  the  desired  Grashof  type  could  be  selected. 
From  this  subset,  mechanisms  meeting  the  specified 
precision  conditions  could  then  be  selected,  yielding 
an  even  smaller  subset.   Following  this,  the  link 
length  ratio  requirements  could  be  applied  as  a  third 
requirement,  and  so  forth.   This  approach  seems  practical 
enough,  and  yet,  to  the  author's  knowledge,  it  has 
never  been  applied. 

The  reason  for  using  precision  position  synthesis 
as  the  first  step  in  mechanism  design  is  not  that  it 
is  the  most  important  consideration.   Rather,  it  is 
used  as  a  first  step  because/ among  all  the  requirements 
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listed  at  the  beginning  of  this  section,  only  the 
precision  position  specifications  are  generally  in  the 
form  of  equality  constraints.   Thus,  it  is  often  possible 
to  use  the  elimination  of  variables  technique  discussed 
in  section  2.4.2.   In  addition,  it  may  be  possible 
for  the  designer  to  choose  the  number  of  equality 
constraints  by  choosing  the  number  of  specified  precision 
positions.   This  gives  the  designer  control  over  the 
number  of  design  variables  he  will  have  to  work  with. 
Although  used  intuitively  for  years,  the  concept  of 
precision  position  specifications  to  reduce  the  number 
of  free-choice  variables  was  formally  advanced  by  Tesar 
(26) ,  Eschenbach  and  Tesar  (27) ,  Tesar  and  Spitznagel 
(51,52)  and  Sutherland  (63),  among  others.   Of  these 
works,  only  Sutherland's  addresses  the  possibility  of 
approximately  satisfying  additional  motion  requirements. 

Two  principal  drawbacks  may  exist  to  using  precision 
condition  equality  constraints  to  eliminate  variables. 
First,  as  discussed  in  section  2.4.2,  this  technique 
introduces  nonlinearities  which  tend  to  distort  the 
objective  function  and  the  remaining  constraints.   Thus, 
for  example,  using  a  three  precision-position  solution 
in  which  the  unknown  variables  are  nonlinear  and  must  be 
solved  for  numerically,  even  though  one  more  variable  is 
eliminated  in  the  latter  case.   The  second  drawback  is 
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that  the  designer  may  not  have  a  need  to  satisfy  any 
precision  conditions  at  all,  but  rather,  may  require 
some  approximate  motion  specification  to  be  met.   In  this 
case,  a  great  many  possible  solutions  are  lost  by 
arbitrarily  setting  up  precision  position  requirements. 
If  no  acceptable  design  can  be  found  from  this  reduced 
solution  set,  the  designer  may  be  forced  to  eliminate 
some,  or  all,  of  the  precision  conditions. 

Although  not  so  obvious  or  as  well  developed  in  the 
literature,  it  may  sometimes  be  possible  to  use  other 
design  requirements  as  equality  constraints  to  eliminate 
variables.   For  example,  one  variable  could  be  eliminated 
by  specifying  the  coupler  link  to  be  twice  as  long  as  the 
input  crank.   Or  perhaps  the  minimum  and  maximum  trans- 
mission angles  could  be  specified  to  have  a  certain 
numerical  value.   Again,  the  designer  must  be  careful 
not  to  introduce  extreme  nonlinearities  in  the  remaining 
design  equations.   Although  not  explored  in  detail  here, 
this  appears  to  be  a  novel  and  potentially  useful  area 
of  research. 

Eschenbach  and  Tesar  (27)  have  suggested  dividing 
the  design  requirements  into  two  groups :   necessary 
requirements  and  desirable  requirements.   In  their  work, 
the  necessary  requirements  are  the  branching,  Grashof 
and  order  conditions,  plus  the  requirement  of  satisfying 
four  precision  conditions.   The  remaining  requirements 
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are  treated  as  desirable  conditions.   The  necessary 
requirements  are  considered  to  be  go  -  no  go  conditions, 
because  they  are  either  acceptable  or  they  are  unaccept- 
able.  In  other  words,  these  are  the  constraint  conditions. 
The  desirable  conditions  collectively  correspond  to  the 
objective  function,  the  idea  being  to  find  the  most 
desirable  value  for  a  given  combination  of  these  conditions. 

It  is  doubtful  that  the  division  of  necessary  and 
desirable  conditions  suggested  by  Eschenbach  and  Tesar 
(27)  will  apply  to  every  mechanism  design  problem.   For 
example,  satisfying  four  precision  positions  may  well 
be  a  desirable  condition,  while  fixed-pivot  location  may 
be  a  necessary  condition.   However,  their  work  is 
indicative  of  the  importance  placed  on  satisfying  the 
first  four  conditions  listed  at  the  beginning  of  this 
section,  namely,  motion  specification,  branch  avoidance, 
order  of  positions  and  Grashof's  condition. 

3. 3   Observations  and  Trends  Affecting  Mechanism 
Optimization 

The  development  of  a  general  philosophy  for  mechanism 
optimization  undoubtedly  must  involve  some  amount  of  sub- 
jective opinion  and  speculation.   Some  of  the  concepts 
presented  in  this  section  break  away  from  established 
trends  in  mechanism  optimization.   This  is  partially 
due  to  a  reevaluation  of  mechanism  optimization  procedures 
based  on  past  literature,  and  partially  due  to  the  rapidly 
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changing  role  of  the  computer  in  engineering  design. 
Although  it  is  difficult  to  make  broad  generalizations 
about  any  subject,  the  following  observations  seem  to 
apply: 

(1)   The  need  to  find  a  globally  optimum  solution 
has,  at  times,  been  overplayed  in  the 
literature  on  mechanism  optimization.   Most, 
if  not  all,  practical  mechanism  design 
problems  are  of  such  a  complex  and  multi- 
faceted  nature  that  it  is  impossible  to 
precisely  define  what  is  meant  by  optimum. 
Furthermore,  finding  the  global  optimum  is 
not  really  necessary  for  many  practical 
mechanism  design  problems.   Therefore,  the 
objective  of  this  dissertation  is  to  use 
optimization  theory  to  find  designs  which 
are  workable  solutions,  rather  than  emphasizing 
finding  the  optimum  solution.   For  this  reason, 
care  has  been  used  to  call  the  present  work, 
"Optimization  of  Spatial  Mechanisms," 
rather  than,  for  example,  "Optimum  Design  of 
Spatial  Mechanisms." 
(2)   Centralized  large-scale  digital  computers 
are  being  replaced  in  many  applications  by 
smaller  and  more  local  ones.   While  the 
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capabilities  of  these  smaller  machines 
have  increased  enormously  in  recent  years, 
storage  capacity  often  limits  their  useful- 
ness to  smaller-scale  problems.   The  trend, 
therefore,  should  be  toward  more  compact 
iteratively-based  programs  which  require 
less  storage. 

(3)  Since,  designers  will  more  often  be  using 
local,  and  perhaps  personal,  computers, 
the  algorithms  developed  should  be  easy  to 
program  and  readily  adaptable  to  a  variety 
of  problems.   This  means  that  the  design 
theories  should,  so  far  as  possible,  be 
based  on  more  easily  understood  concepts. 

(4)  Based  on  points  (2)  and  (3)  above,  the 
optimization  method  employed  should  be 
simple  in  concept,  require  little  computer 
storage,  and  should  be  capable  of  solving  a 
wide  variety  of  problems.   The  Hooke  and 
Jeeves'  method  described  in  Chapter  Two 
possesses  these  qualities. 

3. 3   Development  of  a  General  Mechanism  Optimization 
Philosophy 

In  the  past,  many  authors  have  developed  optimization 
methods  which  do  not  include  a  parameter  reduction  step, 
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for  example  (25,54).   A  simple  flow  chart  for  this 
method  is  shown  in  figure  3.1.   These  methods  might 
be  called  design  by  analysis,  because  there  is  no 
synthesis  step  involved.   The  advantage  of  this  approach 
is  that  no   potential  solutions  are  lost  through  parameter 
reduction.   Given  enough  time  and  computer  resources, 
this  would  be  the  preferred  method.   However,  this 
approach  will  generally   be  too  inefficient  for  use  on 
complex  mechanism  design  problems  where  a  large  number 
of  variables  are  present,  such  as  spatial  mechanism 
design. 

Undoubtedly,  the  most  popular  approach  to  mechanism 
optimization  has  been  the  "standard"  precision  position 
approach  shown  in  figure  3.2.   This  method  is  generally 
more  efficient  than  the  design  by  analysis  technique, 
because  the  number  of  free-choice  parameters  has  been 
reduced.   However,  this  method  is  still  lacking  in  the 
respect  that  no  attempt  has  been  made  to  distinguish 
nonparametric  constraints  from  parametric  ones.   In 
addition,  this  method  does  not  recognize  the  possibility 
of  eliminating  variables  other  than  by  precision  position 
synthesis.   Also,  no  provision  is  generally  included  to 
allow  additional  approximate  motion  specifications. 

A  flowchart  depicting  the  method  of  sequential 
filters,  as  presented  by  Spitznagel  and  Tesar  (51,52), 
is  shown  in  figure  3.3.   It  is  an  improvement  over  the 
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figure    3.1      Design  by  analysis 
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figure  3 . 2   Standard  approach  to  mechanism  optimization 
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figure  3.3   The  method  of  sequential  filters 
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figure  3.4   The  mechanism  optimization  method 
used  in  this  dissertation 


NOTE:   Feedback  may  occur  to  any  previous  block  in  all  figures. 
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standard  approach  because  it  recognizes  many  of  the 
special  properties  which  result  when  dealing  with  a 
planar  four-bar  linkage.   However,  no  distinction  is 
made  between  parametric  and  nonparametric  constraints. 
Also,  only  precision  position  synthesis  is  considered  as 
a  parameter  reduction  tool  and  no  provision  is  made  for 
satisfying  additional  approximate  motion  specifications. 

Sutherland  (63)  was  the  first  to  recognize  the 
potential  need  for  and  the  usefulness  of  satisfying  both 
precision  and  approximate  motion  specifications  simul- 
taneously.  Although  not  a  complete  optimization  scheme 
in  itself,  this  feature  is  both  a  desirable  and  a 
necessary  component  of  any  general  optimization  scheme. 

The  most  general  and  efficient  mechanism  optimization 
scheme  should  utilize  the  best  features  from  all  of  the 
above  methods.   It  should  incorporate  parameter  reduction, 
either  using  precision  position  synthesis  or  perhaps 
using  other  equality  constraints.   It  should  also  incor- 
porate approximate  motion  specification  in  addition  to 
precision  conditions.   Finally,  it  should  make  a 
distinction  between  parametric  constraints  and  nonpara- 
metric ones,  and  should  use  special  properties  of  the 
mechanism  to  eliminate  parametric  constraints,  when 
possible.   One  such  method  is  illustrated  by  the  flow- 
chart, figure  3.4.   While  simple  (and  perhaps  obvious) 
in  concept,  application  of  the  above-described  method 
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will  often  require  a  deep  understanding  of  the  specific 
problem  at  hand.   Spatial  mechanism  design,  the  topic 
of  the  next  four  chapters,  demonstrates  this  point. 
Chapter  Four  discusses  parameter  reduction  by  means  of 
closed-form  precision  position  solutions  for  various 
spatial  dyads.   Chapters  Five  and  Six  demonstrate 
methods  for  formulating  some  of  the  constraints  for  the 
RCCC  and  the  RSSR  mechanisms  in  nonparametric  form. 
These  are  then  used,  along  with  the  precision  position 
synthesis  methods  of  Chapter  Four,  to  demonstrate  the 
optimization  of  these  mechanisms. 


CHAPTER  4 
PRECISION  POSITION  SYNTHESIS  OF  SPATIAL  MECHANISMS 

4 . 1   Introduction  to  Precision  Position  Synthesis 

The  preceding  chapters  emphasized  the  need  for 
parameter  reduction  in  mechanism  optimization.   Clearly, 
precision  position  synthesis  is  often  of  importance,  and 
offers  a  readily  available  means  for  reducing  the  number 
of  parameters. 

In  this  chapter,  closed-form  solutions  are  developed 
for  rigid-body  guidance  problems  using  the  RS ,  CS ,  CC  and 
RC  dyads.   Since  the  function  generation  problem  can 
generally  be  converted  to  a  rigid-body  guidance  problem 
by  a  process  known  as  inversion  (19) ,  and  since  the  path 
generation  problem  can  generally  be  treated  as  an 
incompletely  specified  rigid-body  guidance  problem,  the 
solutions  developed  in  this  chapter  have  a  broad  spectrum 
of  applications.   Furthermore,  as  shown  in  the  following 
section,  the  dyads  treated  here  can  be  assembled  into  a 
variety  of  mechanisms,  giving  the  procedures  an  even 
wider  applicability. 

Synthesis  procedures  for  some  of  the  dyads  treated 
in  this  chapter  have  been  discussed  elsewhere  (76,102, 
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110,  111,  112).   However,  these  works  generally  do  not 
emphasize  finding  closed-form  linear  solutions,  a  topic 
of  considerable  importance  when  the  synthesis  process  is 
to  be  used  in  an  optimization  loop.   A  detailed  discussion 
of  these  synthesis  procedures  was  also  felt  to  be  needed 
in  order  to  provide  a  uniform  notation  for  later  reference 
and  to  provide  a  better  understanding  of  the  available 
free-choice  parameters.   All  of  the  synthesis  procedures 
given  here  are  for  the  maximum  number  of  precision 
positions  which  result  in  closed-form  linear  solutions. 

4 . 2   Dyadic  Synthesis  of  Mechanisms 

For  the  purposes  of  this  work,  dyads  may  be  thought 
of  as  building  blocks  for  mechanism  synthesis.   A  dyad 
is  a  two-link  kinematic  chain  composed  of  a  grounded 
link  and  a  floating  link,  and  having  two  degrees  of 
freedom.   The  grounded  link  is  joined  to  ground  through 
one  kinematic  pair  and  joined  to  the  floating  link 
through  another.   For  simplicity  and  clarity,  the 
concept  of  dyadic  mechanism  synthesis  will  be  explained 
using  the  planar  dyad  shown  in  figure  4.1.   Extension  of 
these  concepts  to  spatial  dyads  will  readily  follow. 

Begin  by  considering  a  set  of  discrete  planar 
precision  positions,  defined  by  f (x  .  ,  v.),  j  =  l,2,3,  as 
shown  in  figure  4.1.   Any  dyad  whose  tracer  point  can 
physically  reach  all  of  these  positions  can  generate 
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figure  4.1   A  planar  dyad 
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these  points.   In  fact,  it  is  possible  to  generate  an 
infinite  number  of  such  positions.   If,  however,  rotations 
of  the  floating  link  or  of  the  grounded  link  are  specified 
as  9  •  ,  j  =1 ,2,3  at  each  position,  the  problem  becomes  much 
more  difficult —  so  much  so,  in  fact,  that  it  is  now 
possible  to  generate  a  maximum  of  only  five  such  positions. 

When  the  rotations  of  the  floating  link  are  specified, 
the  problem  becomes  identical  to  the  planar  rigid-body 
guidance  problem.   Classical  Burmester  theory  is  a  well- 
known  method  for  solving  this  problem.   If  two  different 
dyads  can  be  found  which  satisfy  the  given  motion 
requirements,  their  floating  links  may  be  rigidly 
connected  to  form  a  constrained  planar  four-bar  linkage, 
as  shown  in  figure  4.2. 

It  should  be  apparent  that  spatial  mechanisms  can 
be  constructed  from  spatial  dyads  in  a  similar  fashion. 
Of  course ,  the  problems  become  more  complex  since  the 
motions  are  now  in  three  dimensions,  and  because  a 
greater  number  of  joint  types,  or  pairs,  must  be  consi- 
dered.  For  a  discussion  of  the  various  types  of  joints 
used  in  spatial  mechanisms  see  Harrisberger  (17) . 

When  designing  single-input  mechanisms,  dyads 
usually  must  be  connected  in  combinations  which  result 
in  a  single  degree  of  freedom,  although  a  number  of 
exceptions  to  this  rule  exist  (see,  for  example,  Shigley 
and  Uicker  (113)).   In  any  case,  several  single-input 
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:igure  4.2  Two  dyads  with  identical  floating  link 
rotations  joined  to  produce  a  four-bar 
linkage 
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spatial  mechanisms  which  can  be  synthesized  from  the 
dyads  treated  in  this  chapter  are  shown  in  figure  4.3. 
Thus,  for  example,  if  the  RC  and  CC  dyads  can  be  inde- 
pendently synthesized  for  three  rigid-body  positions, 
they  can  be  joined  together  to  form  an  RCCC  mechanism 
also  capable  of  generating  these  three  positions. 

4 . 3   Position  and  Orientation  of  a  Body  in  Space 

One  of  the  fundamental  tools  necessary  for  designing 
spatial  mechanisms  is  the  ability  to  describe  the  motion 
of  rigid  bodies  in  space.   Since  the  examples  presented 
in  this  work  deal  only  with  finitely  separated  positions, 
the  present  discussion  is  concerned  with  describing 
finite  displacements;  extension  of  this  concept  to  include 
higher  order  properties  can  be  found  in  many  standard 
texts,  for  example  (76). 

The  total  displacement  of  a  rigid  body  in  space  can 
always  be  considered  to  be  the  sum  of  an  angular  rotation 
and  a  linear  displacement  of  a  reference  point  fixed  in 
the  moving  body.   The  displacement  of  a  point  in  space 
is  easily  described  by  a  single  three-component  vector. 
However,  describing  angular  displacements  in  space  is 
not  so  easy  nor  is  the  method  so  obvious.   Some  of  the 
most  popular  methods  are    (1)  Euler  angles;  (2)  angular 
rotation  about  an  axis  in  space;  (3)  prescribed  order 
of  rotations  about  a  right-hand  set  of  Cartesian  axes; 
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"igure  4.3   Some  of  the  dyad-based  spatial  mechanisms 
which  can  be  synthesized  using  the  methods 
of  this  chapter 
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and  (4)  the  direction  cosines  of  a  pair  of  independent 
unit  vectors  fixed  in  the  body.   Simultaneous  translation 
and  rotation  can  be  specified  in  a  single  quantity  in 
several  ways,  such  as,  for  example,  (1)  quaternions  (96); 

(2)  3x3  matrices  and  tensors  with  dual  numbers  (98) ; 

(3)  4x4  matrices  using  homogeneous  coordinates  (101) ; 
and  (4)  dual  quaternions  (97),  etc.   Coaxial  translations 
and  rotations  (so-called  screw  displacements)  can  also 

be  expressed  by  way  of  any  of  these  methods. 

Regardless  of  the  specific  method  being  used,  it 
is  important  to  realize  that  only  three  independent 
parameters  are  needed  to  describe  spatial  angular 
displacements.   Nevertheless,  it  is  almost  always 
most  convenient  to  work  with  a  nine-component  three-by- 
three  matrix  when  describing  spatial  angular  displace- 
ments; this  is  the  so-called  rotation  matrix.   For  a 
given  angular  displacement  of  a  body,  all  of  the  methods 
listed  above  will  lead  to  rotation  matrices.   For  this 
reason,  the  rotation  matrix  will  be  the  basic  tool  used 
in  this  dissertation  to  describe  spatial  angular  displace- 
ments.  For  detailed  discussions  on  finding  the  rotation 
matrix  from  the  angular  displacement  description  methods 
listed  above,  see  (76,114). 

The  rotation  matrix  is  extremely  convenient  for 
specifying  finite  rotations  in  space.   For  example,  if 
v,  and  v?  represent  two  arbitrarily  oriented  spatial 


82 


positions  of  the  vector  v,  then  the  rotation  from 
position  1  to  position  2  can  be  expressed  in  the  form 


y2  =  [R]  y1 


(4.1) 


where  [R]  is  a  rotation  matrix  (76).   Equation  4.1  can 
be  expanded  in  terms  of  components  to  give 
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(4.2) 


4. 4   Synthesis  of  the  Revolute-Spheric  (RS)  Dyad 

The  revolute-spheric  dyad  (RS  dyad  for  short)  is 
shown  schematically  in  the  jth  position  with  associated 
vectors  in  figure  4.4.   It  is  one  of  the  simplest  dyads 
to  synthesize  and  is  perhaps  the  most  useful.   The  RS 
dyad  can  be  synthesized  for  up  to  three  precision  positions 
of  rigid-body  guidance  using  linear,  closed-form  solution 
procedures  (103) . 

The  following  vectors  are  used  in  figure  4.4  to 
describe  the  RS  dyad: 


a^    locates  the  revolute  joint  center  relative  to 

the  fixed  coordinate  system 
A.    locates  the  spheric  joint  center  in  the  jth 
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the  moving  body 
in  the  jth  position 


figure  4.4   The  RS  dyad  and  associated  vectors  in  the 
jth  position 
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position  relative  to  the  fixed  coordinate 

system 
o.    locates  the  center  of  the  moving  coordinate 

system  relative  to  the  fixed  coordinate  system 
r.    locates  the  spheric  joint  center  in  the  jth 

position  relative  to  the  moving  coordinate 

system 
s,    unit  vector  (denoted  by  the  hat)  along  the 

fixed  revolute  axis 

The  synthesis  procedure  can  be  viewed  intuitively 
as  follows.    Select  an  arbitrary  point  fixed  in  the 
moving  rigid  body  as  the  location  of  the  spheric  joint. 
The  location  of  this  point  at  the  three  prescribed 
positions  defines  a  plane  (three  points  define  a  plane) , 
and  a  circle  lying  in  this  plane  (three  points  also  define 
a  circle) .   Since  the  revolute  joint  physically  constrains 
points  to  lie  on  a  circle,  its  axis  must  pass  through  the 
center  of  this  circle,  and  must  also  be  along  the  normal 
to  the  plane  containing  the  three  points.   The  mathemat- 
ical formulation  of  this  easily  visualized  procedure  is 
developed  below. 

The  positions  of  a  body  in  space  are  specified  by 
giving  the  location  and  orientation  of  the  moving  coordi- 
nate system  embedded  in  the  body  at  each  position,  as 
described  in  section  4.3.   Thus,  for  three  prescribed 
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positions,  the  given  quantities  are 


oj  and  [R  ]      j  =  1,2,3  (4.3) 


where  o .  is  a  three  component  vector,  and  [R.]  is  the 
rotation  matrix  rotating  the  moving  body  from  a  reference 
position  to  position  j.   Displaced,  rather  than  absolute, 
positions  of  the  moving  body  are  of  importance  here. 
Therefore,  the  starting  position,  defined  by  o,  and  [R, ] , 
may  be  selected  arbitrarily,  and  the  other  two  positions 
measured  relative  to  the  first.   Choosing  [R,  ]  to  be  a 
3x3  identity  matrix  has  the  effect  of  making  the  fixed 
and  moving  coordinate  systems  parallel  in  the  initial 
position.   Since  this  assumption  generally  simplifies 
the  design  equation,  it  will  be  employed  throughout 
the  remainder  of  this  chapter.   The  initial  position 
vector,  o, ,  may  be  taken  to  be  any  convenient  point  in 
the  moving  body. 

Referring  to  figure  4.4,  begin  the  RS  dyad  synthesis 
by  assuming  a  value  of  r, ,  the  vector  locating  the 
spheric  joint  relative  to  the  origin  of  the  moving 
coordinate  system  in  the  initial  position.   Since  the 
locations  of  the  moving  origin,  o.,  and  rotation  matrices, 
[R.],  j  =  1,2,3  have  been  specified,  the  vector  A., 
locating  the  spheric  joint  center,  can  be  expressed  as 


A.  =  o .  +  r . ,      j  =  1,2,3  (4.4a) 

~d   ~:   ~: 
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Aj  =  o j  +  [Rj]r1#      j  =  1,2,3  (4.4b) 


It  has  already  been  noted  that  the  three  positions 
of  the  spheric  joint  define  a  plane  whose  normal  is 
parallel  to  the  fixed  revolute  axis  s, .   Since  the 
vectors  A_-A  and  A  -A_  lie  in  this  plane,  the  unit 
vector  s,  along  the  revolute  axis  is  defined  by 


(A2  -  A,)  x  (A3  -  A2) 

s1  =  — : : : —  (4.5) 

I  (A2  -  A1)  x  (A   -  A2)  i 


The  plane  containing  the  vectors  A- -A.  and  A  -A_ 
and  normal  to  s,  is  shown  in  figure  4.5.   It  is  necessary 
to  determine  the  vector  an  locating  the  intersection  of 
this  plane  with  the  fixed  revolute  axis.   This  is 
accomplished  by  first  determining  the  unit  vectors  p_ 
and  p_  which  are  respectively  perpendicular  to  the  vectors 
A--A,  and  A^-A„ 

p,  =  S,  x  (A-  -  A,)/|A,  -  A,  |  (4.6) 


P3  =  s,  x  (A3  -  A2)/|A3  -  A2 |  (4.7) 


The  perpendicular  bisectors  of  the  vectors  A_-A, 
and  A_-A_  intersect  at  the  tip  of  position  vector  aQ 
(as  shown  in  figure  4.5).   Denoting  the  perpendicular 
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figure  4.5   The  plane  of  the  RS  dyad 


distances  from  A„-A,  and  A-.-A_  to  the  tip  of  a   by  A- 
and  A,,  the  vector  a„  can  be  expressed  in  the  alternate 
forms 


?0  =  X2?2  +   h    + 


h      h 


(4.8) 


?0  =  A3?3  +  ^2  + 


~3  "  ~2 


:4.9) 


Equating  the  right  sides  of  equations  4.8  and  4.9  and 
forming  the  vector  product  with  p,-,  eliminates  A   and 
yields 


X3?2  X  ?3  +  ?2  x  ~2  +  ?2  x 


~3  "  .2 


(4.10) 


Equation  4.10  can  now  be  solved  explicitly  for  A.,  by 
forming  the  scalar  product  with  (p-  x  p..)  which  gives 


A3  = 


(p2  x  p3)  •  (C±    +  C2  -  C3  -  C4) 
(p2  x  p3)  •  (p2  x  p3) 


(4.11) 


where 


Si  =  ?2  X  bl 


(4.12) 


?2  =  ?2  X 


~2    bl 


(4.13) 
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C3  =  P2  x  A2  (4.14) 


~3  ~  A2 
C.  =  p9  x  —  (4.15) 


The  value  of  X-    thus  obtained  can  be  back  substi- 
tuted into  equation  4.9  to  determine  a_.   The  starting 
position  of  the  grounded  link,  represented  by  the  vector 
a, ,  can  be  obtained  from  the  expression 


tl   =   Al  "  ~0  (4.16) 


This  completely  determines  the  dimensions  and  the 
starting  position  of  the  RS  dyad. 

4. 5   Synthesis  of  the  Cylindric-Spheric  (CS)  Dyad 

The  procedure  for  synthesizing  the  CS  dyad  is 
similar  to  the  RS  dyad  synthesis  procedure.   Additionally, 
the  motion  of  the  grounded  cylindric  joint  along  its  axis 
must  be  considered.   This  scalar  variable  is  denoted  by 
S.  in  figure  4.6,  which  shows  the  1st  and  jth  position  of 
the  CS  dyad. 

The  CS  dyad  synthesis  procedure  can  be  visualized 
as  follows.   Referring  to  figure  4.6,  once  again  assume 
the  vector  r,  locating  the  spheric  joint  relative  to 
the  origin  of  the  moving  xyz  coordinate  system  in  its 
initial  position.   Also,  assume  the  orientation  of  the 
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jth  position 


1st  position 


figure  4.6   CS  dyad  in  the  1st  and  the  jth  positi 


position 
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cylindric  joint  axis,  s,,  in  the  fixed  XYZ  coordinate 
system.   Since  the  spheric  joint  is  physically  constrained 
to  lie  on  a  cylinder  about  s, ,  projections  of  its 
location  onto  a  plane  normal  to  s,  must  lie  on  a  circle. 
The  procedure  is  thus  to  project  the  three  locations  of 
the  spheric  joint  onto  a  plane  normal  to  s, .   These  three 
points  define  a  circle  within  the  plane,  and  the  remainder 
of  the  procedure  becomes  the  same  as  for  the  RS  dyad. 
This  procedure  is  described  mathematically  below. 

As  before,  three  positions  of  the  moving  body  are 
specified 


Oj,  [Rj],      j  =  1,2,3  (4.17) 


Assuming  the  vector  r,  arbitrarily  again  leads  to  the 
result 


A.  =  o.  +  [R.]r,  ,      j  =  1,2,3  (4.18) 


The  orientation  of  the  cylindric  joint  axis  is  now 
assumed  arbitrarily  by  specifying  two  of  its  components. 
The  third  component  can  be  found  from  the  unit  vector 
identity 

2      2      2 
sf   +  s^   +  sf   =  1  (4.19) 

lx     ly     lz 

Figure  4.7  shows  the  vector  s,  in  true  length.   The 
relationship  between  the  vector  s,  and  the  scalar 


92 


figure  4.7   The  CS  dyad  showing  vector  s,  in  true  length 
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displacements  S2  and  S3  can  be  determined  from  this 
figure  to  be 

S2  =  §1  '  (^2  "  ^1}  (4.20) 

s3  =  ?1  *  (A3  "  Ax)  (4.21) 

At  this  point  it  is  worthwhile  to  point  out  that 
the  arbitrarily  prespecified  variables,  r,  and  s, ,  are 
by  no  means  the  only  possible  choices.   In  fact,  for 
three  precision  position  synthesis,  any  five  scalar 
variables  may  be  selected  arbitrarily.   However,  the 
variables  selected  here  are  believed  to  be  the  most 
convenient  to  use  in  solving  the  kinematic  synthesis 
problem. 

Some  argument  could  be  made  for  preselecting  S_  and 
S3  rather  than  s,  (recall  that  only  two  independent 
scalar  quantities  comprise  s, ) ,  and  solving  equations 
4.19  through  4.21  for  the  components  of  s, .   This  is 
beneficial  because  the  designer  has  a  "feeling"  for 
what  values  to  assign  to  S~  and  S,  (constraints  would 
probably  be  placed  on  their  minimum  and  maximum  values) , 
whereas  selection  of  s,  would  be  arbitrary.   The  benefit 
of  having  feasible  starting  values  for  mechanism 
optimization,  however,  might  be  outweighed  by  the  addi- 
tional computation  required  to  find  s,  at  each  iteration. 


94 


The  author's  opinion  is  that  these  effects  would  roughly 
balance,  although  no  test  of  this  has  been  made. 

Figure  4.8  shows  section  A-A  taken  from  figure  4.7. 
This  is  the  plane  defined  by  the  unit  normal  s, ,  and 
passing  through  the  point  A  .   The  points  A'  and  A'  are 
the  projections  of  points  A2  and  A^  cnto  this  plane. 
The  vectors  locating  these  points  are  given  as  follows : 


A^  =  A2  -  S2s1  (4.22) 


A'  =  A3  -  S3s  (4.23) 


It  should  be  noted  that  the  right-hand  sides  of  these 
equations  are  now  known,  and  that  A'  =  A..  . 

In  the  CS  dyad,  the  cylindric  joint  constrains  the 
spheric  joint  to  move  on  the  surface  of  a  cylinder 
whose  axis  is  defined  by  s, .   The  projections  of  points 
on  the  cylinder  onto  a  plane  normal  to  s,  will,  therefore, 
lie  on  a  circle.   Thus,  points  A'  A'  and  A'  define  a 
circle  with  its  center  on  the  line  defined  by  s, . 

The  problem  of  determining  the  initial  location  of 
the  cylindric  joint,  a_,  is  thus  identical  to  the  RS 
dyad  synthesis.   Equations  4.6-4.16  apply  directly  to 
the  CS  dyad  synthesis,  except  that  A,,  A_  and  A^  should 
be  replaced  everywhere  by  A',  A'  and  A'.   The  dimensions 
of  the  CS  dyad  are  thus  completely  determined. 
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A'  -A' 


figure  4.8   Section  A- A  from  figure  4.7 
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4. 6   Synthesis  of  the  Cylindric-Cylindric  (CC)  Dyad 

A  schematic  representation  of  the  CC  dyad  in  the 
jth  position  is  shown  in  figure  4.9.   The  following 
vector  and  scalar  parameters  are  used  in  figure  4 . 9 
to  describe  the  CC  dyad.   All  vectors  are  expressed 
relative  to  the  fixed  coordinate  system  unless  other- 
wise stated. 

a„    locates  the  grounded  cylindric  joint  in  its 

initial  position 
an  .   locates  the  grounded  cylindric  joint  in  its 

jth  position 
s,    unit  vector  along  the  fixed  axis  of  the 

grounded  cylindric  joint 
SI .   scalar  displacement  of  the  grounded  cylindric 

joint  along  its  axis 
s„  .   unit  vector  along  the  ungrounded  cylindric 

joint  axis  in  the  jth  position 
a, ~   unit  vector  along  the  common  normal  from  the 

fixed  to  the  moving  joint  axis 

a'.    locates  the  jth  position  of  the  intersection 
-3 

of  the  ungrounded  joint  axis  and  the  common 
normal  to  the  joint  axes 
a.    locates  the  ungrounded  cylindric  joint  in 
the  jth  position 


97 


the  moving  body  in  the  jth  position 


figure  4.9   The  CC  dyad  and  associated  vectors  in  the 
jth  position 
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S2  .   scalar  displacement  of  the  ungrounded  cylindric 

joint  along  its  axis 
r.    locates  the  ungrounded  cylindric  joint  relative 

to  the  moving  (body-fixed)  coordinate  system 
o.    locates  the  origin  of  the  moving  coordinate 

system 

As  before,  three  prescribed  rigid-body  positions  are 


given 


Oy     [Rj]  ,      j  =  1,2,3  (4.24! 


From  figure  4.9,  the  following  vector  relationships  can 
be  deduced : 


*j  =  °j  +  Ej  (4.25) 


r.    =    [RjJr1  (4.26) 


§2j  =  [Rj]§21  (4"27> 


a'.  =  a.  -  S2  .s-  ■  (4.28 


*0j  =  *0  +  slj?i  (4-29) 


It  should  be  noted  that  the  joint  displacements  Si.  and 

1 

S2  .  may  be  set  equal  to  zero  in  the  initial  position 
without  loss  of  generality. 
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It  should  be  pointed  out  that,  unlike  the  two 
previous  synthesis  cases,  it  is  quite  difficult  to 
visualize  the  CC   dyad  synthesis.   For  this  reason,  the 
approach  that  will  be  taken  is  to  solve  a  set  of 
simultaneous  equations  derived  by  considering  the 
physical  constraints  that  exist  between  the  links  of 
the  dyad.   These  constraint  equations  are  enumerated 
below.  ■ 

(1)   Plane  equations:   these  require  the  link 
defined  by  unit  vector  a, ~  to  be  the 
perpendicular  bisector  of  the  joint  axes. 


s^  •  (a!  -  aQj)  =0,      j  =  1,2,3      (4.30) 


s         •  (a'  -  a   )  =0,     j  =  1,2,3      (4.31) 
~zj     ~j    ~uj 


(2)   Constant  twist  equations:   express  the 

requirement  that  the  twist  angle  between 
the  moving  axis,  s_  .  ,  and  the  fixed  axis, 


s, ,  remains  constant. 


s„ .  •  s,  =  s2,  ■  s,,     j  =  2,3         (4.32) 


(3)   Constant  moment  equations:   express  the 

requirement  that  the  moment  of  vector  s_ . 

about  the  s,  axis  must  be  constant. 
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fl  '  (!j  ~  ?0j)  X  ?2j 


§1  '  (?1  "  ~0)  x  ?21'     ^  =  2'3        (4.33) 

Condition  (3)  above  could  be  replaced  by  a  constant 
link-length  equation  of  the  form 


(!j  -  *oj]   '    (!j  "  ?0j 


(?1  "  ?01}  '  (*1  "  *01}  '         j  =  2'3        (4.34) 

However,  equation  4.3  3  leads  to  an  easier  solution 

because  it  is  linear  in  the  components  of  the  vectors 

a.,  a.-.,  j  =  1,2,3,  whereas  equation  4.34  is  quadratic 

in  these  components. 

For  three  position  synthesis  of  the  CC  dyad, 

fourteen  unknown  parameters  exist:   a.,  a, ,  s, ,  s~, , 

Sl„,  Si-,,  S2_,  S2_  (recall  that  s,  and  sn,  are  unit 
2     3     2     3  -1      -21 

vectors  containing  two  independent  parameters  each) . 

Equations  4.30  through  4.33  represent  a  total  of  ten 

scalar  equations.   This  means  that  four  scalar  parameters 

may  be  chosen  arbitrarily.   The  solution  strategy  is  as 

follows : 

(1)   Assume  s, ,  then  solve  for  s_,  from  equation 

4.32,  noting  that  s„ .  =  [R.]s_,  and  that 
3       -2]      j  -21 

=  1. 
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(2)  Expand  the  remaining  eight  equations 
(equations  4.30,  4.31  and  4.33)  in  terms 

of  the  ten  unknown  quantities  (an,  a,,  Sl_  , 
Si,,  S2_,  S2.,)  by  substituting  for  other 
unknown  quantities  from  equations  4.25 
through  4.29. 

(3)  Next,  arbitrarily  assume  Sl~  and  Si-,,  and 
solve  the  resulting  set  of  eight  linear 
equations  in  eight  unknowns. 

The  details  of  this  procedure  are  given  below. 

Assume  the  orientation  of  the  fixed  cylindric  joint 

axis  s,  =  s,  i  +  s,  j  +  s,  k.   Rewriting  equations  4.32 
-1     lx~     ly^     lz~  ^   ^ 

with  j  =  2 , 3 


s-j^  =  s21  ■  s-^  (4.35; 


!23  ■  !l  =  !21  '  !l  (4'36) 


From  equation  4.27,  sL,,  and  s„.  can  be  expressed  in  terms 
of  s21 

s22  =  [R2]s2]_  (4.37) 


s23  =  [R3]s21  (4.38) 


where   [R-,]  and   [R-,]  are  3x3  rotation  matrices  of  the 
form 
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[R: 


all     a12 


a31      a32 


'13 


a21      a22      a23 


'33 


(4.39) 


[R3]  = 


bll      b12      b13 


b21      b22      b23 


31 


32 


33 


(4.40) 


Substituting  equation  4.37  into  4.35  and  substituting 
equation  4.38  into  4.36  gives 


[R2  J!2i  '  !i  "  !2i  '  !i 


(4.41) 


[R,  ]s 


3  J~21  '  ?1    ?21  '  ?1 


(4.42) 


Expanding  these  equations  into  components  using 

s~,  =  s_,  i  +  s-t  j  +  sn,  k   and  s,  =  S-,  i  +  s,  j  +  s,  k 
-21    21x~     21ys     21z~      ~1     lx~     ly~     lz~ 

and  substituting  for  [R_]  and  [R-,]  from  equations  4.39 
and  4.40  yields,  after  some  manipulation 


As-,   +  Bsni   +  Cs„,   =  0 
21x      21y      21z 


(4.43) 


Ds~,   +  Es0,   +  Fs_,   =  0 
21x      21y      21z 


(4.44) 


where 


A  =  (a,.  -  l)s,   +  a_,s,   +  a-,,s, 
11       lx     21  ly     31  Iz 


(4,45) 
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B  =  a12Slx  +  (a22  '  1)sly  +  a32Slz  (4'46) 


C  =  a,_,S-   +  a„,s-   +  (a,.,  -  l)s.  (4.47' 

13  lx     23  ly      33       lz 


D  =  (bll  ~  1)Slx  +  b21Sly  +  b313lz  (4'48) 


E  =  b123lx  +  (b22  "  1)sly  +  b32Slz  (4'49) 


F  =  b13Slx  =  b23Sly  +  {b33  "  1] Slz  (4'50) 


Since  the  a's  and  b's  are  prescribed  components  of  the 

rotation  matrices,  and  since  s,  =  s,  i  +  s,  j  +  s,  k  has 
'  ~1     lx-     ly~     lz~ 

been  assumed  arbitrarily,  the  values  of  the  scalar 
constants  A  through  F  may  be  calculated.   Therefore, 
equations  4.43  and  4.44  are  linear  in  the  unknown  com- 
ponents of  L,.   A  third  equation  is  obtained  using  the 
expression 


§21  •  !21  "  1  (4-51) 

or,  in  terms  of  components 

si*    +  s?,   +  si.       =1  (4.52) 

21x     21y     21z 


Solving  equations  4.43,  4.44  and  4.52  for  s2lz  <3ives 


j 

s     =±1 (4.53! 

21z      \  I  2  9 

EC  -  FB  DC  -  FA     +  1 

DB  -  EA  EA  -  DB 
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and  s_,   and  sn,   are  given  by 
21x       21y 


'21x 


EC  -  FB 


DB  -  EA 


'21: 


(4.54) 


2xy 


DC  -  FA 


EA  -  DB 


'21z 


(4.55) 


The  direction  of  the  moving  cylindric  joint  axis  in  its 

initial  position  is  now  completely  determined.   It  should 

be  noted  that  s0,   and  hence  s«,   and  s~1      will  have  two 
21z  21x       21y 

values,  corresponding  to  the  two  roots  of  equation  4.53. 
Either  one  of  these  two  roots  can  be  used  because  the 
joint  axis  is  determined  by  the  direction  of  the  vector 
and  its  sense  is  of  no  importance. 

The  direction  of  the  moving  joint  axis,  s_,  ,  has 
been  determined  using  the  assumed  direction  s   together 
with  the  rotation  matrices  [R  ]   and  [R  ]  .   This  is 
important  in  the  design  of  the  RCCC  mechanism  (see 
Chapter  Five)  since  Grashof  and  branching  conditions  for 
this  mechanism  depend  only  on  the  directions  of  the  joint 
axes . 

Returning  to  the  CC  dyad  synthesis,  the  following 
vector  relations  are  obtained  from  equations  4.25  through 
4.29. 


*1  =  Si  +  ll 


*2    =  °2  +  £2  =  °.2    +  [R2]£l 


(4.56) 
(4.57) 
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-3  =  ?3  +  £3  =  ?3  +  [R3]:i  (4.58) 

~1  =  ~1  =  ?1  +  El  (4.59) 

a2  =  a2  "  S22^2  =  °2  +  [R2^rl  ~  S22^R2^1      (4.60) 


a^  =  a3  -  S23s2  =  o3  +  [R-^]^  -  S23[R3]s1      (4.61) 


a01  =  !o  (4'62) 


?02  =  a0  +  S12!l  (4'63: 


?03  =  a0  +  Sl3?l  (4"64) 


Substituting  these  results  into  the  remaining  constraint 
equations  (equations  4.30,  4.31  and  4.32)  gives 


!l  '  (2l  +  £l  "  to]    =  °  (4.65) 


!l    •     (22    +    fR2^1    -    S22^R2^21    '    ?0    "    S12!l}    =    °  (4'66) 

!l    *     (?3    +    [R3^1   "    S23[R3^21   -    to   ~   S13!l>    =    °  (4'57) 


(ol  +  rl  "  ~0}    =  °  (4.68) 


^^l*     '     {°2    +    [R2]£l    "    S22[R2]!21    "    to    ~    S12!i}    -   °    (4'69) 
([R3]!21}     •     (?3    +    [R3^1    "    S22[R3]!21    "    to    ~    S13!l)    =    °    (4-70) 
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!i  •  (e2  +  tR2]5i-  S22tR2]!2i  -  e0  -sl2!i}  x 

[R2]!2i  =  !i  •  (2i +  5i  - 10]  x  !2i  (4-71) 

!i  '  (23  +  fR3^i  -  S23tR3^2i  -  e0  -  513h]  x 

[R3]!2i  =  !i  •  (2i  +  ^i  -  to]  x!2i  (4-72) 

If  Sl~  and  SI,  are  assumed  arbitrarily,  equations 
4.65  through  4.72  form  a  set  of  eight  scalar  equations 
in  the  eight  unknown  parameters  r, ,  an ,  32^  and  S2- 
(all  other  parameters  are  either  prescribed  or  have 
already  been  calculated) .   Expanding  these  equations 
using  the  rules  of  dot-product  and  vector  cross-product 
multiplication  will  show  them  to  be  linear  in  the 
unknown  quantities.   The  result  of  this  expansion  is 
shown  in  matrix  form  in  figure  4.10.   Although  this 
appears  on  the  surface  to  be  somewhat  complex,  it 
should  be  remembered  that  the  8x8  coefficient  matrix  and 
the  8x1  column  vector  on  the  right  of  the  equality  sign 
will  simply  be  arrays  of  numbers.   The  solution  vector 
can,  therefore,  be  easily  found  using  any  of  the  standard 
methods,  such  as  Gaussian  elimination  or  matrix  inversion, 

Upon  solving  the  above  derived  set  of  equations, 
the  dimensions  and  starting  position  of  the  CC   dyad  are 
completely  determined. 
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4. 7   Synthesis  of  the  Revolute-Cylindric  (RC)  Dyad 

Synthesis  of  the  RC  dyad  follows  simply  and  directly 
from  the  foregoing  synthesis  of  the  CC  dyad.   It  is 
apparent  that  the  RC  dyad  is  identical  to  the  CC  dyad 
with  one  degree  of  freedom  removed,  namely,  the  trans- 
lation along  the  grounded  cylindric  joint  axis. 
Referring  to  the  schematic  of  the  CC  dyad  shown  in 

figure  4.8,  it  is  seen  that  the  parameter  SI.  will  be 

D 

equal  to  zero  for  all  possible  positions  of  the  RC  dyad. 
Fortunately,  in  synthesizing  the  CC  dyad,  SI-  and  Sl_ 
are  free-choice  (arbitrarily  assumed)  parameters.   As 
a  result,  the  RC  dyad  synthesis  equations  are  the  same 
as  the  CC  dyad  synthesis  equations,  except  that  Si-  and 
Sl^  must  be  set  equal  to  zero.   Also  recall  that  SI, 
is,  by  definition,  equal  to  zero. 

It  is  important  to  realize  that  setting  Sl?  and  S1-. 
equal  to  zero  only  insures  zero  displacement  along  the 
cylindric  joint  axis  at  these  two  positions.   It  is, 
therefore,  necessary  to  physically  constrain  the  revolute 
joint  from  moving  along  its  axis. 

4 . 8   Conclusions  of  Precision  Position  Synthesis 

Methods  have  been  presented  for  synthesizing 
RS ,  CS ,  CC  and  RC  dyads  to  meet  three  prescribed 
precision  positions  of  spatial  rigid-body  guidance. 
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In  each  case,  three  is  the  maximum  number  of  such  posi- 
tions which  result  in  linear  closed-form  solutions. 

An  extremely  important  point  to  be  made  here  is 
that  the  use  of  these  closed-form  synthesis  techniques 
does  not  limit  the  number  of  prescribed  positions  to 
three.   For  example,  if  only  two  precision  positions  are 
specified,  the  third  position  can  be  considered  to  be 
a  set  of  six  variable  parameters  in  the  optimization 
scheme  (notice  that  six  parameters  will  thus  be 
reintroduced  in  the  design  equations) .   If  more  than 
three  precision  positions  are  specified,  the  requirement 
that  these  be  satisfied  can  be  incorporated  in  the 
objective  function  or  the  constraint  equations  of  the 
optimization  procedure.   Of  course,  only  a  limited 
number  of  precision  conditions  can  be  satisfied. 
Beyond  this  number,  only  approximate  solutions  can 
generally  be  found. 

The  advantages  of  the  vector-based  approach  presented 
here  include  generality  and  extendability  and,  of  course, 
the  widespread  familiarity  engineers  have  with  vector 
methods.   The  generality  and  extendability  of  this 
approach  can  be  seen  from  the  work  of  Hernandez  (115) , 
who  used  vector  constraint  equations  in  synthesizing 
spatial  mechanisms  containing  higher  pair  joints,  such 
as  sphere-plane  and  cylinder-plane  pair  joints. 
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The  next  two  chapters  demonstrate  the  use  of  the 
synthesis  techniques  developed  in  this  chapter  as  a 
parameter  reduction  tool  in  the  optimization  orocess. 


CHAPTER  5 
OPTIMIZATION  OF  THE  RCCC  MECHANISM 
5. 1   Problem  Definition 

Although  the  RCCC  mechanism  is  well  known  to 
kinematicians ,  few  examples  can  be  found  of  its  practical 
application.   Potentially,  it  is  useful  for  the  generation 
of  path  and  rigid-body  motions  and  also  for  the  generation 
of  variable-pitch  screws  (73) . 

A  great  deal  of  work  remains  to  be  done  before 
design  theories  for  the  RCCC  mechanism  are  complete  to 
the  point  where  all  the  design  requirements  discussed  in 
Chapter  Three  can  be  included.   The  example  of  this 
chapter  does,  however,  demonstrate  the  design  of  an  RCCC 
mechanism  to  meet  rigid-body  motion  requirements  while 
being  free  from  branch,  order,  link-length  ratio  and 
fixed  pivot  location  problems  and  having  complete  input 
crank  rotatability .   It  is  believed  that  the  solution  of 
this  problem  represents  a  significant  step  forward  in  the 
state-of-the-art  of  design  of  this  mechanism. 
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5-2   Satisfying  Additional  Motion  Requirements 

Sections  4.6  and  4.7  presented  linear  closed-form 
methods  for  synthesizing  CC  and  RC  dyads  to  satisfy 
three  positions  of  rigid-body  guidance.   Once  synthesized, 
these  dyads  may  be  joined  together  to  form  an  RCCC 
mechanism  capable  of  satisfying  the  three  specified 
positions.   Often,  however,  it  is  desirable  to  satisfy 
more  than  three  positions.   Unfortunately,  the  RC  dyad 
cannot  in  general  be  synthesized  for  more  than  three 
precision  positions  of  rigid-body  motion  generation  (76) . 
The  CC  dyad  can  be  synthesized  for  a  maximum  of  five 
exact  rigid-body  positions.   However,  the  synthesis 
equations  are  nonlinear  for  both  four  and  five  positions 
and,  as  a  result,  numerical  solution  procedures  are 
required  in  each  case  (76) . 

Although  more  than  three  exact  positions  of  rigid- 
body  guidance  will  not  generally  be  obtainable  with  the 
RCCC  mechanism,  a  dyad-based  procedure  will  now  be 
described  which  allows  additional  positions  to  be  satis- 
fied in  an  approximate  sense.   Since  the  RC  dyad 
synthesis  procedures  follow  directly  from  the  CC  dyad 
synthesis  procedures  (see  Chapter  Four) ,  the  method 
described  below  will  be  for  the  CC  dyad. 

Suppose  that  four  rigid-body  precision  positions 
have  been  specified  as 
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Oj/  [Rj]»      J  =  1,2,3,4  (5.1) 


This  four-position  motion  specification  is  now  broken 
into  two  three-position  motion  specifications 

Oj/  [RjL      j  =  1,2,3  (5.2) 


Oj,  [Rj J,      j  =  1,2,4  (5.3) 


Now  refer  to  the  three-position  synthesis  procedure 
for  the  CC  dyad  described  in  section  4.6.   The  direction 
of  the  grounded  joint  axis,  s,,  and  the  scalar  displace- 
ments Sl2  and  SI   are  assumed  arbitrarily.   Substituting 
these  assumed  values  into  the  synthesis  equations  yields 
the  complete  dimensions  and  starting  position  of  the  CC 
dyad.   Notice  that,  for  the  same  values  of  s  ,  Sl„  and 
Si.,,  the  motion  specifications  given  by  equations  5.2 
and  5.3  will  generally  result  in  different  solution 
dyads.   If,  however,  a  set  of  values  for  s,  ,  SI-  and  Si., 
can  be  found  which  yield  the  same  CC  dyad  from  both 
motion  specifications,  the  resulting  dyad  will  exactly 
generate  all  four  prescribed  positions  of  equation  5.1. 
Referring  to  figure  4.9,  it  can  be  seen  that  the  two 
dyads  will  be  the  same  provided  the  vectors  an,  s_,  and 
r,  are  the  same  for  each  dyad.   Using  the  superscripts 
5.2  and  5.3  to  denote  which  motion  specification  the 
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above  vectors  correspond  to,  the  constraint  equations 
for  four-position  synthesis  may  be  written  as 


a5-2     5-3 

*0   "  ?0   =  °  <5-4) 


?21   "  !21   =  °  (5-5) 


5.2     5.3 
fl    "  5l    =  °  (5.6) 


If  one  of  the  prescribed  positions  is  an  approximate 
rather  than  a  precision  motion  specification,  the  above 
equalities  need  not  strictly  apply.   Instead,  the  solution 
can  be  cast  in  the  form  of  a  function  whose  value  is  to  be 
minimized.   One  possibility  for  such  a  function  is 

F(!l'S12'S13>  "  Wll?o'2  "  ?o'3! 

+  W2l!2i2  "  !2i3|  +W3lfl'2  "  it3!  (5-7) 

where  the  absolute  value  signs  denote  the  magnitude  of 
the  vector,  i.e.,  |a|  =  (a  •  a)  %  and  where  W. ,  W   and 
W3  are  positive  weighting  factors  which  determine  the 
relative  importance  of  each  term. 

Clearly,  when  F  equals  zero  both  CC  dyads  are  the 
same,  and  all  the  prescribed  positions  are  satisfied 
exactly.   When  F  is  greater  than  zero,  the  dyads  are 
not  the  same,  and  only  three  of  the  four  positions  are 
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satisfied  exactly.   If  a^ ' 2 ,  s^2  and  r^2  are  taken  as 
the  solution,  positions  1,2  and  3  will  be  exactly 
satisfied  and  position  4  will  be  approximately  satisfied. 
On  the  other  hand,  if  a^3,  s^3  and  r^ "  3  are  taken  as 
the  solution,  positions  1,2  and  4  will  be  exactly  satisfied 
and  position  3  will  be  approximately  satisfied. 

The  synthesis  procedure  described  above  will 
hereafter  be  referred  to  as  the  method  of  position 
stacking.   As  far  as  the  author  is  aware,  no  published 
work  exists  on  the  use  of  this  method  for  the  partially- 
exact,  partially-approximate  synthesis  of  either  spatial 
or  planar  mechanisms.   The  method  of  position  stacking 
has  several  advantages.   First,  it  retains  the  closed- 
form  parameter  reduction,  which  greatly  reduces  the 
number  of  components  in  the  design  vector.   Second,  it 
can  easily  be  extended  to  include  any  number  of  positions. 
For  example,  nine  positions  could  be  specified  in  three 
groups  of  three.   A  third  advantage  is  that  the  maximum 
number  of  precision  positions  can  be  exactly  satisfied 
without  deriving  and  numerically  solving  a  set  of  non- 
linear equations.   A  fourth  advantage  of  this  method  is 
that  it  can  be  applied  to  any  dyad  for  which  a  closed- 
form  solution  is  available. 

The  position  stacking  method  differs  from  the 
partially-exact,  partially  approximate  synthesis  method 
proposed  by  Sutherland  (63)  in  several  respects.   First, 


116 


it  can  be  applied,  without  special  formulation,  to  any 
mechanism  for  which  a  closed-form  solution  is  available. 
Sutherland's  work  deals  exclusively  with  planar  four-bar 
linkages  (although' extension  to  other  mechanisms  is 
possible) .   A  second  difference  lies  in  the  method  of 
approximation.   Sutherland  determines  the  moving  pivot 
locations  which  result  in  the  smallest  least-squares 
deviation  in  the  length  of  the  grounded  link.   The 
position  stacking  method,  when  applied  to  a  planar  four- 
bar  linkage,  would  also  include  the  difference  in  the 
locations  of  the  fixed  pivot. 

The  least  squares  method  proposed  by  Sutherland 
does  offer  the  possibility  of  finding  an  optimal 
solution  in  closed  form  -  provided  the  synthesis  equations 
are  linear.   Such  a  solution  method  has  been  discussed 
in  greater  detail  by  Gmerek  and  Matthew  (116) .   Solutions 
to  the  position  stacking  problems  presented  in  this  work 
are  found  using  optimization  methods.   It  should  be 
possible  to  find  closed-form  solutions  using,  for  example, 
the  methods  presented  by  Gmerek  and  Matthew  (116),  since 
the  synthesis  equations  are  linear  in  the  unknown  quantities 
However,  such  a  solution  probably  would  not  satisfy  the 
constraint  conditions  and  hence  would  not  be  an  acceptable 
design. 

Mixed  sets  of  exact  and  approximate  positions  are 
believed  to  have  many  practical  applications  (63) .   The 
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common  transfer  operation,  where  an  object  is  picked  up 
at  one  position,  moved  past  an  intermediate  obstacle  and 
placed  at  another  position,  is  one  example.   The  pick 
and  place  positions  must  be  satisfied  exactly,  whereas 
intermediate  positions  only  need  to  be  satisfied 
approximately . 

The  principal  drawback  of  the  position  stacking 
method  is  that  the  designer  does  not  have  direct  control 
over  the  error  at  the  approximate  positions.   The 
positions  are  approximated  in  the  sense  that  the  dyad 
or  dyads  that  generate  the  approximate  positions  are 
similar  to  the  dyad  being  used,  i.e.,  the  one  generating 
the  exact  positions.   It  is  not  possible  to  evaluate  the 
position  error  until  the  dyads  have  been  assembled  into 
a  mechanism. 

The  effectiveness  of  the  position  stacking  method 
has  been  demonstrated  in  working  the  examples  of  this 
dissertation  and  also  in  verifying  the  numerical  example 
of  reference  (103) .   However,  more  detailed  research  is 
needed  to  determine  the  potential  usefulness  of  this 
novel  method. 

5. 3   The  Grashof  Condition 

Grashof  analysis,  also  known  as  type,  limit  position 
or  mobility  analysis,  attempts  to  determine  the  relative 
rotation  of  links  within  a  given  mechanism.   A  fully 
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rotating  ground-pivoted  link  is  called  a  crank.   A  link 
which  oscillates  about  a  fixed  pivot  between  limit  positions 
is  called  a  rocker.   The  present  analysis  seeks  to  deter- 
mine whether  or  not  the  input  link  of  an  RCCC  mechanism 
of  known  dimensions  is  a  crank. 

One  obvious  way  to  find  this  out  is  to  perform  a 
kinematic  analysis  of  the  relative  displacements  of  the 
links  within  the  mechanism  for  a  series  of  closely 
spaced  positions  of  the  input  link.   A  mobility  limit 
position  is  found  when  the  solution  of  the  loop-closure 

equation  becomes  imaginary.   Unfortunately,  the  approach 
just  described  involves  the  independent  input  of  the 

mechanism.   Any  optimization  based  on  this  approach  must 

necessarily  be  parametric. 

Fortunately,  the  works  of  Duffy  and  Gilmartin  (83-85) 

have  eliminated  the  need  for  involving  the  independent 

input  in  determining  the  limit  positions  of  the  RCCC 

mechanism.   What  follows  is  a  brief  summary  of  portions 

of  their  work. 

The  RCCC  mechanism  is  shown  in  the  notation  of 

Gilmartin  and  Duffy  (84)  in  figure  5.1.   The  parameters 

used  to  specify  it  are  as  follows: 

s.    the  unit  vector  along  the  i    pair  axis 
a. .   the  unit  vector  along  the  common  perpendicular 
between  s .  and  s . 
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figure  5.1   The  RCCC  mechanism  and  its  spherical  RRRR 
equivalent 
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aij  the  len<?th  of  the  common  perpendicular  which 
is  taken  as  the  kinematic  link  length 

S.  the  distance  between  links  a.  .  and  a.,  aloncr 
the  pair  axis  s .  (a  double  subscript  such  as 
S,-  denotes  an  invariant  distance) 

aii   the  P°sitive  angle  from  s.  to  s .  measured 

clockwise  looking  in  the  direction  of  a. ., 

i.e.,    in   the   positive    (or    right-hand)    direction 

about   a . . 
-ID 

®j_        the  positive  angle  from  a.  .  to  a.   measured 
clockwise  looking  in  the  direction  of  s . 

It  can  be  shown  (see  for  example  the  text  by  Duffy 
(117))  that  the  spatial  RCCC  has  an  equivalent  spherical 
RRRR  representation  in  which  all  angular  relationships 
within  the  mechanism  are  maintained.   This  equivalent 
spherical  mechanism  is  also  shown  in  figure  5.1.   Notice 
that  the  cu.'s  become  arcs  of  great  circles  physically 
representing  the  link  lengths.   Now,  since  angular 
relationships  are  maintained,  the  spherical  mechanism  and 
the  spatial  mechanism  must  be  of  the  same  type,  i.e., 
crank-rocker,  double-crank  or  double-rocker.   The  input 
link  rotatability  may  be  determined  by  first  replacing  all 
obtuse  link  lengths  by  their  supplements  so  that 
0°  1  ou.  <  90°  for  all  a.   Now,  the  input  link,  a,-,  will 
be  a  crank  if  it  is  the  shortest  link  and  if 

a12  +  90°  <_  h(a12    +   a23  +  a34  +  a41)  (5.8) 
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The  above  expression  provides  a  simple,  nonpar ame trie 
means  for  determining  the  input  link  rotatability  in 
the  RCCC  mechanism. 

5. 4   The  Branch-Avoidance  Condition 

The  mechanism  branching  problem  was  discussed 
briefly  in  section  3.2.   It  occurs  when  the  prescribed 
positions  lie  on  more  than  one  physical  closure  of  the 
mechanism.   This  problem  generally  renders  the  solution 
unsuitable  for  the  task  at  hand. 

Up  until  now,  the  only  way  to  avoid  this  problem 
in  designing  the  RCCC  mechanism  was  to  perform  a  position 
analysis  on  every  candidate  solution.   By  checking  a 
mechanism  at  a  large  number  of  positions  of  the  input, 
the  designer  could  see  that  the  mechanism  did  or  did  not 
pass  through  the  desired  positions.   This  analysis 
necessarily  involved  the  independent  input  and  was, 
therefore,  parametric  optimization. 

An  easy-to-apply  nonparametric  method  for  determining 
the  branching  characteristics  of  the  RCCC  mechanism  will 
now  be  developed.   Begin  by  considering  the  two  branches 
of  the  planar  four-bar  mechanism  as  shown  in  figure  5.2. 
It  is  evident  that,  for  the  upper  branch 


-23  X  ?34  '  ?3  <  °  (5-9 
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figure  5 . 2   The  two  branches  of  the  planar  four-bar 
linkage 
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and,  for  the  lower  branch 

2     2         n 
-23  X  a34  '  !3  >  °  (5.10) 

where  s^    is  the  joint  axis  unit  vector  pointing  out  of 
the  plane  of  the  paper  and  where  the  superscripts  1  and 
2  refer  to  the  first  and  second  position  of  the 
mechanism.   The  mechanism  can  only  move  from  one  branch 
to  the  other  by  passing  through  a  special  configuration, 
where  the  link  vectors  a23  and  a^.    are  collinear.   The 
above  logic  leads  to  the  conclusion  that,  either 

?23  X  ?34  *  §3  K    °      for  a11  3  (5.11) 

or 

-23  x  a34  *  §3  >  °      for  a11  3  (5.12) 

if  the  positions  j  are  to  lie  on  the  same  branch. 

Now  returning  to  the  spherical  mechanism  of  figure 
5.1,  a  similar  condition  can  be  seen  to  apply.   The 
vectors  a23  and  a34  represent  chords  of  the  great  circles 
from  ?2  to  ^3  and  from  ^3  to  s4,  respectively.   The  two 
branches  of  the  mechanism  can  be  distinguished  by  the 
direction  of  the  vector  product  a„3  x  a...   The  resulting 
vector  will  always  be  directed  toward  the  inside  of  the 
sphere  for  one  branch  and  toward  the  outside  of  the 
sphere  for  the  other  branch.   A  special  configuration 
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is  reached  when  the  resultant  vector  is  tangent  to  the 
sphere  at  joint  axis  s3-   Since  s3  is  always  an  outward 
pointing  normal  to  the  sphere,  the  branch  avoidance 
conditions  for  the  spherical  mechanism  are  the  same  as 
for  the  planar  four-bar,  i.e.,  equations  5.11  and  5.12 
apply  directly. 

It  is  convenient  to  express  this  result  in  terms 
of  the  joint  axes  unit  vectors  only.   The  vector 
difference  s3  -  s2  produces  a  vector  in  the  direction 
of  ?2  3'  althou9h  not  necessarily  of  the  same  magnitude 
(see  figure  5.3).   However,  since  only  the  directions 
of  the  vectors  a12  and  a23  are  important  in  equations 
5.11  and  5.12,  the  RRRR  spherical  mechanism  will  be  free 
from  branching  problems  if,  in  all  prescribed  positions, 
either 

(?3  "  !2}  X  (?4  "  h]     *  !3  =  U    X  §2  '  !3  <  °  (5'13) 


(!3  ~  ?2]  X  (?4  "  !3)  '  §3  =  !4  x  !2  '  ?3  >  °  (5.14) 

This  constraint  condition  applies  directly  to  the  RCCC 
mechanism  since  it  is  an  angular-motion  equivalent  to 
the  spherical  RRRR  in  terms  of  relative  rotations  at 
each  joint.   Equations  5.13  and  5.14  are  nonparametric 
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?23    ?3    ?2 


figure  5.3   The  vector  chord  of  the  spherical  mechanisi 
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constraints  which  determine  the  branching  characteristics 
of  the  RCCC  mechanism.   Since  the  branching  condition 
depends  only  on  the  directions  of  the  joint  axes,  it  is 
possible  to  find  a  suitable  set  of  joint  axes  before 
dimensional  synthesis  of  the  mechanism  is  complete 
(see  Chapter  Four,  section  4.6). 

5.5   The  Order  Condition 

The  order  problem  was  introduced  briefly  in  Chapter 
Three.   The  order  of  positions  of  a  moving  body  is 
dependent  on  both  the  sequence  in  which  the  positions 
are  traversed  and  the  sense,  or  direction,  of  traversal 
(88).   For  example,  positions  1,2,3  and  3,2,1  are  in 
the  same  sequence  but  have  opposite  sense.   Since  the 
sense  can  be  reversed  by  reversing  the  direction  of 
the  mechanism's  input  crank  rotation,  only  the  sequence 
of  positions  will  be  of  concern  here.   Notice  that  only 
one  sequence  exists  for  three  positions,  i.e.,  1,2,3; 
2,3,1;  3,1,2;  etc.,  are  all  in  the  same  sequence.   Thus 
the  order  problem  can  only  exist  when  more  than  three 
positions  are  specified. 

Most  often  in  mechanism  design  problems  the 
specified  positions  naturally  occur  in  a  smooth, 
well-behaved  sequence.   In  such  cases  it  is  usually 
quite  easy  to  satisfy  the  sequence  requirement,  since 
the  output  of  most  mechanisms  also  tends  to  be  smooth 
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and  well-behaved.   However,  it  is  still  necessary  to 

test  if  the  sequence  of  output  positions  is  correct. 

This  can  be  done  by  finding  the  rotation  of  the  input 

link  at  each  of  the  prescribed  positions.   If  these 

rotations  are  in  a  successively  increasing  or  decreasing 

sequence,  the  order  condition  can  be  satisfied. 

Referring  to  figure  5.1,  6   is  the  input  angle 

measured  from  a41  to  a±2    in  a  right-hand  sense  about  s  . 

If  the  vector  (a,_) .  represents  the  jth  position  of  a,0, 
~  J-^  J  ~  12 

then  the  jth  position  of  6,  9.,  is  given  by 


and 


i  .  =  arccos{ (§12)  .  ■  a41> 


if  (a12^  i  *  ^41  x  s-,  >  0 


I.  =  2tt  -  arccos{(a12)  .  •  a41)  } 


if  (a]_2)  -;  *  a4i  x  s,  <  0 


(5.15) 


(5.16) 


Let  91'62'93  and  94  be  the  values  of  9  at  four  prescribed 
positions.   The  sequence  of  these  four  positions  will 
be  correct  if,  after  arranging  the  four  angles  in  a  series 
from  the  smallest  to  the  largest,  any  of  the  following 
eight  combinations  result: 
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Dl'a2'°3'°4  y4'S3'e2'9l 

e2,e3,e4,e1        e3,e2,e1,e4 

(5.17) 

e3,e4,e1,e2        62/6i'e4'e3 

G4'ei'92'e3  91,64,83,e2 

Notice  that  the  left-hand  set  above  is  in  the  correct  order, 
i.e. ,  the  positions  are  traversed  in  the  correct  sequence 
and  sense.   The  right-hand  set  will  be  in  the  correct 
order  if  the  direction  of  input  crank  rotation  is  reversed. 

The  procedure  just  described  is  a  nonparametric 
method  for  determining  if  the  prescribed  positions  are 
satisfied  in  the  correct  order. 

5- 6   Fixed-Pivot  Location  and  Link-Length  Ratio  Conditions 

The  fixed-pivot  locations  are  the  points  in  the 
global  reference  system  where  the  mechanism  is  physically 
connected  to  the  fixed  machine  frame.   For  both  the  RC 
and  the  CC  dyad,  such  fixed  points  are  located  by  the 
vector,  aQ,  as  shown  in  figure  4.9.   Restrictions  on 
the  value  of  aQ  will  be  strongly  dependent  on  the  specific 
problem  at  hand.   However,  these  constraints  will  gener- 
ally be  quite  easy  to  apply.   For  example,  the  constraint 

a0x  <  °  <5-18) 
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could  be  used  to  restrict  the  location  of  the  fixed 
pivot  to  the  left  of  the  plane  defined  by  the  Y  and  Z 


axes 


Link-length  ratio  conditions  will  also  be  quite 
easy  to  apply  in  practice.   Referring  to  figure  5.1, 
a12'  a23'  a34  and  a41  are  the  kinematic  link  lengths 
of  the  RCCC  mechanism.   One  possible  requirement  might  be 
that  the  ratio  of  the  longest  to  the  shortest  link 
length  be  less  than  ten,  or 


longest    . n         n 
10  <  0  (5.19) 

shortest 


It  should  be  noted  that  the  RCCC  mechanism  may  be  viable 
even  when  ashortest  is  equal  to  zero.   It  may  therefore  be 
desireable  to  restrict  only  the  maximum  allowable  link 
length.   Both  of  the  above-mentioned  constraints  will 
generally  be  nonparametric. 

5.  7   The  Objective  Function 

As  discussed  in  Chapter  Two,  the  objective  function 
serves  as  a  measure  of  the  relative  merits  of  a  particular 
design.   Since  each  of  the  requirements  discussed  in  this 
chapter  must  be  included  in  determining  the  relative 
merits  of  a  particular  design,  the  objective  function 
(O.F.)  will  be  of  the  following  form: 
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O.F.  =  M.E.T.  +  C.R.T.  +  B.A.T.  +  O.P.T. 

+  F.P.L.T.  +  L.L.R.T.  (5.20) 

where 

M.E.T.  =  Motion  Error  Term 
C.R.T.  =  Crank-Rotatability  Term 
B.A.T.  =  Branch-Avoidance  Term 
O.P.T.  =  Order-of-Positions  Term 
F.P.L.T.  =  Fixed-Pivot  Location  Term 
L.L.R.T.  =  Link-Length  Ratio  Term 

Motion  Error  Term.   The  first  term  in  the  above 
expression,  M.E.T.,  will  result  from  application  of  equation 
5.7  to  the  CC  dyad  and  the  RC  dyad.   The  total  error  will  be 
defined  as  the  sum  of  the  absolute  values  of  the  resulting 
errors  and  will  provide  a  direct  measure  of  the  success  in 
approximating  the  fourth  prescribed  position.   The  remaining 
terms  are  all  constraint  conditions  which  will  be  incorpo- 
rated using  penalty  function  methods.   In  many  cases,  the 
most  difficult  problem  is  finding  a  suitable  measure  of  the 
amount  by  which  a  constraint  is  violated.   A  brief  discussion 
of  each  constraint  condition  follows: 

Crank-Rotatability  Term.   Part  of  this  term  follows 
directly  from  equation  5.8  and  the  penalty  function  methods 
of  section  2.4.2.   However,  an  additional  constraint  is 
needed  to  ensure  that  the  input  link,  a   ,  is  the  shortest 
link.   This  can  be  accomplished  by  assessing  the  following 
penalty  when  a,2  is  not  the  shortest  link: 
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Aa   =    largest   of 

{(ai2    -   a23)'     (ai2    "   a34>'     (ai2    "   a4l)}  (5'21> 

where  the  angular  brackets  Indicate  the  singularity 
function  as  defined  by  equation  2.30.   It  can  be  seen 
that  Aa  will  equal  zero  when  a,-  is  the  shortest  link. 
Otherwise  Aa  will  be  a  positive  number  which  indicates 
the  amount  by  which  a,,  differs  from  the  shortest  link. 
With  this,  the  crank  rotatability  term  becomes 

C-R-T-  =  r{  <a12  +  90°  -  i,(o12  +  a23  +  a34  +  a^  +  Aa}     (5.22) 

where  r  is  the  penalty  parameter  defined  in  section  2.4.2. 
Equation  5.22  will  equal  zero  when  the  input  link  of  the 
RCCC  mechanism  is  a  crank.   Any  result  greater  than  zero 
serves  as  an  indication  of  the  amount  by  which  the 
mechanism  fails  to  meet  this  requirement. 

Branch-Avoidance  Term.   The  branch-avoidance  condi- 
tions are  given  by  equations  5.13  and  5.14.   For 
convenience,  the  function  F,  is  defined  to  be  equal  to 
the  left-hand  side  of  these  equations,  i.e. 

Fb  =  (!3  "  !2}  x  (!4  "  !3>  '  §3  (5-23) 

so  that  F   must  be  either  greater  than  zero  or  less  than 
b 

zero  in  all  positions  if  branching  problems  are  to  be 
avoided.   When  branching  does  occur,  it  will  be  necessary 
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to  have  some  measure  of  the  magnitude  of  the  problem. 
One  such  measure  is  given  by 


B.A.T.  =  r{smallest  of  (EFbg , | ZF^ \ ) }  (5.24) 


where  ZF   is  the  sum  of  the  Ffa '  s  greater  than  zero  at 
the  prescribed  positions  and  |£Fb£|  is  the  absolute  value 
of  the  sum  of  all  F.  's  less  than  zero  at  the  prescribed 
positions.   When  all  Fb's  are  of  the  same  sign,  B.A.T. 
will  equal  zero.   Notice  that  this  term  will  only  be  an 
approximation  for  positions  which  are  approximately 
satisfied. 

Order-of-Positions  Term.   Equation  5.17  gives  all 
possible  combinations  of  the  four  values  of  6  which  are 
in  the  correct  sequence.   Again  the  problem  arises  of 
what  penalty  to  assess  when  this  constraint  is  violated. 
Regardless  of  the  arrangement  of  the  four  6  values,  a 
correct  sequence  can  always  be  obtained  by  interchanging 
two  values  of  8.   However,  two  possibilities  always 
exist  for  this  interchange.   Therefore,  it  is  necessary 
to  scan  the  possible  choices  to  find  the  most  closely 
spaced  pair  of  angles  which  can  be  interchanged  to  give 
the  correct  sequence.   These  angles  will  be  called  9 
and  9,  .   Now  the  difference  between  these  two  angles 
can  serve  as  a  measure  of  the  error  in  satisfying  the 
order  condition,  i.e. 
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O.P.T.  =  r 


'a  "  9b'  <5-25) 


When  the  sequence  is  correct,  9   and  9,  will  be  set 

a      d 

equal  to  zero. 

Fixed-Pivot  Location  Term.   Application  of  this 
constraint  condition  will  usually  be  simple  and  direct. 
The  components' of  the  fixed  pivot  location  vector  a 
will  most  often  be  restricted  to  lie  within  a  specified 
range,  for  example 

xa  K   a0x  K   xb  (5-26) 

This  contributes  the  following  term  to  the  objective 
function 

F.P.L.T.  -  r{(xa  -  aQx>  +  (aQx  -  xb  >  }      (5.27) 

where,  again,  r  is  the  penalty  parameter. 

Link-Length  Ratio  Term.   Application  of  this 
constraint  will  also  usually  be  simple  and  direct.   For 
example,  equation  5.19  would  result  in  the  following 
penalty  term: 


L.L.R.T.  =  r  <  (a,      Va  ,        )  -  10  )       (5  2! 
vv  longest7  shortest'       '  ,-.^< 


It  can  be  seen  from  the  foregoing  discussion  that 
formulating  the  objective  function  is  far  from  an  exact 
science.   This  work  represents  the  first  known  attempt 
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to  apply  many  of  the  nonparametric  constraint  conditions 
discussed  in  this  chapter  to  the  optimization  of  the 
RCCC  mechanism.   As  a  result,  no  basis  exists  for 
comparing  the  formulations  presented  here.   Future 
research  will  undoubtedly  lead  to  improved  and  more 
sophisticated  methods  for  assessing  penalties  when 
constraints  are  violated. 

5. 8   Numerical  Example 

The  following  positions  and  orientations  of  the 
spatial  moving  body  are  specified: 


o,  =  {Oi  +  Oj  +  Ok} 
o2  =  {4i  +  2j  +  Ok}' 
o,  =  {8i  +  2j  +  3k}r 


o.  =  {10i  +  3j  +  5k} 


.966 

-.258 

.008 

.258 

.962 

-.087 

.015 

.086 

.996 

.868 

-.494 

.045 

.489 

.834 

-.255 

.089 

.243 

.966 

.824 

-.564 

.059 

.548 

.765 

.337 

.145 

.310 

.940 
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Recall  from  section  4.4  that  [R  ]  is  assumed  to  be  a  3x3 
identity  matrix. 

Problem  Statement.   Design  an  RCCC  mechanism  which 
exactly  satisfies  the  first  three  of  the  above  positions 
and  orientations.   The  fourth  position  and  orientation 
are  to  be  satisfied  in  an  approximate  sense  using  the 
position-stacking  method.   The  solution  mechanism  is 
required  to  be  of  the  crank-rocker  type,  free  from 
branching  problems.   The  fixed  pivots  should  be  located 
within  one-hundred  units  of  the  global  reference  origin 
and  the  ratio  of  the  longest  to  the  shortest  link  length 
should  not  exceed  ten  units. 

Solution.   Solution  mechanisms  are  obtained  by 
minimizing  the  objective  function  of  equation  5.20  using 
the  Hooke  and  Jeeve's  optimization  method.   Appendix  1 
provides  a  complete  listing  of  all  the  APL-language 
computer  programs  used  in  solving  this  problem  along  with 
a  sample  run  of  the  program  package. 

Figure  5.4  shows  the  RCCC  mechanism  and  associated 
design  parameters.   Notice  that  the  fixed  and  moving 
coordinate  systems  are  coincident  in  the  initial  position. 
Tables  5.1  and  5.2  give  the  results  of  two  separate 
program  runs  which  converged  to  feasible  solutions. 
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/ 


~4 


figure  5.4   The  RCCC  mechanism  and  associated  design 
parameters 
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Table  5.1 
RCCC  Mechanism  Optimization  -  Example  1 
Assumed  quantities  (refer  to  figure  5.4) 

§x  =  {0|  +  lj  +  Ok}  si   =  si   =  3 

S   =  (.9811  +  .1733  -  .086k}  (displacements  along  s 


Z4' 


After  35  iterations  of  the  APL  program  HOOKE  and  a  total  of 
601  evaluations  of  the  objective  function,  all  constraints  were 
satisfied  and  the  solution  had  converged  to  the  results  given  below. 


i 

J 

k 

?0 

= 

-22 

.634 

13 

.635 

-9 

.229 

^ 

~S2 

= 

.281 

- 

.363 

.889 

r 
~a 

= 

-18 

.595 

14 

.848 

-10 

.012 

*0 

= 

-12 

.188 

29 

.798 

3 

.518 

h 

= 

- 

.307 

- 

.201 

.930 

r~h 

= 

-  7 

.233 

17 

.979 

2 

.594 

L 

I 

k 

S0 

= 

-23 

349 

13 

.601 

-9 

.864 

A 

-S2 

= 

281 

- 

363 

889 

r 
~a 

= 

-19 

142 

14 

865 

-10 

680 

*0 

= 

-5 

656 

16 

266 

2 

341 

53 

= 

- 

469 

-. 

241 

850 

?b 

= 

-6. 

859 

19. 

197 

2 

510 

mechanism 
satisfying 
positions 
1,2  and  3 


mechanism 
satisfying 
positions 
1,2  and  4 


For  both  mechanisms 

S   =  {-314i  -  .542j  +  -780k}     SI   =  -.089      SI   =  .781 
S   =  {.869'i  +  .336j  +  .363k/ 
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Table  5.2 

RCCC  Mechanism  Optimization  -  Example  2 

Assumed  quantities  (refer  to  figure  5.4) 

l±   =   {Oi  +  lj  +  Ok}  SI   =  SI   =  -3 

f4  =  tOi  +  l3  +  Ok}  (displacements  along  s  ) 

After  25  iterations  of  the  APL  function  HOOKE  and  a  total  of 
541  evaluations  of  the  objective  function,  all  constraints  were 
satisfied  and  the  solution  had  converged  to  the  results  given  below. 


i 

J 

k 

to 

= 

-15 

.271 

6 

.477 

-62 

.412 

§2 

= 

.725 

.457 

.515 

r 

~a 

= 

-12 

.704 

2 

.653 

-62 

.633 

bo 

= 

-10 

495 

10 

.246 

-20 

.034 

§3 

= 

137 

- 

.061 

.989 

5b 

= 

-12 

092 

14 

134 

-19 

574 

i 

] 

k 

-10 

.833 

7 

.950 

-42 

.64 

.720 

.447 

.529 

-9 

041 

5 

.229 

-42 

.771 

-13 

179 

22 

411 

-12 

.152 

130 

- 

063 

.989 

-10 

393 

15 

642 

-12 

947 

mechanism 
satisfying 
positions 
1,2  and  3 


~0 

!2 

mechanism 
satisfying 

bn   =   -13.179     22.411     -12.152  positions 

~°  1,2  and  4 

?3 

5b 

For  both  mechanisms 


s   =  {.564i  +  .335]  +  .755k}     SI   -  .375     SI 


-.019 


s4  =  {.6681  +  .352]  +  -.656k} 


CHAPTER  6 
OPTIMIZATION  OF  THE  RSSR-SC  AND  RSSR-SS  MECHANISMS 

6 . 1   Problem  Formulation 

Chapter  Five  presented  what  is  believed  to  be  the 
most  complete  discourse  to  date  on  designing  the  RCCC 
mechanism  for  spatial  rigid-body  guidance.   However, 
several  important  problems  remain  to  be  solved  in  the 
design  of  this  mechanism.   First,  nonparametric  condi- 
tions are  needed  to  determine  its  transmission 
characteristics.   These  conditions  must  include  consid- 
eration of  the  forces  which  cause  binding  of  the 
cylindric  joints  along  their  axes.   Second,  the  problem 
of  cylindric  joint  axes  becoming  collinear  during  the 
motion  cycle  must  be  addressed.   When  this  occurs,  the 
displacements  along  these  joint  axes  go  to  infinity,  and 
the  mechanism  becomes  useless.   Until  these  problems  are 
solved,  the  potential  uses  of  the  RCCC  mechanism  will 
remain  quite  limited. 

This  chapter  discusses  the  design  of  the  five-link 
RSSR-SC  and  RSSR-SS  mechanisms  for  spatial  rigid-body 
guidance.   These  mechanisms  are  shown  in  figures  6.1  and 
6.2.   As  a  result  of  the  problems  discussed  above,  the 
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figure    6.1      The    RSSR-SC   mechanisi 
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figure  6.2   The  RSSR-SS  mechanism 
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author  believes  these  mechanisms  are  currently  more  useful 
than  the  RCCC  mechanism.   Also,  since  the  five-link  mecha- 
nisms presented  in  this  chapter  have  three  grounded  joints, 
they  will  tend  to  be  more  stable  in  holding  a  certain  position 
than  the  RCCC  mechanism,  which  has  only  two  grounded  joints. 

The  objective  of  this  chapter  is  to  develop  procedures 
for  designing  RSSR-SC  and  RSSR-SS  mechanisms  for  three 
positions  of  rigid-body  guidance.   Constraints  dealing  with 
branch  avoidance,  Grashof  type,  transmission  characteristics, 
link-length  ratios  and  fixed  pivot  locations  will  be 
developed  and  applied.   Also  discussed  are  the  use  of  the 
position  stacking  method  and  the  order  problem  when  more 
than  three  positions  are  specified. 

6. 2   Method  of  Design 

The  RSSR-SC  and  RSSR-SS  mechanisms  may  be  synthesized 
for  three  positions  of  rigid-body  guidance  using  the 
methods  of  Chapter  Four.   This  requires  finding  two  RS 
dyads  and  either  a  CS  or  an  SS  dyad  whose  floating  links 
satisfy  the  specified  motion.   The  floating  links  of 
these  dyads  are  then  connected  to  form  a  single-degree-of- 
freedom  mechanism. 


The  RSSR-SS  mechanism  possesses  an  additional  idle  degree 
of  freedom  which  is  the  free  rotation  of  the  grounded  link 
of  the  SS  dyad  about  its  own  axis. 
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Ideally,  the  remaining  requirements,  such  as  branch- 
avoidance,  Grashof  type,  etc.,  could  be  formulated  non- 
parametrically  for  each  mechanism,  as  was  done  for  the 
RCCC  mechanism  in  the  previous  chapter.   Unfortunately, 
the  relative  complexity  of  the  RSSR-SC  and  the  RSSR-SS 
mechanisms  renders  this  approach  impractical,  at  least 
at  the  present  time.   An  alternate  approach,  based  on 
the  design  formulation  of  the  RSSR  mechanism  is  now 
proposed. 

The  RS  dyad  may  be  thought  of  as  a  special  case 
of  both  the  CS  and  SS  dyads.   The  RS  dyad  is  obtained 
from  the  CS  dyad  by  removing  the  axial  translation, 
and  obtained  from  the  SS  dyad  by  removing  two  rotational 
freedoms  from  the  grounded  S  joint.   Three  RS  dyads 
may  be  synthesized  for  three  positions  of  rigid-body 
guidance  and  assembled  into  an  RSSR-SR  structure,  as 
shown  in  figure  6.3.   This  structure  can  be  assembled 
at  the  prescribed  positions,  but  cannot  move  between 
them.   However,  replacing  any  of  the  grounded  revolute 
joints  by  either  a  cylindric  joint  or  a  spheric  joint 
results  in  viable,  single-degree-of-f reedom  mechanisms 
capable  of  generating  the  prescribed  motion.   By 
designing  the  RSSR-SC  and  RSSR-SS  mechanisms  in  this 
way,  they  may  be  considered  to  be  a  pair  of  RSSR 
mechanisms  with  one  dyad  in  common,  at  least  at  the 
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figure  6.3   The  RSSR-SR  structure 
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precision  positions.   This  will  greatly  facilitate 
finding  nonparametric  conditions  for  the  various  design 
requirements,  although  some  of  these  conditions  will  be 
approximations.   For  the  remainder -of  this  chapter,  the 
design  of  the  RSSR  and  the  RSSR-SR  mechanism  will  be 
used  as  an  intermediate  step  in  the  design  of  the  RSSR-SC 
and  the  RSSR-SS  mechanisms.   However,  the  constraint 
conditions  developed  for  the  RSSR  mechanism  are  extremely 
important  in  their  own  right,  since  this  is  one  of  the 
most  often  used  spatial  mechanisms  for  function  generation. 

6. 3   The  Branch-Avoidance  Condition 

The  mechanism  branching  problem  has  been  discussed 
in  sections  3.2  and  5.3.   It  is  evident  that  kinematic 
position  analysis  is  often  too  inefficient  to  incorporate 
in  an  optimization  loop. 

Fortunately,  the  recent  works  of  Gupta  and  Tinubu 
(92)  and  Sandor  and  Zhuang  (93,94)  have  eliminated  the 
need  for  this  time-consuming  analysis  by  devising  non- 
parametric  tests  to  determine  if  the  branching  problem 
is  present  in  the  RSSR  mechanism.   Since  the  method  of 
Sandor  and  Zhuang  (93,94)  applies  directly  to  motion- 
generating  mechanisms,  and  since  it  can  be  applied  to  a 
variety  of  other  mechanisms  as  well,  it  will  be  demon- 
strated here.   It  is  important  to  realize  that  the 
present  analysis  applies  exactly  to  the  RSSR-SC  and 
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RSSR-SS  mechanisms.  This  is  true  because  the  branching 
conditions  are  applied  only  at  the  precision  positions, 
where  the  configuration  of  the  RSSR-SR  mechanism  is  the 
same  as  these  other  mechanisms. 

Figure  6.4  shows  a  schematic  representation  of  the 
RSSR  mechanism.   The  two  branches  can  be  determined 
by  holding  input  link  AQA  fixed  and  finding  the  inter- 
sections of  the  loci  of  spheric  joint  B  about  spheric 
joint  A  and  of  spheric  joint  B  about  revolute  joint  B- . 
These  loci  are  a  sphere  and  a  circle,  respectively,  as 
shown  in  figure  6.5. 

It  has  been  shown  that  when  the  two  branches 
coincide,  links  AB ,  BQB  and  axis  IL  are  coplanar  (93); 
this  is  a  so-called  special  configuration  of  the 
mechanism.   When  the  branches  are  separated,  the  two 
intersections  will  be  located  on  opposite  sides  of  this 
plane,  and  cannot  change  sides  without  passing  through 
a  special  configuration  of  the  mechanism.   Using  the 
subscript  j  to  denote  the  prescribed  positions  of  the 
mechanism,  the  above  logic  leads  to  the  conclusion  that 
either 

A.B.  •  s~  x  bTb .  >  0      for  all  j  (6.1) 

j  j    ~  *■  u  j 


A.B.  •  L  x  BnB.  <  0      for  all  j  (6.2 

]  3    ~2     0  j 

if  the  positions  j  are  to  lie  on  one   branch  only. 
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figure  6.4   The  RSSR  mechanism 
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figure  6.5   The  two  branches  of  the  RSSR  mechanism 
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It  should  be  noted  that  the  vector  A.B .  is  simply 
the  vector  difference  in  the  locations  of  the  two 
spheric  joints  for  a  pair  of  independently  synthesized 
RS  dyads.   For  example,  if  rA  and  rB  are  used  to  denote 
the  location  of  the  spheric  joint  relative  to  the 
moving  coordinate  system  for  dyads  A  and  B,  then 


AJBJ  =  Sj  +  [Rj]~rB  "   (?j  +  [rj]5a} 

=  [Rj^B  "  tRj^A 

=  [Rj] (rB  -  rA)  (6.3) 

The  vector  representing  the  output  link,  B~B  .  ,  is 
given  by 


B0Bj  =  ?j  +  ([Rj^B}  "  !?0  (6'4> 

where  bQ  locates  the  grounded  revolute  joint  center  of 
dyad  B.   It  can  be  seen  that  constraint  conditions  6.1 
and  6.2  are  simple,  nonparametric  expressions  that  can 
be  used  to  identify  RSSR  mechanisms  without  branching 
problems . 

As  has  already  been  discussed,  the  RSSR-SR  mechanism 
can  be  considered  to  be  two  RSSR  mechanisms  with  one  dyad 
in  common.   The  forgoing  analysis  may  therefore  be  applied 
to  the  two  RSSR  mechanisms  to  determine  the  branching 
characteristics  of  the  overconstrained  RSSR-SR  mechanism, 
and  hence,  the  RSSR-SC  and  RSSR-SS  mechanisms. 
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6.  4   The  Grashof  Condition 

Sections  3.2  and  5.3  presented  brief  discussions 
on  the  nature  and  importance  of  the  Grashof  condition. 
The  goal  of  this  section  is  to  determine  input-link 
rotatability  conditions  for  the  RSSR-SC  and  the  RSSR-SS 
mechanisms.   Unfortunately,  these  conditions  will  not 
be  exact,  since  they  will  be  based  on  the  input-link 
rotatability  conditions  for  an  independent  pair  of  RSSR 
mechanisms.   When  the  coupler  links  of  the  RSSR  mechanisms 
are  joined,  and  one  of  the  revolute  joints  is  replaced  by 
a  spheric  or  a  cylindric  joint,  the  motion  of  that  joint 
will  no  longer  be  pure  rotation.   Because  of  this,  the 
Grashof  condition  for  the  RSSC  or  RSSS  loop  of  the  mechanism 
will  be  an  approximation.   While  this  analysis  will  not 
guarantee  the  correct  mechanism  type,  it  is  likely  to 
yield  an  acceptable  solution  in  far  fewer  iterations  of 
position  analysis  than  would  otherwise  be  required. 
Additionally,  it  should  be  noted  that  either  revolute 
joint  which  is  not  the  input  may  be  replaced  by  a 
cylindric  or  spheric  joint.   It  is  quite  possible  that 
one  case  will  yield  an  acceptable  solution  when  the 
other  case  does  not.   Therefore,  both  possibilities 
should  be  explored  when  searching  for  an  acceptable  design. 

Several  authors  have  devised  nonparametric  methods 
to  identify  crank  rotatability  of  the  RSSR  mechanism 
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(86,118,119,120).   Adopting  the  notation  of  Nolle  (119), 
the  RSSR  mechanism  is  shown  in  figure  6.6.   The  input  link 
axis  is  coincident  with  the  Z  axis,  and  the  origin  of  the 
coordinate  system  is  located  at  the  intersection  of  this 
axis  and  the  common  perpendicular  to  the  two  revolute 
joint  axes.   The  input  angle,  9,  and  the  output  angle,  <j> , 
are  measured  from  lines  parallel  to  the  X  axis  in  a  right- 
hand  sense  about  the  respective  joint  axes.   In  terms  of 
the  notation  of  Chapter  Four,  the  link  lengths  a,b  and  c 
are  simply  the  magnitudes  of  the  previously  determined 
vectors  representing  these  links.   The  length  of  the 
common  normal,  e,  is  given  by 

e  =  (bQ  -  aQ)  •  (  {s1   x  §2)  4-  \s^   x  s2|)  (6.5) 

where  a  ,  b_,  s..  and  s„  are  the  revolute  joint  locations 
and  the  joint  axis  unit  vectors,  as  shown  in  figure  6.4. 
The  scalar  quantity,  f,  can  be  determined  from 

e(s1  x  s2)  +  gs2  +  (aQ  -  bQ)  =  fs,  (6.6) 

Cross-multiplying  both  sides  by  s_  eliminates  the  term 

containing  g  and  gives 

o 

(6.7) 

Finally,  taking  the  scalar  product  of  both  sides  of  equation 
6.7  with  the  vector  (s,  x  s~)  and  solving  for  f  gives 


(an  -  bn)  '  s9  x  (s   x  s9) 

f  =  — ~ — — -£—  (6.8) 

(§1  x  §2)  •  (s.^  x  §2) 
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figure  6.6   The  RSSR  mechanism  in  the  coordinate  system 
of  Nolle  (119) 
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Similarly,  g  can  be  shown  to  be 


(§0  "  b_)  x  s   •  (s   x  s„) 

g  =  : : - —  (6.9) 

(s1  x  s2)  •  (s.  x  |2) 


The  angle  between  the  revolute  joint  axes,  ty ,    is  given 
by 

i|)  =  arccos  (s^  •  s"2)  (6.10) 

It  should  be  noted  that  negative  values  of  ij;  need  not 
be  considered,  because  the  Grashof  type  is  a  symmetric 
property  about  the  value  ijj  =  0 . 

The  relationship  between  the  input  angle,  0,  and 
the  output  angle,  §,    has  been  shown  to  be  (118) 

1    ±     [1    +    f!{Q)    -    F2(0)  }h 

tan    <t>/2    = (6.11) 

F1(9)    -    F2(0) 

where 


k      +  k„cos0    +  k^sin© 

F1(9)    »   — (6.12) 

k,    +   k_sin0 


k.    +   k-cosO 

F2(9)    =    — (6.13) 

kfi    +   k_sin0 


and 


154 


k,  =  fg  cos^  -  ^(f2  +  a2  +  e2  +  g2  +  c2  -  fc,2) 


(6.14; 


k2    =    ae  (6.15) 

k3  =  ag  sintjj  (6.  16) 

k4  =  -ce  (6.17) 

k5  =  ca  (6.18) 

kg  =  cf  simjj  (6.19) 

k   =  -ca  cos\\i  (6.20) 

The  two  roots  of  equation  6.11  correspond  to  the 

two  possible  values  of  the  output,  <j> ,  for  any  given  value 
of  the  input,  6.   However,  when  the  input  link  reaches  a 

limit  position,  the  two  values  of  <j>  become  the  same  (118). 

For  this  to  happen,  the  quantity  within  the  radical  of 
equation  6.11  must  equal  zero,  i.e., 

1  +  PjO)  +  fJ(8)  =  0  (6.21) 

2  2 

Since  F1(6)  and  F_(0)  are  quadratic  in  cosG  and  sine , 

the  following  substitutions  are  made: 

cosG  =  x  (6.22) 

sine  =  y  (6.23) 

y2  =  1  -  x2  (6.24) 
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With  these,  equation  6.21  may  be  rewritten  as 

4      3      2 
x   +  Bx   +  Cx   +  Dx  +  E  =  0  (6.25) 


where 


B  =  (2K2K5  +  2K3K4)/(K^  +  Kg)  (6.26) 

C  =  (K2  +  K3  -  K2  +    2K1K5)/(K^  +  K2)           (6.27) 

D  =  (2KXK5  -  2K3K4)/(kJ  +  K2)  (6.28) 

E  =  (K2  -    K2)/(K2  +  Kg)  (6.29) 


and 


Kl  =  k4  +  k6  +  k7  "  kl  "  k3  (6'3°) 

K2  =  2k4k5  "  2klk2  (6.31) 

K3  =  2k6k7  "  2klk3  (6.32) 

K4  =  "2k2k3  (6.33) 

2     2     2     2 

K5  =  k3  +  k5  "  k2  "  k7  (6.34) 

Solving  equation  6.25  will  generally  yield  four 
values  of  9  corresponding  to  the  four  possible  limit 
positions.   However,  if  the  input  link  is  a  crank,  no 
limit  positions  should  exist,  and,  consequently, 
equation  6.25  should  have  no  real  roots.   It  is  possible 
to  find  the  roots  of  a  quartic  equation  in  closed  form. 
However,  to  determine  whether  or  not  the  input  link  is 
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a  crank,  it  is  only  necessary  to  determine  whether  or 
not  real  roots  of  the  quartic  equation  exist.   Necessary 
and  sufficient  conditions  for  the  absence  of  real  roots 
in  a  quartic  equation  are  given  by  Dickson  (121) .   To 
apply  these  conditions,  it  is  first  necessary  to  replace 
x  by  z  -  B/4  in  equation  6.25  to  obtain  the  reduced 
quartic  equation 

4     ,         2 
z      +qz      +rz+s=0  (6.35) 

where 


q  =  -3B2/8  +  C  (6.36) 

r  =  B3/8  -  CB/2  +  D  (6.37) 

s  =  -3B4/256  +  CB2/16  -  DB/4  +  E  (6.38) 


From  the  reduced  quartic,  the  discriminant,  A,  and  a 
quantity  L  are  defined 

A  =  -4P3  -  27Q2  (6.39) 

3      2 

L  =  8qs  -  2q   -  9r  (6.40) 

where 

P  a  -4s  -  q2/3  (6.41) 

Q  -  8qs/3  -  r2  -  2q3/27  (6.42) 
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The  absence  of  real  roots  is  assured  when  the  following 
conditions  are  satisfied  (121,  p.  100): 

A  >  0  (6.43) 

q  >  0       or      L  <  0  (6.44) 

Although  somewhat  complex  to  derive,  the  above  expressions 
are  simple,  nonparametric  constraint  conditions  which 
determine  the  input  link  rotatability  of  the  RSSR  mech- 
anism.  It  should  be  noted  that  these  same  conditions  may 
also  be  applied  to  the  output  link  to  determine  its 
rotatability. 

Again  it  should  be  emphasized  that  the  RSSR-SC  and 
RSSR-SS  mechanisms  may  not  have  fully  rotating  input 
links  even  though  the  RSSR  mechanisms  may  satisfy  these 
conditions.   On  the  other  hand,  these  mechanisms  may 
have  fully  rotating  input  links  even  though  one  of  the  RSSR 
mechanisms  does  not.   In  this  case,  the  non-input  revolute 
joint  on  the  RSSR  mechanism  that  does  not  have  a  fully  rotat- 
ing input  link  must  be  the  joint  replaced  by  a  spheric  or 
cylindric  joint.   Finally,  if  neither  RSSR  mechanism  has 
a  fully  rotating  input  link,  the  input  links  of  the 
RSSR-SC  and  RSSR-SS  mechanisms  will  also  not  rotate  fully. 

6 . 5   The  Transmission  Characteristic  Condition 

As  was  the  case  for  the  Grashof  condition  in  the 
previous  section,  the  transmission  characteristic 
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conditions  developed  in  this  section  will  not  be  exact 
for  the  RSSR-SC  and  RSSR-SS  mechanisms.   This  is  true 
because  the  transmission  characteristics  are  of  concern 
over  a  full  cycle  of  the  mechanism  rather  than  at  the 
precision  positions  alone.   The  approach  taken  will  be 
to  optimize  the  transmission  characteristics  of  two 
RSSR  mechanisms  with  one  dyad  in  common.   These  mechanisms 
are  than  combined  to  give  either  the  RSSR-SC  or  RSSR-SS 
mechanism.   Once  again,  this  approach  will  not  guarantee 
an  acceptable  design.   It  does,  however,  provide  a 
nonparametric  basis  for  decision-making  by  which  many 
unacceptable  solutions  may  be  rejected. 

Following  the  work  of  Sutherland  and  Roth  (122)  and 
Soylemez  and  Freudenstein  (78) ,  the  transmission  ratio  in 
the  RSSR  mechanism  is  defined  as  the  fraction  of  force  in 
the  coupler  which  acts  to  produce  output  rotation.   Using 
the  notation  of  figure  6.6  and  of  the  previous  section,  the 
transmission  ratio,  TR,  may  be  expressed  as 

TR  =  [1  -  cos  m  -  sin2a]  2  (6.45) 

where  y  is  the  angle  between  the  coupler  link,  b,  and 
the  output  link,  c,  and  a  is  the  angle  between  the 
coupler  link,  b,  and  the  plane  normal  to  the  unit  vector 
along  the  output  link  axis,  s~. 

Using  the  geometric  relations  within  the  mechanism 
to  replace  p  and  a  in  equation  6.45  leads  to  the  result  (122) 
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12  9 

TR  =  [4c    (fsirnji   -   acos^sine) 

2bc 

+   4c2  (e   +   acos9)  2 


-  (R  +  2aecos0  -  2agsinii;sin6)  y-  (6.46) 


where 


2     2     2     2     2     2 

R=a   -b   +  c   +e   +  f   +  g   -  2fgcos^       (6.47) 

The  maximum  and  minimum  values  of  TR  may  be  determined 
by  taking  the  derivative  of  the  above  expression  with 
respect  to  6,  and  equating  the  result  with  zero,  to 
give 

Acos6  +  BsinG  +  Csin20  +  Dcos28  =  0  (6.48) 

where  A,  B,  C  and  D  are  now  defined  by 

2 

A  =  RgsiniJ;  -  2c  f sinijjcosijj  (6.49) 

B  =  Re  -  2c2e  (6.50) 

2     2       2     2       2    2 
C  =  c  acos  i|)  +  ae   -  c  a  -  ag  sin  ^  (6.51) 

D  =  2  (aegsintjj)  (6.52) 

At  this  point,  Soylemez  and  Freudenstein  (78) 
consider  only  special  cases  of  the  RSSR  mechanism; 
these  greatly  simplify  finding  the  roots  of  equation 
6.48.   To  extend  their  procedure  to  the  general  case, 
it  is  necessary  to  find  the  roots  of  equation  6.48. 
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given  any  mechanism  dimensions.   This  is  done  with  the 
help  of  the  following  trigonometric  identities: 


(6.53) 


sin6 

2t 

1 

+  t2 

COS0 

1 

-  t2 
+  t2 

1 

(6.54 


tane  =  — ^—j  (6.55 

1  -  t 


sin29  =  ^M^_  .    4t(l  -  t2)  (6  56. 

1  +  tan  6      (1  -  t  )z+  4t 

cos29  -  1  '  ^\*  -   (1  -  *?>'  -    ^        ,6.57: 
1  +  tan  6      (1  -  t  )   +  4t 


where  t  =  tan  0/2.   Substituting  these  into  equation 
6.47  yields,  after  some  manipulation, 

(D  -  A)t4  +  (2B  -  4C)t3  +  (-6D)t2 

+  (2B  +  4C)t  +  (A  +  D)  =0  (6.58) 

The  real  roots  of  the  above  quartic  will  correspond 
to  the  extreme  values  of  the  transmission  ratio.   The 
roots  of  equation  6.58  may  be  determined  using  any  of 
the  standard  quartic  root-finding  procedures  (121) .   For 
the  purpose  of  working  the  example  at  the  end  of  this 
chapter,  optimum  force  transmission  will  be  defined  as 
the  case  when  the  minimum  value  of  TR  has  been  maximized. 
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Since  two  RSSR  mechanisms  are  being  considered  simul- 
taneously, the  minimum  value  of  TR  must  be  selected 
between  the  two  mechanisms. 

6. 6   Fixed-Pivot  Location  and  Link-Length  Ratio  Conditions 

Once  again  the  constraint  equations  for  these 
conditions  will  usually  be  self-evident  and  easy  to  apply. 
The  fixed-pivot  locations'  for  the  RSSR-SC  and  RSSR-SS 
mechanisms  will  be  the  same  as  for  the  RSSR-RS  mechanism 
which  in  turn  will  be  the  same  as  for  the  pair  of  RSSR 
mechanisms  with  one  dyad  in  common.   A  typical  constraint 
condition  may  require  the  fixed  pivots  to  be  within  a 
given  radius  of  the  origin  of  the  global  reference  system. 
This  could  be  expressed  as 

|aQ|  <  RQ  (6.59) 

|b0l  -  R0  (6.60) 

where  R_  is  the  given  radius. 

For  the  purposes  of  this  section,  the  link  lengths 
will  be  defined  as  the  distances  between  consecutive 
joints  in  their  starting  position.   Thus,  for  example, 
the  ternary  coupler  link  will  result  in  three  link 
lengths,  and  the  entire  mechanism  will  have  a  total  of 
nine  link  lengths.   As  in  the  previous  chapter,  the 
ratio  of  the  longest  link  length  to  the  shortest  link 
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length  should  not  exceed  a  certain  value,  in  this  case 
ten 


length  of  longest  link 

— 10  <_   0  (6.61) 

length  of  shortest  link 


6. 7   Satisfying  Additional  Motion  Requirements 

Using  the  position  stacking  method  presented  in 
Chapter  Five,  it  is  possible  to  design  the  RSSR-SC  and 
RSSR-SS  mechanisms  to  satisfy  up  to  four  precision 
positions  plus  any  additional  number  of  approximate 
positions.   The  procedure  is  analogous  to  the  procedure 
for  the  RC  and  CC  dyads.   The  specified  precision 
positions  are  divided  into  groups  of  three.   The  optimi- 
zation process  then  seeks  to  find  the  most  similar  dyads 
satisfying  each  three-position  motion  group.   The  RS 
dyad  constraint  equations  for  position  stacking  must 
ensure  similarity  of  the  location,  a_ ,  and  the  orientation, 
s,,  of  the  revolute  joint.   Thus,  for  two  groups  of  three 
precision  positions,  denoted  by  the  superscripts  1  and  2, 
the  function  to  be  minimized  is 

F(r)  -  W]_|aJ  "  ajj|  +  W2  |  s^  -  s^  |  (6.62) 

where  F(r)  shows  a   and  s,  to  be  functions  of  the  location 
of  the  spheric  joint,  r,  in  the  moving  frame  of  reference 
and  where  W,  and  W_  are  weighting  factors. 
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Of  course,  this  procedure  must  be  used  in  designing 
all  three  RS  dyads  and  the  same  solution  must  be  chosen 
in  each  case  (i.e.,  either  aQ ,  sj  or  a2      sh  .       The 
procedure  just  described  will  add  some  complexity  to  an 
already  complex  problem.   Nevertheless,  it  is  at  least 
conceptually  quite  straightforward  to  apply.   It  also 
offers  a  solution  to  what  seems  to  be  a  practical 
problem,  namely,  partially-exact  partially-approximate 
synthesis.   As  shown  in  the  next  section,  however,  any 
attempt  to  satisfy  four  or  more  positions  using  the  RS 
dyad  must  be  done  with  careful  consideration  of  the 
order  problem. 

6.8   The  Order  Condition 

Section  5.5  discussed  the  general  order  problem  in 
some  detail.   It  was  demonstrated  that  if  the  sequence 
of  positions  is  correct,  the  order  requirement  can  always 
be  satisfied  by  choosing  the  correct  direction  of  input 
link  rotation.   Furthermore,  it  was  pointed  out  that  the 
sequence  of  positions  will  usually  be  correct  for  a  well- 
formulated  problem. 

This  is  extremely  fortunate  because,  based  on  the 
findings  of  Sun  and  Waldron  (123) ,  it  may  be  difficult 
or  impossible  to  find  RSSR  mechanisms  which  satisfy  a 
given  set  of  rigid-body  positions  in  more  than  one 
sequence.   This  conclusion  was  reached  by  noting  that 
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points  of  the  guided  body  which  lie  on  a  circle  for  four 
finitely  separated  spatial  positions  are  distributed  on 
a  sixth-order  spatial  curve  (110,111).   Changes  of 
sequence  can  only  occur  at  points  along  this  curve  where 
the  spheric  joint  has  the  same  location  for  the  two 
positions  of  the  moving  body  experiencing  a  change  in 
sequence  (123) .   Since  a  body  moving  in  space  generally 
contains  no  points  which  remain  stationary  between  two 
positions,  sequence  appears  to  be  an  invariable  property 
of  the  design  positions  on  a  given  branch  of  the  sixth- 
order  curve. 

6. 9   The  Objective  Function 

As  before,  the  objective  function  must  include  a 
term  for  each  design  requirement.   Thus  the  objective 
function  (O.F.)  for  the  three  precision  position  design 
of  the  RSSR-SC   and  RSSR-SS  mechanisms  may  be  expressed 
in  symbolic  form  as 

O.F.  =  B.A.T.  +  C.R.T.  +  T.C.T. 
+  F.P.L.T.  +  L.L.R.T. 
where 

B.A.T.  =  Branch-Avoidance  Term 

C.R.T.  =  Crank-Rotatability  Term 

T.C.T.  =  Transmission  Characteristic  Term 
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F.P.L.T.  =  Fixed-Pivot  Location  Term 

L.L.R.T.  =  Link-Length  Ratio  Term 

The  last  two  of  these  terms,  F.P.L.T.  and  L.L.R.T.,  are 
handled  in  the  same  way  as  the  similar  conditions  developed 
for  the  RCCC  mechanism  in  Chapter  Five.   Additional  motion 
requirements  could  be  included  in  equation  6.63,  as  dis- 
cussed in  section  6.7.   Notice  that  this  would  necessitate 
including  an  order  of  positions  term  in  equation  6.63. 

Transmission  Characteristic  Term.   The  third  term  of 
equation  6.63,  the  transmission  characteristic  term,  can  be 
taken  to  be  the  negative  of  the  minimum  value  of  the 
transmission  ratio  found  in  section  6.5,  i.e., 

T.C.T.  =  - (TR   .  .    ) 
minimum 

The  negative  sign  is  necessary  because  the  best  transmission 
corresponds  to  the  largest  value  of  TR.   The  first  two  terms 
of  equation  6.6  3  are  constraint  conditions  which  must  be 
handled  using  penalty  function  methods,  as  discussed  below. 

Branch-Avoidance  Term.   Equations  6.1  and  6 . 2  of 
section  6.3  give  a  nonparametric  means  for  determining 
whether  or  not  all  the  prescribed  positions  of  the  RSSR 
mechanism  lie  on  one  branch  only.   However,  these  conditions 
do  not  directly  provide  a  measure  of  the  amount  by  which 
a  given  mechanism  fails  to  meet  the  branching  condition. 
One  such  measure  is  now  proposed. 
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First,  define  a  new  variable,  BR.,  equal  to  the  left- 
hand  sides  of  equations  6.1  and  6.2  (which  are  the  same), 

BRj  =  Vj    *    ~2  x  ¥j  <6-64> 

Now  if  branching  occurs,  some  of  the  BR.  will  be  greater 

than  zero  and  some  will  be  less  than  zero.   In  this  case 

the  objective  should  be  to  either  force  the  negative  values 

in  the  positive  direction  or  force  the  positive  values  in 

the  negative  direction.   The  decision  of  which  direction 

to  move  is  based  on  which  set  has  the  largest  sum  magnitude. 

The  branch-avoidance  term  thus  becomes 

n 
B.A.T.  =  r{the  lesser  of  ( |  £   negative  BR . I , 

j  =  l  3 

n 

I   positive  BR . ) |     (6.65) 
j=l  3 

where  r  is  the  penalty  parameter.   It  can  be  seen  from 

the  above  expression  that  when  all  the  BR.'s  are  of  the 

J 

same  sign,  B.A.T.  will  be  equal  to  zero. 

Crank-Rotatability  Term.   Equations  6.4  3  and  6.44  give 
conditions  which  assure  crank  rotatability  in  the  RSSR 
mechanism.   Once  again,  however,  it  will  be  necessary  to 
have  some  measure  of  the  amount  by  which  these  conditions 
are  violated.   Using  the  logic  symbols  A  and  V  which 
respectively  denote  "and"  and  "or",  the  following  equations 
satisfy  the  above-mentioned  requirement: 
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C.R.T.  =  rA  if  A  <  0  A  (q  >  0  V  L  <  0)  (6.66) 
C.R.T.  =  r{the  lesser  of  (|q|,  |L|)}   if 

A  >  0  A  (q  <_   0  v  L  >  0)  (6.67) 

C.R.T.  =  r{A  +  the  lesser  of  (|q|,  |  L  | )  }   if 

A  <_   0  A  (q  <  0  V  L  >  0)  (6.68) 

where,  again,  r  is  the  penalty  parameter. 

The  effect  of  these  equations  is  to  assess  one  penalty 
for  violation  of  equation  6.43,  and  a  second  penalty  for 
violation  of  equation  6.44.   Since  only  one  of  the  two 
expressions  of  equation  6.44  must  be  satisfied,  the 
magnitude  of  the  penalty  is  proportional  to  the  lesser  of 
the  amounts  by  which  the  two  terms  are  in  violation. 

A  special  problem  arises  in  the  application  of  the 
crank-rotatability  conditions  which  merits  additional 
discussion.   When  conditions   6.43   and   6.44   are 
satisfied,  no  real  limit  positions  will  exist  for  the  input 
link  of  the  RSSR  mechanism.   This  situation  occurs  when 
either  the  input  link  of  the  mechanism  is  a  crank,  or  when 
the  mechanism  cannot  be  assembled  at  all.   The  precision 
position  synthesis  procedure  removes  the  latter  possibility 
because  it  guarantees  the  mechanism  can  be  assembled  at 
the  precision  positions.   However,  this  does  prevent  the 
solution  from  moving  in  this  direction.   In  other  words, 
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C.R.T.,  given  by  equations  6.66,  6.67  and  6.68,  has 
a  minimum  value  of  zero  when  the  input  link  is  fully 
rotatable  and  when  the  mechanism  cannot  be  assembled. 
Somewhere  between  these  possibilities,  it  has  a  maximum 
value  greater  than  zero.   The  optimization  process  may 
drive  the  solution  toward  the  unassemblable  minimum,  even 
though  it  can  never  actually  reach  a  value  of  zero. 
Because  of  this,  it  may  be  necessary  to  try  several 
different  starting  points  before  an  acceptable  solution 
is  found. 

6. 10   Numerical  Example 

The  following  positions  and  orientations  of  the 
moving  body  are  specified: 


o1  =  {Ox  +  Oj  +  Ok} 
o0  =  {li*  +  1]  +  l£}T 


o3  =  {li  +  2j  +  3k} 


.868 

-.494 

.045 

.489 

.834 

-.255 

.089 

.243 

.966 

.804 

-.590 

.071 

.577 

.745 

-.335 

.145 

.310 

.937 

Again  recall  that  [R  ]  is  defined  to  be  a  3x3  identity 
matrix  without  loss  of  generality. 
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Problem  Statement.   Design  an  RSSR-SR  mechanism 
which  satisfies  the  above  precision  positions  while  having 
acceptable  transmission  qualities.   The  solution  mechanism 
is  required  to  have  a  fully  rotating  input  link  and  must 
be  free  from  branching  problems.   Furthermore,  the 
fixed  pivots  of  the  mechanism  should  be  located  within 
ten  units  of  the  global  reference  origin,  and  the  ratio 
of  the  longest  to  the  shortest  link  length  should  not 
exceed  ten. 

Solution.   Solution  mechanisms  are  obtained  by 
minimizing  the  objective  function  of  equation  6.6  3  using 
the  Hooke  and  Jeeves  optimization  method.   Appendix  2 
provides  a  complete  listing  of  all  the  APL-language 
computer  programs  used  in  solving  this  problem  along 
with  a  sample  run  of  the  program  package. 

Figure  6.7  shows  the  RSSR-SR  mechanism  and  associated 
design  parameters  in  the  initial  position.   Again  it 
should  be  noted  that,  without  loss  of  generality,  the 
fixed  and  moving  coordinate  systems  have  been  assumed 
to  be  coincident  in  the  initial  position.   Tables  6.1  and 
6.2  give  the  results  of  two  separate  program  runs  which 
converged  to  feasible  solutions. 
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figure  6.7   The  RSSR-SR  mechanism  and  associated  design 
parameters 
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Table  6.1 
RSSR-SR  Mechanism  Optimization  -  Example  1 

Assumed  quantities  (refer  to  figure  6.7) 

?a  =  {8l  ~  33  +  3k> 
rb  =  {7i  +  53  -  2k} 
rc  =  {8i  -  6:  +  3k} 


After  nine  iterations  of  the  APL  program  HOOKE  and  a 
total  of  320  evaluations  of  the  objective  function,  all 
the  constraints  were  satisfied  and  the  solution  converged 
to  the  results  given  below. 


J  k 


ra  =  a1   =      8.355      -1.520      -1.400 
rb  =  b1   =      6.200       2.080        .200 


rQ  =  c±    =      7.200      -6.538       -.200 


aQ  =   6.546        .125       4.994 


bQ  =   5.335        .761       6.605 


cQ   =   6.850      -2.605       4.357 


-.173        .309 

.135        .160 

-.451        .449 


with  a  minimum  value  of  the  transmission  ratio 
TR  =  .663 


!i  = 

.935 

?2    = 

.978 

§3    = 

.771 
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Table  6.2 
RSSR-SR  Mechanism  Optimization  -  Example  2 

Assumed  quantities  (refer  to  figure  6.7) 
r   =  {3i  +  2j  -  2k} 

■v  a.         ~        ^,        ^ 

rb  =  {5i  -  43  +  2k} 

~c  =  {4  ~  X3  "  48} 

After  ten  iterations  of  the  APL  program  HOOKE  and  a 
total  of  312  evaluations  of  the  objective  function,  all 
the  constraints  were  satisfied  and  the  solution  converged 
to  the  results  given  below. 

<■  ->  /\ 

1  D  k 

ra   =   a1    =      3.940  2.925  -1.173 

~b  =  *?1  =  5'985  -3.924  -1.136 

rQ  =  c1  =  6.116  -1.000  -4.000 

aQ  =  4.195  .394  5.340 

bQ  =  5.143  -1.212  3.415 

cQ  =  4.327  .243  2.286 

s±   =  .952  .295  .077 

s2  =  .824  -.406  .394 

s3  =  .937  -.176  .301 


with  a  minimum  value  of  the  transmission  ratio 
TR  =  .551 


CHAPTER  7 
CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 
7. 1   Conclusions 

This  dissertation  has  attempted  to  lay  a  groundwork 
for  efficient  mechanism  design  and  optimization  in  the 
form  of  a  general  philosophy  of  mechanism  optimization. 
The  two  most  important  points  of  this  philosophy  are 
(a)  the  need  to  express  design  requirements  in  the  form 
of  nonparametric  constraint  equations,  and  (b)  the  need 
to  reduce  the  number  of  design  parameters. 

Design  and  optimization  of  mechanisms  to  meet  all  of 
the  various  real-world  requirements  have  been  emphasized, 
i.e.,  a  new,  "wholistic"  approach  was  attempted.   This 
is  in  contrast  to  many  recent  works  dealing  with  mechanism 
optimization  which  consider  only  precision  position  synthe- 
sis and  perhaps  one  or  two  additional  design  requirements. 
The  role  of  precision  position  synthesis  as  a  parameter 
reduction  tool  was  discussed,  along  with  the  possibility 
of  using  other  design  requirements  for  parameter  reduction. 
The  formulation  of  nonparametric  constraint  equations 
representing  the  various  design  requirements  was  addressed 
in  a  number  of  places  throughout  the  text.   This  process 
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was  seen  to  be  highly  problem  dependent,  often  requiring 
a  great  deal  of  insight  into  the  special  characteristics 
of  a  particular  mechanism.   The  formulation  of  these 
nonparametric  constraints  is  often  the  most  critical  step 
in  the  mechanism  design  process. 

Development  of  the  aforementioned  philosophy  was 
based,  to  a  large  extent,  on  a  study  of  past  research 
work  and  on  a  careful  examination  of  current  optimization 
theory.   The  optimization  method  of  Hooke  and  Jeeves  was 
selected  for  working  the  examples  of  this  dissertation 
because,  among  other  things,  it  is  simple,  requires 
minimal  computer  storage,  can  solve  a  large  number  of 
problem  types  and  requires  no  derivatives  of  the  objective 
function. 

The  philosophy  developed  in  the  first  three  chapters 
of  this  dissertation  has  been  applied  in  Chapters  Four, 
Five  and  Six.   The  types  of  mechanisms  treated  in  the 
examples  of  these  chapters  are  necessarily  limited. 
Nevertheless,  this  work  is  believed  to  contain  the  most 
complete  theories  yet  developed  for  the  design  of  the 
RCCC,  the  RSSR-SC  and  the  RSSR-SS  spatial  mechanisms. 
More  importantly,  perhaps,  the  examples  serve  to  illus- 
trate the  previously-discussed  philosophy,  particularly 
the  concept  of  parameter  reduction  and  the  formulation 
of  nonparametric  constraints. 

Many  of  the  theories  developed  in  this  dissertation 
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can  be  used  directly  to  solve  related  problems.   For 
example,  in  Chapter  Six,  constraint  conditions  were 
developed  for  the  RSSR  mechanism  as  an  intermediate  step 
in  the  design  of  the  RSSR-SC  and  the  RSSR-SS  mechanisms. 
These  constraint  conditions  are,  however,  directly  appli- 
cable to  the  RSSR  function-generating  mechanism.   Also, 
the  branching  condition  for  the  spherical  RRRR  mechanism 
developed  in  Chapter  Five  applies  directly  to  that  mechanism 
as  well  as  to  the  equivalent  RCCC  spatial  mechanism. 
The  use  of  the  position  stacking  method  for  the 
partially-exact  partially-approximate  synthesis  of  spatial 
mechanisms  has  been  developed  and  applied.   This  novel 
approach  allows  the  designer  to  satisfy  up  to  the  maximum 
number  of  precision  positions  and,  at  the  same  time, 
approximate  any  number  of  additional  positions.   One  of 
the  primary  benefits  of  this  method  is  the  parameter 
reduction  which  accompanies  precision  position  synthesis. 

Many  of  the  penalty  function  expressions  developed 
in  Chapters  Five  and  Six  are  believed  to  be  unique  to 
this  work.   An  example  of  this  is  the  branch-avoidance 
term  for  the  RCCC  mechanism,  given  in  equation  5.24. 
This  expression  serves  as  a  measure  of  the  amount  by 
which  the  branching  constraint  is  violated. 

7. 2   Recommendations  for  Future  Research 

It  is  believed  that  the  material  in  this  dissertation 
represents  a  significant  step  forward  in  the  theory  and 
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the  practice  of  mechanism  optimization.   Yet  it  must  be 
recognized  that  all  the  research  to  date  in  this  area 
represents  only  a  small  fraction  of  what  can  potentially 
be  accomplished.   This  section  should  by  no  means  be 
considered  a  complete  listing  of  the  unsolved  mechanism 
optimization  problems.   Rather,  it  is  a  brief  summary  of 
a  few  of  the  more  obvious  problems  which  would  follow  as 
an  extension  and  an  improvement  of  this  work. 

One  of  the  most  important  and  most  difficult  topics 
which,  to  a  large  extent,  remains  to  be  studied,  is  mechanism 
dynamics.   This  area  includes  balancing  and  kineto-elasto- 
dynamics  as  well  as  force  and  torque  analysis.   Improved 
theories  are  needed  for  both  planar  and  spatial  mechanisms 
which  will  allow  the  nonparametric  prediction  of  critical 
dynamic  properties.   At  present,  dynamic  analysis  generally 
employs  a  position,  velocity  and  acceleration  analysis  for 
closely-spaced  values  of  the  input  parameter.   This  type 
of  parametric  approach  is  very  inefficient  in  an  optimi- 
zation loop. 

Another  important  consideration  in  the  design  of 
mechanisms  is  the  effect  of  tolerance  and  clearance  on 
the  specified  performance.   No  precision  position  can 
truly  be  precise,  and  no  mechanism  dimension  can  be 
totally  accurate.   Obviously,  the  effect  of  these  on 
mechanism  performance  must  be  held  within  a  certain  range 
if  the  mechanism  is  to  be  acceptable.   But  what  is  this 
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range  and  how  can  it  be  determined  nonparametrically? 
Notice  that  when  tolerance  or  clearance  is  present,  all 
of  the  design  requirements  contain  some  degree  of  uncer- 
tainty.  This  points  towards  the  exploration  of  a 
statistical  or  stochastic  approach. 

Theories  for  the  design  and  optimization  of  spatial 
mechanisms  are  in  their  infancy.   It  is  not  yet  clear 
which  spatial  mechanisms  are  useful  and  for  what  purposes. 
The  questions  about  the  transmission  characteristics  of 
the  RCCC  mechanism  leave  some  doubt  about  its  practical 
usefulness.   Although  single-degree-of-freedom  seven-link 
spatial  mechanisms  are  well  known  (116)  ,  design  procedures 
are  available  for  only  the  simplest  three-,  four-  and  five- 
link  mechanisms.   Even  these  theories  are  quite  limited  in 
their  scope  of  application,  and  often  only  apply  to  special 
cases  of  these  mechanisms. 

To  date,  very  little  consideration  has  been  given  to 
mechanism  workspace  requirements.   Since  these  mechanisms 
will  inevitably  be  part  of  a  machine,  interference  with 
surrounding  components  must  be  avoided.   This  is  likely 
to  be  a  much  greater  problem  in  spatial  mechanism  design 
than  in  planar  mechanism  design  since  the  motions  are 
usually  much  more  difficult  to  visualize  and  to  represent 
on  graphic  computer  displays. 

The  author  realizes  that  many  of  the  problems  that 
have  been  so  casually  mentioned  here  will  be  extremely 
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difficult  to  solve.   And,  of  course,  countless  other 
important  problems  exist  which  have  not  been  mentioned 
here.   It  is  doubtful  that  the  goal  of  truly  complete 
mechanism  design  theories  will  be  attained  in  the 
foreseeable  future.   There  will  always  be  higher-order 
properties  and  unique  situations  to  challenge  the 
ingenuity  of  the  designer.   Still,  the  demands  of 
competitive  industry  and  of  an  increasingly  mechanized 
world  make  improved  design  theories  an  urgent  necessity. 
The  author  hopes  this  dissertation  has  made  at  least  a 
small  contribution  toward  the  goal  of  improved  mechanism 
design  theories. 


APPENDIX  1 

RCCC  MECHANISM  OPTIMIZATION  PROGRAMS  AND  SAMPLE  RUN 

This  appendix  contains  a  complete  listing  of  the  APL- 
language  computer  programs  required  for  optimizing  the 
RCCC  mechanism  as  discussed  in  Chapter  Five.   A  sample 
computer  output  of  a  typical  run  of  the  program  package 
follows  the  program  listings.   A  brief  description  of 
each  program  is  given  below. 


HOOKE     This  program  seeks  the  minimum  value  of  the 

function  FUN  given  an  arbitrary  starting  point. 
It  is  based  on  the  Hooke  and  Jeeve ' s  optimization 
procedure  described  in  section  2.4.4. 

FUN      This  program  calculates  the  objective  function  as 
given  by  equation  5.20. 

UNITVEC  This  program  calculates  a  unit  vector  given  two 
angles  (in  degrees) .  These  two  angles  serve  as 
the  design  variables. 

MET      This  program  calculates  the  motion  error  term 
(M.E.T.)  of  equation  5.20. 

RCCCSYN   This  program  is  called  by  MET.   It  calculates 
the  dimensions  of  an  RCCC  mechanism  given  the 
two  fixed-axis  directions  and  a  pair  of  scalar 
displacements  (see  Chapter  Four) . 

CCDYAD    This  program  is  called  by  RCCCSYN.   It  calculates 
the  dimensions  of  an  RC  or  CC  dyad  for  three 
precision  positions  of  rigid-body  guidance. 

AXIS      This  program  is  called  by  CCDYAD.   It  calculates 
the  moving  dyad  joint  axis  direction  given  the 
fixed  axis  direction. 
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CROSS 


This  program  calculates  the  cross-product  of 
two  vectors. 


CRT 


BAT 


OPT 


FPLT 


LLRT 


FBCALC 


EULER 


PRINT 


This  program  calculates  the  crank-rotatability 
term  (C.R.T.)  of  equation  5.20. 

This  program  calculates  the  branch-avoidance 
term  (B.A.T.)  of  equation  5.20. 

This  program  calculates  the  order-of-positions 
term  (O.P.T.)  of  equation  5.20. 

This  program  calculates  the  fixed-pivot  location 
term  (F.P.L.T.)  of  equation  5.20. 

This  program  calculates  the  link-length  ratio 
term  (L.L.R.T.)  of  equation  5.20. 

This  program  is  called  by  BAT.   It  performs  the 
repetitive  calculation  of  equation  5.23. 

This  program  is  not  directly  used  in  the  optimi- 
zation package.   It  calculates  the  rotation 
matrix  given  a  set  of  three  Euler  angles  (see 
section  4.3). 

This  program  is  called  by  HOOKE.   It  prints 
output  information  at  prespecified  intervals 
(see  the  sample  program  run  which  follows) . 
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RCCC  Mechanism  Optimization  Programs 


V  HOOKE 

I]  'INPUT  STARTING  POINT' 
!2]    N-*-pX0«-n 

,3]  'INPUT  INITIAL  STEP  SIZE  VECTOR' 

:  4  ]    A-HD 

5]    'INPUT  MINIMUM  STEP  SIZE' 

6]    AMIN«-Q 

7]    'INPUT  PENALTY  PARAMETER  VALUE' 

8]    PEN«-D 

9]    'INPUT  MAXIMUM  NUMBER  OF  ITERATIONS' 

10]   MMAX«-Q 

II]  'INPUT    THE    NUMBER    OF    ITERATIONS    BETWEEN    PRINTINGS' 
12]       MPRINT-HD 

13]       AX«-10*"6 

14]       M«-0 

15]       FE«-0 

16]  L4:X«-X0 

17]  L6  :Q«-1  ,0  +  M«-M  +  l 

18]       -»-(L8,L13  ,Ll4,Ll4)[l  +  (0=MPRINT|M)+2xMMAX  =  M] 

19]  L13 :PRINT    X,A 

!20]  L8:X1«-X0 
!21]       Xl[Q>Xl[Q]+A[Q] 
!22]       -»-(L2  ,Ll  )[1  +  (FUN    X1)<FUN    XO  ] 

!23]  Ll:X0-<-Xl 

24]      A[Q>1  .2xACQ] 

25]      +L7 

26]  L2  :Xl[Q]«-Xl[Q]+2xA[Q>-A[Q] 
!27]       ->(L3  ,L1)  [1  +  (FUN    Xl)<FUN    XO  ] 

:28]  L3:A[Q]+-A[Q]*2 

!29]  L7  :-*(L8,L5)[l+N<Q+Q  +  l] 

30]  L5:*(L9.L11)[1+a/( |X0-X)<AX] 

31]  L9:X1«-(2*X0)-X 
!32]   -+(L4,L10)[1+(FUN  Xl)<FUN  XO  ] 

33]  L10  :X«-X0«-X1 

34]   -+L6 

!35]  L11:*(L6,L12)[1+a/( |A)<AMIN] 

36]  L14 : '*****ITERATION  LIMIT  REACHED*****' 

37]   -»L15 

38]  L12 : '*****SOLUTION  HAS  CONVERGED*****' 

.39]  L15 :PRINT  XO.A 
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Programs  (continued) 


V  0F«-FUN    SD 

[1]  FE+FE+l 

[2]  Sl+UNITVEC    2+SD 

[3]  S4-HJNITVEC    2  +  2<t>SD 

[4]  SS-<-Sl  ,S4,"2  +  SD 

[5]  MET-f-MET    SS 

[6]  AXES<-(6  +  Vl)  ,  (3tl7<t>Vl)  ,3  +  14<fVl 

[7]  CRT-<-CRT    AXES 

[8]  BAT>B_A.T    AXES 

C9]  0PT-<-0_P_X(6tAXES)  ,    3+AXES 

[10]  FPLT«-I£Ll(3t9<t>Vl)  ,3  +  23<t>Vl 

[11  ]  LLRT^LLM(3  +  9<t>Vl)  ,  (01+3  +  64>Vl)  ,  (  01 +3  +  2O0V1 )  ,3  + 

[12]  OF-<-MET  +  PENx(100xCRT)+BAT  +  OPT+(l  +  FPLT)+l  +  LLRT 

[13]  0F«-0F  +  FPLT[2]+LLRT[2] 


V  UO-KJNITVEC    A;UOX 

[1]  A-*-(ol*180)xA 

[2]  UOX«-(1*1  +  ((30A[1])*2)  +  (30A[2])*2)*.5 

[3]  U0-<-D0X,(U0Xx3OA[l])  ,U0Xx3OA[2] 


23<t>Vl 


V 


V  MET-«-MIl    SS 

[1]  R13^RA 

[2]  Vl+RCCCSYN    SS 

[3]  R13^-RB 

[4]  V2-«-RCCCSYN    SS 

[5]  V1P^(3  +  90V1)  ,  (3  +  3<J>Vl)  .  (3  +  6<|>Vl) 

[6]  V1P^V1P,(3  +  234>V1)  ,(3  +  17<J)Vl)  ,3t20<J>Vl 

[7]  V2P^(3+9$V2)  ,(3  +  34>V2)  ,3  +  6<*>V2 

[8]  V2P*-V2P,(3  +  2  3  4>V2)  ,  (  3  +  1  7<t>V2  )  ,  3  +  20<t>V2 

[9]  MET>+/((+/(6    3pVlP-V2P)*2)*.5)xl     100     1     1     100    1 


V  V-f-RCCCSYN    SS;S1RC;S12;S13;  S1CC 

[1]  SlRC-<-AXIS    3tSS 

[2]  S12-K>0    S13-<-0 

[4]  V«-CCDYAD(3tSS)  ,0,0,S1RC 

[5]  V-<-(3tSS)  ,S1RC,V 

[6]  S1COAXIS    3  +  3<}>SS 

[7]  S12«-l"f6<t>SSO    S13-<-l  +  7<t>SS 

[9]  V^V,(3  +  34>SS)  ,SlCC,CCDYAD(3t3ct)SS)  ,S12,S13,S1CC 


Programs  (continued)  ^ 


V  RAS-CCDYAD  S S ; U2 ; Ul ; U3 ; M ; UO ; C ; UC1 ; UC2 ; UC3 

Cl]  U0^3  +  SSO  Ul«-  3iSSO  S12+SS[4]0  S13+SS[5] 

[5]  U2^Rl2  +  .xU10  U3-»-R13  +  .xui 

[7]  UC1-HJ1  CROSS  DOO  UC2+U2  CROSS  UOO  UC3+U3  CROSS  UO 

[10]  M«-U0.(-U0)  ,0,0,((S?R12)  +  .xU0)  ,(-UC)  ,(-  +  /U0xU2)  ,0 

[11]  M^M,((4)Rl3)  +  .xU0),(-U0),0,(-  +  /U0xU3),Ul,(-Ul),0,0 

[12]  M«-M.U1.(-U2).  1,0,U1,(-U3) ,0,"1 

[13]  M+M.(((^Rl2)+.xUC2)-UCl) ,(UC1-DC2),0,0 

[14]  M-<-M,(((i?R13)  +  .x0C3)-UCl)  ,  (UC1-DC3  )  .0  ,0 

[15]  M-*-8  8pM 

[16]  C-<-(-  +  /U0xOl)  ,(S12-  +  /U0xO2)  ,(S13-  +  /U0x03) 

[17]  OC,(-+/Ulxoi) , ((+/D2xS12xU0)-+/U2xO2) 

[18]  Cx-C,((+/D3xS13xU0)-  +  /U3xO3) 

[19]  C«-C,((-  +  /U0xO2  CROSS  U2)+  +  /U0xOl  CROSS  Dl) 

[20]  C«-C,((-+/U0xO3  CROSS  U3)+  +  /U0xOl  CROSS  Ul) 

[21]  RAS«-Cl]M 


V  Dl+AXIS  U0;A;B;C;D;U1Z 
[1]     ftGIVEN  THE  1ST  AXIS  OF  AN  R- C  OR  C- C  DYAD 
[2]     n  AND  3  PRESCRIBED  POSITIONS,  CALC  THE  2ND  AXIS 
[3]     ^ROTATIONS  R12  AND  R13  MUST  BE  GLOBAL  VARS 

[4]  A«-(($R12)-(i3)o.  =  i3)  +  .xU0 

[5]  B+(($R13)-(i3)-.=i3)+.xU0 

[6]  C^((B[2]xA[3])-B[3]xA[2])f(B[l]xA[2])-B[2]xA[l] 

[7]  D^((B[l]xA[3])-B[3]xA[l])f(B[2]xA[l])-B[l]xA[2] 

[8]  UlZ-<-(l*(C*2)  +  (D*2)+l)*.5 

[9]  U1-<-(CxD1Z)  ,  (DxUlZ)  ,U1Z 


V    C«-A    CROSS    B 
[1]  C-f-le(Axl$B)-BxleA 

V 


7  CRT-f-CjLX    U 

[1]  ALPHA^"20  +  /Uxl4>[l  ]u-<-4    3pU 

[2]  ALPHA-f-(ALPHAxALPHA^O,5)  +  (  (ol  )  -  ALPHA)  xALPHA>0  .  5 

[3]  AALPHA-<-(l  +  ALPHA)  -  1  +  ALPHA[  AALPHA] 

[4]  AALPHA+10xAALPHA+AALPHA>0 

[5]  CRT-<-.01+(  (o.5)+l  +  ALPHA)-.5x  +  /ALPHA 

[6]  CRT«-(  CRTxCRT>0)+AALPHA 

[7]  CRT«-CRT  +  CRT>0 


Programs  (continued)  ^84 


V  BATUMI    U;FB;FBG;FBL 

[I]  S-HJ+-4    3pU 
[2]  FB-<-FBCALC    U 

[3]  S[2;>R12  +  .xu[2;]0    S  [3  ;  >R1  2  +  .  xu[3  ;  ] 

[5]  FB+-FB.FBCALC    S 

[6]  S[2;>RA+.xu[2;]0    S [ 3 ; >RA  + . xU[ 3 ; ] 

[8]  FB-<-FB,FBCALC    S 

C9]  S[2;>RB+.xU[2;]0    S[3 ; >RB  + . xu[3 ; ] 

[II]  FB«-FB,FBCALC    S 

[12]  BAl>(+/FBxFB>.01)L l+/FBxFB<.01 

[13]  BAT«-BAT  +  BAT*0 


V    OPT>.Q£X    S;A;A41  ;TH;0RD;  IF;B 

[I]  Slf-3+SO    S21«-3  +  3<t>SO    S41«-~3tS 

[4]  A«-(S1     CROSS    S21),S1     CROSS    R12+.xS21 

[5]  A«-A.(S1     CROSS    RA+.xS21),Sl     CROSS    RB  +  .xS21 

[6]  A-«-4    3pA 

[7]  A41+-S41     CROSS    SI 

[8]  IF"-A+IxA41     CROSS    SI 

[9]  TH«-((    20A+.XA41  )xiF^0)  +  ((O2)-"2OA+.xA41  )xiF<0 

[10]  ORD+ATH 

[II]  B+-1    234143213241423 
[12]  B«-6    4pB,l    2_4    3    13    4    2 

[13]  T«-(+/(6  4p(  l+ORDil)<|)ORD)=B)  i4 

[14]  -*(L1,L1,L2,L2,L3.L3)[T] 

[15]  Ll:OPT«-0 

[16]  +LA 

[17]  L2  :OPT«-(  |  (  TH[  ORDi  2  ]  )  -  TH[  ORDi  3  ]  ) 

[18]  OPH-OPTL I TH[ORDil]-TH[ORDi4] 

[19]  +LA 

[20]  L3:OPT«-( | TH[ ORDi 1 ]- TH[ ORDi 2 ] ) 

[21]  OPT«-OPTL  I  (TH[ORDi3])-TH[ORDi4] 

[22]  LA:OPT-«-OPT  +  OPT>0 


V    FPLT-<-F_EIj;    AB;T 
[1]  T«-((+/(2    3pAB)*2)*.5)-100 

[2]  FPLT«-(+/TxT>0)  .+/T  +  100 

V 


V  LLRT-*-LLRT    JL;LL 

[1]  LL«-LL[ALL«-+/4    3p|JL-3<|)JL] 

[2]  LL«-LL[4]rLL[l] 

[3]  LLRT«-(  (LL>10)xLL)  ,LL 
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Programs  (continued) 


V  T-«-FBCALC  S 
[1]    T^((S[3;]-S[2;])CR0SS  S C 4 ; ]- S [ 3 ;  ]  )  +  .  xs [ 3  ;  ] 


V  RJ+EULER    ANG;C;S 

[1]  ANG«-ANGxO*180 

[2]  02  0ANG 

[3]  S^lOANG 

[4]  RJ+((C[l]xc[3])-S[l]xs[3]xc[2]) 

[5]  RJ«-RJ.(  Ix(c[l]xs[3])  +  S[l]xc[3]xc[2])  ,S[l]xS[2] 

[6]  RJ^RJ,((SCl]xC[3])+C[l]xC[2]xS[3]) 

[7]  RJ-<-RJ,(("lxSCl]xS[3])  +  C[l]xC[2]xcC3])  ,(~lxc[l]xs[2]) 

C8]  RJ^(3  3)pRJ,(S[2]xs[3]),(SC2]xC[3]),C[2] 


V  PRINT  X 
CI] 

[2]  'ITERATION  NUMBER  ';M 

[3]  'CURRENT  DESIGN  VECTOR  ';6+X 

[4]  'CURRENT  STEP  VECTOR  ' ;~6+X 

[5]  'OBJECTIVE  FUNCTION  VALUE  ' ; FUN  6+X 

[6]  'CONSTRAINT  CONDITION  VALUES' 

[7]  'CRT=';CRT 

C8]  'BAT=';BAT 

[9]  'OPT=';OPT 

[10]  'FPLT=' ;ltFPLT 

[11]  'LLRT=' ;ltLLRT 

[12]  'TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS= ' ; FE 

[13]  " 
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RCCC  Mechanism  Optimization  Sample  Run 


01 

0  0  0 

02 

4  2  0 

03 

8  2  3 

04 

10  3  5 

R12 

.96598  " 

.25849 

.00759 

.25816 

.96219 

".08682 

.01513 

.08583 

. 99619 

RA 

.86805 

.49444 

.04494 

.48852 

.83449 

.25489 

.08852 

.24321 

.96593 

RB_ 
.82358   .56409   .05939 
.54848   .76533  ".33682 
.14454   .30998   .93969 


HOOKE 
INPUT  STARTING  POINT 
D: 

270  45  190  175  3  3 
INPUT  INITIAL  STEP  SIZE  VECTOR 
D: 

25  25  25  25  3  3 
INPUT  MINIMUM  STEP  SIZE 
D: 

.01 
INPUT  PENALTY  PARAMETER  VALUE 
D: 

1000 
INPUT  MAXIMUM  NUMBER  OF  ITERATIONS 
D: 

25 
INPUT  THE  NUMBER  OF  ITERATIONS  BETWEEN  PRINTINGS 

5 
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Sample  Run  (continued) 


ITERATION  NUMBER.  1 

CURRENT  DESIGN  VECTOR  270  45  190  175  3  3 

CURRENT  STEP  VECTOR  25  25  25  25  3  3 

OBJECTIVE  FUNCTION  VALUE  2.5769E6 

CONSTRAINT  CONDITION  VALUES 

CRT=11 . 923 

BAT  =  0 

OPT  =  0 

FPLT=1381.8 

LLRT=0 

TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=l 


ITERATION  NUMBER  5 

CURRENT  DESIGN  VECTOR  302.5  62.5  202.5  200  0  .6 

CURRENT  STEP  VECTOR  9  ~9  3.75  ~3.75  .45  "1.08* 

OBJECTIVE  FUNCTION  VALUE  148.01 

CONSTRAINT  CONDITION  VALUES 

CRT  =  0 

BAT  =  0 

OPT  =  0 

FPLT=0 

LLRT=0 

TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS= 96 


ITERATION  NUMBER  10 

CURRENT  DESIGN  VECTOR  304.75  63.625  202.03  200.47 

.02812  .6 
CURRENT  STEP  VECTOR  .675  ".675  .28125  ".28125   03375 

.03375 
OBJECTIVE  FUNCTION  VALUE  132.37 
CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
OPT  =  0 
FPLT=0 
LLRT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=2 1 3 
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Sample  Run  (continued) 


ITERATION  NUMBER  15 

CURRENT  DESIGN  VECTOR  303.91  66.506  201.58  201  98 
.05017  .74404 

CURRENT  STEP  VECTOR  .2916  .69984  ".05062  .2916  "  0145 

.03499 
OBJECTIVE  FUNCTION  VALUE  113.67 
CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
OPT  =  0 
FPLT=0 
LLRT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=3 1 9 


ITERATION  NUMBER  20 

CURRENT  DESIGN  VECTOR  300.28  68.13  201.25  202  36 

.06912  .74842_ 
CURRENT  STEP  VECTOR  ".72559  .30233  ".05248  .05248 

.00262  .00262 
OBJECTIVE  FUNCTION  VALUE  107.14 
CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
OPT  =  0 
FPLT=0 
LLRT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=42 8 

*****ITERATION  LIMIT  REACHED***** 

ITERATION  NUMBER  25 

CURRENT  DESIGN  VECTOR  300.01  68.111  201.19  202  55 

.080  94  .7649 
CURRENT  STEP  VECTOR  .7523  ".02267  .00944  ".02267 

.00272  .00113 
OBJECTIVE  FUNCTION  VALUE  105.3 
CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
OPT  =  0 
FPLT=0 
LLRT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATION S=53 9 


APPENDIX  2 

RSSR-SR  MECHANISM  OPTIMIZATION  PROGRAMS  AND  SAMPLE  RUN 

This  appendix  contains  a  complete  listing  of  the  APL- 
language  computer  programs  required  for  optimizing  the 
RSSR-SR  mechanism  as  discussed  in  Chapter  Six.   A  sample 
computer  output  of  a  typical  run  of  the  program  package 
follows  the  program  listings.   A  brief  description  of 
each  program  is  given  below.   It  should  be  noted  that 
only  the  HOOKE  program  is  the  same  as  was  used  in  the 
RCCC  mechanism  optimization.   Other  programs  with  the 
same  name,  for  example  CRT,  serve  the  same  purpose  as 
before;  however,  the  programs  are  not  the  same. 

HOOKE      This  program  seeks  the  minimum  value  of  the 

function  FUN  given  an  arbitrary  starting  point. 
It  is  based  on  the  Hooke  and  Jeeve's  optimiza- 
tion procedure  described  in  section  2.4.4. 

FUN        This  program  calculates  the  objective  function 
as  given  by  equation  6.63. 

RSSRSRSYN   This  program  calculates  the  fixed-pivot 

locations  and  axis  orientations  of  the  RSSR-SR 
mechanism  given  the  locations  of  the  three 
spheric  joints. 

BAT        This  program  calculates  the  branch-avoidance 
term  (B.A.T.)  of  equation  6.63. 


MECHDIM 


This  program  calculates  the  dimensions  of  the 
RSSR-SR  mechanism  in  the  coordinate  system  of 
figure  6.6. 
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CRT        This  program  calculates  the  crank-rotatabilitv 
term  (C.R.T.)  of  equation  6.63. 

LLRT       This  program  calculates  the  link-length  ratio 
term  (L.L.R.T.)  of  equation  6.63. 

TR         This  program  calculates  the  transmission 

characteristic  term  (T.C.T.)  of  equation  6.6  3. 

FPLT       This  program  calculates  the  fixed-pivot 
location  term  of  equation  6.63. 

RSDYAD     This  program  is  called  by  RSSRSRSYN.   It 

calculates  the  fixed  pivot  location  and  axis 
orientation  of  an  RS  dyad  given  the  initial 
spheric  joint  location. 

REALROOT    This  program  is  called  by  CRT.   It  determines 
whether  or  not  a  given  quartic  equation  has 
real  roots. 

CROSS      This  program  calculates  the  cross-product  of 
two  vectors . 

PRINT      This  program  is  called  by  HOOKE.   It  prints 
output  information  at  prespecified  intervals 
(see  the  sample  program  run  which  follows) . 

EULER      This  program  is  not  directly  used  in  the 
optimization  package.   It  calculates  the 
rotation  matrix  given  a  set  of  three  Euler 
angles  (see  section  4.3). 


The  following  programs  were  authored  by  Dr.  Gary  K. 
Matthew.  They  calculate  the  real  and  imaginary  roots  of 
a  quartic  equation  and  are  called  by  TR. 


QUART 

CUBIC 

QUAD 

CROOT 

CR 

ATAN 
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RSSR-SR  Mechanism  Optimization  Programs 


V  HOOKE 

[I]  'INPUT  STARTING  POINT' 
[2]    N+pXO+Q 

[3]    'INPUT  INITIAL  STEP  SIZE  VECTOR' 

C4]    A+D 

[5]    'INPUT  MINIMUM  STEP  SIZE' 

[6]    AMIN«-Q 

[7]     'INPUT  PENALTY  PARAMETER  VALUE' 

[8]    PEN-HH 

[9]    'INPUT  MAXIMUM  NUMBER  OF  ITERATIONS' 

[10]  MMAX-<-D 

[II]  'INPUT    THE    NUMBER    OF     ITERATIONS    BETWEEN    PRINTINGS' 
[12]  MPRINT-HD 

[13]  AX-t-10*    6 

[14]  M+-0 

[15]  FE+-0 

[16]  L4:X-<-X0 

[17]  L6  :Q«-1  ,0  +  M-f-M+l 

[18]  ->-(L8,L13  ,Ll4,L14)[l  +  (0=MPRINT|M)+2xMMAX  =  M] 

[19]  L13:PRINT    X.A 

[20]  L8:X1«-X0 

[21]  X1[Q]-X1[Q]+A[Q] 

[22]  *(L2,L1)[1+(FUN    X1)<FUN    XO ] 

[23]  LI :X0+X1 

[24]  A[Q]«-1.2xA[Q] 

[25]  +L7 

[26]  L2:Xl[Q>Xl[Q]+2xA[Q]-«-A[Q] 

[27]  -*(L3,Ll)[l  +  (FUN    X1)<FUN    XO  ] 

[28]  L3:A[Q>A[Q]*2 

[29]  L7  :  +  (L8,L5)[l+N<Q-<-Q+l] 

[30]  L5:+(L9,L11)[1+a/( |X0-X)<AX] 

[31]  L9:Xl«-(2xxO)-X 

[32]  f(L4,L10)[l+(FUN  X1)<FUN  XO ] 

[33]  L10:X-^X0-<-Xl 

[34]  -+L6 

[35]  Lll :+(L6,L12)[1+a/( |A)<AMIN] 

[36]  L14 : '*****ITERATION  LIMIT  REACHED*****' 

[37]  +L15 

[38]  L12 : ******SOLUTION  HAS  CONVERGED*****' 

[39]  L15 :PRINT  XO.A 
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Programs  (continued) 


7  OF«-FUN    RRR 

[I]  FE«-FE  +  1 

[2]  SA«-RSSRSRSYN    RRR 

[3]  BATl«-fiAX(3t3<J>RRR),3  +  18<t>SA 

[4]  DIM1«-MECHDIM    18+SA 

[5]  CRTl^CRT    DIM1 

C6]  TR1«-£R.    DIM1 

C7]  LLRTlfLLEI    DIM1 

[8]  BAT2-B_Al(3  +  6<|>RRR)  ,  (  3  +  9+  SA)  ,  "9+  SA 

[9]  DIM2-»-MECHDIM(9tSA)  ,~9+SA 

[10]  CRT2-*-CRI  DIM2 

[II]  TR2-«-XE    DIM2 
[12]  LLRT2^LLRT    DIM2 
[13]  FPLT>FPLT    SA 

[14]  OF^BATl+BAT2+10xCRTl+CRT2 

[15]  0F+(-l+(TRl ,TR2)[iTRl,TR2])+OF 

[16]  0F+0F+LLRT1+LLRT2+FPLT 


V  SA«-RSSRSRSYN    RRR 

[1]  SA-*-RSDYAD    3  + RRR 

[2]  SA-<-SA,RSDYAD    3  +  3<t>RRR 

[3]  SA+SA.RSDYAD    "3  +  RRR 

7 


7    BAT«-B_AX    X;S2;RB;AB;B0;C;T1  ;T2 

[I]  S2-K3t9<t>XO    RB-«-3  +  X 

[3]  AB«-(    3  +  X)-3  +  6<f>XO    B(K3t  12<|>X 

[5]    OAB+.XS2  CROSS(0[l; ]+RB)-B0 

[6]    C^C,(R2+.xAB)+.xS2  CROSS ( 0[ 2 ; ] +R2 + . xRB ) -BO 

[7!    C^-C,(R3  +  .xAB)  +  .xS2  CROS  S  (  0[  3  ;  ]  +  R3  +  .  xRB  )  -BO 

[8]    Tl-«-+/CxC£0 

[9]    T2«-| +/CXCS0 

[10]   +(LA,LB)[l+(Tl=0)vT2=0] 

[II]  LA:BAT-hT1  LT2 
[12]   BAT«-BAT  +  BAT>0 
[13]   +0 

[14]  LB:BAT-<-0 
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Programs  (continued) 

7  DIM^MECHDIM    Al  A2 ; A; B ; C ; E ; F ; G ; PS  I ; SI ; S2  ;  AO  ;  Al  ;  BO  :  Bl  • 
T1;T2 

Cl]  S1+3  +  A1A20   A0+3  +  3<J>AlA2O   Al«-3  +  64>A1  A2 

[4]  S2-3+94>AlA20   BO+3  +  12  <J>A1  A20   B  1+3  +  1  5  <t>Al  A2 

[7]  A«-((Al-A0)  +  .xAl-A0)*.5 

[8]  O((B1-B0)+.xb1-B0)*.5 

[9]  B-«-((Al-Bl)  +  .xAl-Bl)*.5 

[10]  Tl-*-Sl     CROSS    S2 

[11]  E«-|  (A0-B0)  +  .xTli(+/Tl*2)*.5 

[12]  F-<-((A0-B0)  +  .xS2    CROSS    Tl)*Tl+.xTl 

[13]  G«--((A0-B0)  +  .xsi     CROSS-Tl)fTl  +  .xTl 

[14]  PSI-*-~20Sl  +  .xS2 

[15]  PSID-^PSIx180tO1 

[16]  DIM-«-A,B,C,E,F,G,PSI 


V  CRT-£M    D;T1;SK1;SK2;SK3;SK4;SK5;SK6;K1;K2;K3;K4-K5 
;B1;C1;D1;E1;SQ;SR;SS;A;B;C;E;F;G;PSI 

[1]  A+D[1]0    B-<-D[2]0    C«-D[3]0    E-«-D[4] 

[5]  F-<-D[5]0   G«-D[6]0    PSI«-D[7] 

[8]  SK1«-(F*2)  +  (A*2)  +  (E*2)  +  (G*2)  +  (C*2)-B*2 

[9]  SKl-^(FxGx(20PSI))-.5xSKl 

[10]  SK2«-AxEO    SK3-<-AxGxlopSI 

[12]  SK4-<-CxEO    SK5t-CxA 

[14]  SK6-<-CxFx1oPSIO    SK7-<— CxAx20PSI 

[16]  K1+-(SK4*2)  +  (SK6*2)  +  (SK7*2)-(SK1*2)+SK3*2 

[17]  K2^(2xSK4xSK5)-2xSKlxSK2 

[18]  K3^(2xSK6xSK7)-2xSKlxSK3 

[19]  K4-<-  2xSK2xSK3 

[20]  K5"-(SK3*2)  +  (SK5*2)-(SK2*2)  +  SK7*2 

[21]  Bl«-((2xK2xK5)+2xK3xK4)*Tl-<-(K4*2)+K5*2 

[22]  C1-((K2*2)+(K3*2)+(2xk1xK5)-k4*2)tT1 

[23]  Dl  +  (  (2xr1xK2)-2xK3xK4)tT1 

[24]  E1^((K1*2)-K3*2)*T1 

[25]  CRT>REALROOT    1  ,  B  1  ,  Cl  ,  Dl  ,  El 

[26]  CRT-<-CRT+CRT>0 

[27]  TEMP-t-l.Bl.Cl.Dl.El 


V  LLRT^LLRJT  DIM 
[1]    LLRT-«-(4+DIM)[MtDIM] 
[2]    LLRT-<-LLRT[4]tLLRT[1] 
[3]    LLRT-<-LLRTxLLRT>10 

V 
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Programs  (continued) 

V  TH-2B,    D;A;B;C;E;F;G;PSI;CA;CB;CC;TR1;TR2;TR3 

[1]  A*-D[1]0    B*D[2]0    C-<-D[3]0    E«-D[4] 

[5]  F«-D[5]0    G«-D[6]0    PSI+-D[7] 

[8]  R«-(A*2)  +  (C*2)  +  (E*2)  +  (F*2)+G*2 

[9]  R-<-R-  (B*2  )+2xFxGx20PSI 

[10]  CA^(10PSI)x(RxG)-2x(c*2)xFx20PSI 

Cll]  CB-<-ExR-2xC*2 

[12]  CO(2*(02)x(20PSl)*2)+2xE*2 

[13]  CC^.5xAxCC-(2x(c*2))+2x(G*2)x(loPSI)*2 

[14]  CD-«-2xAxExGxlopsi 

[15]  T^(CD-CA),(2xCB-2xCC) ,("6xCD) ,(2xCB+2xCC),CA+CD 

[16]  TH-hQDART    T 

[17]  TH«-(  ,14>TH  =  lp0)/,TH 

[18]  TH-«-2x    30TH 

[19]  TRl-<-(R+(2xAxEx20TH)-2xAxGx(loPSI)x(10TH))*2 

[20]  TR2^4x(c*2)x(e+Ax20TH)*2 

[21]  TR3-^4x(C*2)x(  (Fx10PSI)-Ax(2  0PSI)x1oTH)*2 

[22]  TR-«-(TR3+TR2-TRl)*4x(B*2)xc*2 

[23]  TI*-((TR>0)/TR)*.5 


7  FPLT-t-FP^T    SA 

[1]  A0-»-3  +  34>SAO    B0«-3  +  12<t>SAO    CO+3  +  2  1  <J>SA 

[4]  FP-f-((+/A0*2)  ,  (+/B0*2)  , +/C0*2)*.5 

[5]  FPLT«-+/FPxFP>5 


V  Z-RSDYAD    X;A;P2;P3;C1;C2;S;A0;L3;RJ 

[I]  R«-X 

[2]  RJ«-(3     3)p(Rl+.xR) ,(R2+.xR) ,R3+.xr 

[3]  A-^OJ  +  RJ 

[4]  S«-(A[2;]-A[1;])CR0SS    A[3;]-A[l;] 

[5]  S«-S*(+/S*2)*.5 

[6]  P2-HS    CROSS    P2*(+/(P2«-A[2;]-A[1;])*2)*.5 

[7]  P3«-S  CROSS  P3*(+/(P3«-A[3;]-A[2;])*2)*.5 

[8]  C1«-(P2  CROSS  A[l;])+P2  CROS S ( A[ 2 ; ]- A[ 1 ; ] ) *2 

[9]  C2*(P2  CROSS  A[2;])+P2  CROS S ( A[ 3 ; ] - A[ 2 ; ] ) *2 

[10]  L3-«-(L3  +  .xCl-C2)iL3  +  .xL3-«-P2  CROSS  P3 

[II]  A0^(L3xp3)+(A[3;]+A[2;])^2 
[12]  Z«-S,A0,A[1;] 
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Programs  (continued) 


7  A«-REALR00T  X;Q;R;S;P 

[I]  Q^(("3*X[2]*2)t8)+X[3] 

[2]    R^X[4]+((X[2]*3)*8)-(X[2]xX[3])*2 

[3]    S+X[5]+(-(X[4]xx[2])*4)+((X[3]xx[2]*2)*16) 

[4]    S^S-(3xxC2]*4)t256 

[5]    P«— (4xS)+(Q*2)*3 

[6]    CQ-«-((8xQxs)t3)-(R*2)  +  (2xQ*3)t2  7 

[7]    L^(8xQxS)-(2xQ*3)+9xr*2 

[8]    A-*— (4xp*3  )+27xCQ*2 

[9]    T«-(a>0)a(q;>0)vl<;0 

[10]   -KLA,LB)[T+1] 

[II]  LA:A-<-+/((A£0)  ,(Q<0)vl>0)x  |a,(  |Qxn<o)LLxL>0 
[12]   -0 

[13]  LB:A-<-0 


V    C-«-A    CROSS    B 
[1]  01e(Axl<}>B)-BxleA 


V  PRINT  X 
[1] 

[2]  'ITERATION  NUMBER  ';M 

[3]  'CURRENT  DESIGN  VECTOR  ';9  +  X 

[4]  'CURRENT  STEP  VECTOR  ';"9+X 

[5]  'OBJECTIVE  FUNCTION  VALUE  • ; FUN  9+X 

[6]  ,TR=,;lt(TRl,TR2)UTRl,TR2] 

[7]  'CONSTRAINT  CONDITION  VALUES' 

[8]  'CRT=' ; CRT1+CRT2 

[9]  'BAT=» ;BAT1+BAT2 

[10]  'LLRT= ' ;LLRT1+LLRT2 

[11]  'FPLT=';FPLT 

[12]  'TOTAL  NUMBER  OF  FUNCTIONAL  EVALUAT IONS= » • FE 
[13]" 


V  RJ-<-EULER  ANG;C;S 

[1]  ANG-<-ANGxo*180 

[2]  C«-20ANG 

[3]  S-<-10ANG 

[4]  RJ*((C[l]xc[3])-S[l]xs[3]xc[2]) 

[5]  RJ«-RJ,(  lx(C[l]xS[3])+S[l]xC[3]xC[2])  ,S[llxS[2] 

[6]  RJ^RJ,((S[l]xC[3])+C[l]xc[2]xS[3]) 

[7]  RJ+RJ,((  Ixs[l]xs[3])+C[l]xc[2]xc[3]),("lxc[l]xs[2]) 

[8]  RJ^(3  3)pRJ,(S[2]xS[3]),(S[2]xC[3])  ,C[2] 


[1] 

[2] 

C3] 

C4] 

[5] 

[6] 

[7] 

[8] 

C9] 

[10] 

[11] 

[12] 

[13] 

[14] 

[15] 
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Programs  (continued) 

7    Z«-QUART    X;F;A;D2;D;V;DI0 
QlO-*-l 

X«-Xx1*1+x«-1+XtX[1  ] 
A«-.25*X[l] 

+(<(1+X[4]) ,(0A.=x|x[l     3]),1+X[1])=1+F+0)/STOP,OUT,G 
0 

X«-0,(Al6    3    lxl.xCl    2]),(Al4    3    2    lxl, x[l    2    3]),Axl,X 

*((1+X[4    3])=F^1)/ST0P,0UT 
GO:Z-hCUBIC    1  ,  (  2  xX[  2  ]  )  ,  ( -/X[2    4] xX[ 2 ] , 4) , -X[ 3 ]*2 

D+(D2«-lt(Z[;l]>0)/(l=l+Z[;2])/Z[;l])*.5 

Z«-(QUAD    l,D,.5xX[2]  +  D2-X[3]*D)  ,[1]QUAD    1  ,  (  -  D)  ,  .  5  xX[  2 

]+D2+X[3]*D 

-►OUT  +  2 
STOP:Z«-0,[1]CUBIC    1,X[1    2    3] 

-»-(F.l)/0UT+2    3 
0UT:Z<-2     CROOT    QUAD    1,X[2     4] 

-»-(F,l  )/0UT  +  2    3 

Z+Z+4    2pA,0 


7    Z-t-CUBIC    X;A;B;C;D;Q;R;T;0;DI0 
[1]         -»-9xil  =  (niO*l)+0lX 

[2]         X+U+XCl]) ,(B+X[2]+3) ,(OX[3]i3) ,D«-X[4] 
[3]         -►6x1(q  =  o)vO^T^((Q^(AxC)-Bxb)*3)  +  (R-<-(-B*3)  +  .5x(x/3,A, 

B,C)-x/A,A,D2*2 
[4]  Z«-(3     1p((2xx-.5  +  r^o)x(q*.5)x2o((^3)x-2o(  |  R)  *  (  Q«--Q)  *1 

.5)+0    1    2x02*3)), 0 
[5]         ->-7,,Z«-(Z-3    2pB,0)rA 

[6]  Z«-3    2p(+A)x(-B-+/C),0,Q.R,(Q*(-.5x  +  /c)-B),-R-f-.5x(3*. 

5)x-/C«-(R+T,-  T>T  *  .  5  )  *  t  3 
[7]  Z«-Zxl*l+z 

[8]         ->-0 
[9]  Z«-0,[1]QUAD    X[l     2    3] 


7  Z«-QUAD    X;DlO 

[1]  DI0*1 

[2]  +4xiO*Z«--([]CT<|Z)xZ«-(4xX[3])-(X«-XtX[1])[2]*2 

[3]  +0,,Z-(2    lp.5x(-x[2]-Z*.5).-X[2]+Z*.5),0 

[4]  Z«-2    2p.5xx,(-Z),(X-*--X[2])  ,Z«--(-Z)*.5 

7 
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Programs  (continued) 


V  Z+-FL    CROOT    F;I 

[1]  Z«-0    2pl+~l+ii 

C2]  Z-<-Z,[l]FL    CR,F[]>I  +  1  ;  ] 

[3]  +2XK1  +  PF 


V  Z«-N  £R  X;T 

Cl]    Z«-(((+/X*2)*.5)**N)x$(2,N)p(20T)  . 1 OT-( *N ) x ( ( o2 )  I  (1+X 
)ATAN  l+X)+02x  l+lN 
7 


7  Z^Y  ATAN  X;Xl;S;XO;IZ;QlO 

Cl]  DlO-«-l 

C2]  S-«-pX 

[3]  Xl-«-X+(X-«-,X)=0 

[4]  Z+-   3o(Y-<-,Y)fXl 

C5]  Z[X0/lZ-f-ipZ]-t-.5xox(X0^X  =  0)/Y 

[6]  Z[Xl]f-Z[Xl-HXO/IZ]+Ox(xO^-X<0)/Y 

[7]  Z[(X0AY=0)/IZ]*Ol 

[8]  Z-<-SpZ 
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RSSR-SR  Mechanism  Optimization  Sample  Run 


0 

0  0  0 

1  1  1 
1  2  3 

R2 

.86805  ".49444  .04494 

.48852   .83449  ".25489 

.08852   .24321  .96593 

R3 

.80393  ".59045  .07111 

.57688   .74517  ".33455 

.14454   .30998  .93969 


HOOKE 
INPUT  STARTING  POINT 
D: 

8  "3  375   28  "6  3 
INPUT  INITIAL  STEP  SIZE  VECTOR 
□  : 

111111111 
INPUT  MINIMUM  STEP  SIZE 
D: 

.05 
INPUT  PENALTY  PARAMETER  VALUE 
D: 

100 
INPUT  MAXIMUM  NUMBER  OF  ITERATIONS 
D: 

25 
INPUT  THE  NUMBER  OF  ITERATIONS  BETWEEN  PRINTINGS 
D: 

2 

ITERATION  NUMBER  1 

CURRENT  DESIGN  VECTOR  8  ~3  3  7  5  "2  8  "6  3 

CURRENT    STEP    VECTOR    111111111 

OBJECTIVE  FUNCTION  VALUE  71.056 

TR=. 54516 

CONSTRAINT  CONDITION  VALUES 

CRT=1 .3868 

BAT  =  0 

LLRT=26.087 

FPLT=31 .646 

TOTAL  NUMBER  OF  FUNCTIONAL  EVALUAT ION S = 1 
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Sample  Run  (continued) 


ITERATION  NUMBER  2 

CURRENT  DESIGN  VECTOR  6  "3  1  5  5  "2  6  "6  1 

CURRENT  STEP  VECTOR  "1.2  ".5  "1.2  "1.2  "  5  "  5  "1  2 
1.2  ... 

OBJECTIVE  FUNCTION  VALUE  ".03536 

TR=. 03536 

CONSTRAINT  CONDITION  VALUES 

CRT  =  0 

BAT  =  0 

LLRT=0 

FPLT=0 

TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=3 9 


ITERATION  NUMBER  4 

CURRENT  DESIGN  VECTOR  8.4  ".8  "1.4  6.2  2  8   2  7  2  "6  5 
.2  .   .    . 

CURRENT  STEP  VECTOR  ".72  .72  .72  .72  "72   72   72  "  3 
.72  .... 

OBJECTIVE  FUNCTION  VALUE  ".60477 

TR=. 60477 

CONSTRAINT  CONDITION  VALUES 

CRT  =  0 

BAT  =  0 

LLRT=0 

FPLT=0 

TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=10 2 


ITERATION  NUMBER  6 

CURRENT  DESIGN  VECTOR  8.4  "1.52  "1.4  6.2  2  08   2  7  2 

6.5   .2       _  *" 

CURRENT  STEP_VECTOR  ".18  .432  .18  .18  .432   18   18 

.075   .18 
OBJECTIVE  FUNCTION  VALUE  ".66116 
TR=. 66116 

CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
LLRT=0 
FPLT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATI0NS=1 7 5 
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Sample  Run  (continued) 


ITERATION  NUMBER  8 

CURRENT  DESIGN  VECTOR  8.4  "1.52  "1.4  6  2  2  08   2  7  2 

6.53  7  5  ".2 
CURRENT  STEP_VECTOR  ".045  .108  .045  .045  .108  .045   045 

".045   .045 
OBJECTIVE  FUNCTION  VALUE  ".66207 
TR=. 66207 

CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
LLRT=0 
FPLT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=25 0 

*****SOLUTION  HAS  CONVERGED***** 

ITERATION  NUMBER  9 

CURRENT  DESIGN  VECTOR  8.355  "1.52  "1.4  6.2  2  08   2  7  2 

6.5375  ".2 
CURRENT  STEP  VECTOR  .027  .027  .01125  .01125  .027  .01125 

.01125   .01125  ".01125 
OBJECTIVE  FUNCTION  VALUE  ".66329 
TR=. 66329 

CONSTRAINT  CONDITION  VALUES 
CRT  =  0 
BAT  =  0 
LLRT=0 
FPLT=0 
TOTAL  NUMBER  OF  FUNCTIONAL  EVALUATIONS=32 3 
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